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APPENDIX 1 


Listing of Program and 
Documentation of Use 


44 



USE 


0 F 


S P E C T 


THE FOlLakrlNr, IS ^ LIST OF INPUT PARAMETER'S AND DATA CARDS AND THEIR 
INTERPRET AT 1 JNS . 

JOOM\N-0 JMA I N :iF ANALYSIS* 0 OR rtL ANK = PE R I OD , INFREQUENCY. 

J BANDS- HW OF FREQUENCY OR_ PF R.X0ll_J^ AND S_XLU-_B.E-„OB.TAlit£CLi 

0 OH *^lVnk = GENRRAt'F AfciipniNG TO LIMITS SUBSEQUENTLY READ IN. 

_ l = HEAn IN A VFCTUR UF BANDS Tu BE USED FR!)M CAROS. , 

N&ANDS-M:1W M .^NY THFUJENCY DR PERIOD BANOS WILL BE USED. 

NOBS-THF <NUMrtER OF nnSERVATIONS IN THE DATA VECTOR. 

JANALY-H.nw The OATA SHALL BE ANALYZED. 

0 PH BLANK^nn UNIVAKIAT F SPECTRUMt SEARCH FOR' BEST K-VARlAT£ 

MODEL 

l = S IMPLY DO UNIVARIATF SPFCTRUM .. , 

2=Dfl K-VAPIATF SPECTRUM ACCORDING TO MODEL PARAMETER READ IN 

^JpATA-TElLS SilHOUTINE DF PVAR WHE^E TO GET DATA* 

0 OK BLAN:< = REA0 FROM CARDS ACCORDING TO FORMAT (IOX,7FIO.O) 

1=GFN E RA^E data F ROM MCDE LING PROGRAM CALLED BY DATA. THIS 

MODELING PROGRAM HILL ASK FOR MORE CAROS. 

jcrjtr-sflects cpitepion fop OPTIMIZATION of K- V AR I ATE__SPEC TRUM. 

0 OR BLANK=.MUL T I PLF R-SOUARED 

1 = T-VALUES . 

JPR INT-PR 1 NTOUT control DIGIT 

Q-DO N_r._P R I NT ANY RFSUI TS 

IMPRINT JNl Y FINAL RESULTS 

2=PRINT FINAL ANO INTEPMEOIATE RESULTS . . 

• 3=PRINT as FDR 2 PLUS MATRICES. CORREL AT I ONS , ETC FOR EACH STEP 
JPLOT-PLDT CUNTROL DIGIT, VALUES A5_AB.QVE EXCEPT .FOR,. 3. 

BANDLG-Ln^CS T_F R E Q . OR t ONGFST PERIOD. AS APPR O PRIATE TO DOMAIN, 
PAnT)EL-OFL TAI^reO.) UELTAI PFR IOD I » AS APPROPRIATE TO DOMAIN., 


PANDS-VfC'TOR OF 8ANOS- 

NS IN-NUM.afR of sine WAVES IN SIMULATED DATA. 

^^FR F Q , LM _PF R I O0_p F _TH^L_S I NTH SIMUL ATED SINE WAVF 

S 0 - S -"A N DA R 7) 11 E V i A f Yo N ■ 0 f" T H E L S I N T M SIMULATED SINE WAVE. (AMPLITUDE) 
SDN-STANDARU OFV. OF GAUSSIAN NOISE . „ 


i 


ORIGINAL PAGE IS 

OF POOR QUAIITS! 


THF FaiLUWlM'', A LIST OF INPUT CARDS AND THEIR CONTENTS 


CARD 1 

- 

- - 

MANTATURY 




roi 

UMN 

VAT I A!.U, F 

COMMENTS 

- - - 


1 

- 

5 

- JDUMAN 

R IGHT 

JUSTIFY 


6 

- 

10 

jmands 

R I GHT 

justify 


11 

- 

IS 

- NMANOS 

RIGHT 

JUSTIFY 


16 

- 

2 0 

- NOBS 

R IGHT 

JUSTIFY 


21 

- 

2S 

- JANALY 

R IGHT 

JUSTIFY 


26 

- 

30 

JTATA 

R IGHT 

JUSTIFY 


31 

- 

3S 

- JCRTTR 

R IGHT 

justify 


36 

- 

40 

- JPRINT 

RIGHT 

JUSTIFY 




4S 

- JPLOT 




CARO 2 

- 

- - 

RFDUIRED 

ONLY IF J0ANOS=O» OTHERWISE OMIT. 



COLUMN , VA^IAiiLE COMMENTS - - 

I - 10 - ;UNOLO punch DECIMAL 

I V 2 0. . - , _ 3 A N 0 E L PUN C_H „ D E C..T M .A L 


:ARDS 3 - - - RF Q UIREO ONLY IF JQANDS= 
COLUMN VARIArtLF COMMENTS - - - 


10 - RANDSIII . PUNCH DEC IMAL 1ST BAND IN SPECTRUM, 
11 - 20 - 3ANUS(2» PUNCH DECIMAL - 2ND BAND IN SPECTRUM. 




♦ 

CONTINUE for as many CARPS AS NEEDED TO READ IN NBANOS BAND' 
FORMAT! 3F 10, 0) 


Ca”rOS a - - - PEOtHREO ONLY IF J DAT A = 6, ' . 

. COLUMN VARIABLE COMMENTS - - “ ' 

I - 10 - ID rone may contain anything - not read. 

It - 20 - r> ATM t ) 1 ST DATA PO INT , PUNCH D E C I M AL O R R ! GH T_^ U 5X1 F^.____ ' 

21 - 30 - DATA! 2) 2ND DATA POINT, PUNCH DEC IMAL OR RIGHT JUSTIFY. 

■ 

* 

......CONTINUE SAME FORMAT ON CARpS...A,5 NE.EDE D..FO.R .NOBS D AT A POINTS. 

FORMAT ( lOX, 7FIO.O > i 


CAR05 5 - - - RFJUIkEO ONLY IF J0ATA = 1. 

5A - “ MANDATORY , 

COLUMN VARIABLE COMMENTS - - - 
I • 5 - NSIN RIGHT JUSTIFY, 


5B - - R .E . j ( I I.P.E ON LY__I_^NS.1 N._ t S NOT ZFR Q. ? 

COLUMN VAMABLE COMMENTS - - - 

I - to - RANDin PUNCH DECIMAL,. ... . ; 

ll - 2^ - SO in PUNCH DECIMAL. [ 

CONTINJE same format thru NSIN capos (BANOSI. FORMAT ( 2F 10,0) i 

VANpATt);^ Y ' ! 

COLUMN VA^IAiR,: COMMENTS - * - ■ 

1 - to - SilN PUNCH DECIMAL 

I 

SO - - kr JUT RED UNI. Y if SDN IS NOT 2FP0 | 

COLUMN VA'MAfALE COMMENTS - - - | 

^ lu C'ODF MAY CONTAIN ANYTHING - NOT RFAD, i 

- 10 N)I^rn) PUNCH tJCCIMAI, 1 St”ga‘usS I AN RANDOM NUMBER . ’ ; 

11-16 ^4 1TSE(2) . PUNCH OECIMAL. 2ND GAUSSIAN .RANDOM NUMBER. 

* 

' i 

continue Thj j ^S MANY CAROS AS NEEDED FOR NOBS RANDOM NUMBERS. { 

format ( .vx, I 6 . 4 4 


ORIGINAL PAGE IS 
OF POOR QUALlTYf 




A 


c „ _ 

c 

jC 

c 

0IH6NSI0N BANOSl 50KBAN( 50 ) t PEAKS < 50 1 ,PHASE ( 501 ,AMNAM£ ( 9 » *TNAH6 

U9I .P^NAME^9) ,BANAME ( 5) ,XNAME<5 I ,LVAft* 100 »♦ ISAVEUOO) 


THE" FOLLO'^ING OTMENSION' MUST BE EQUAL TO OR GREATER THAN 
NOaS*( NBANDS^2<-1 > 


SPT 

10 

SPT 

20 

SPT 

30 

SPT 

^0 

SPT 

50 

SPT 

60 

SPT 

70 

SPT 

80 

SPT 

90 

SPT 

100 

SPT 

110 

SPT 

uo 

SPT 

130 

SPT 

“l60 

SPT 

150 

SPT 

160 

SPT 

170 

SPT 

160 

SPT 

190 

SPT 

200 

SPT 

210 

SPT 

220 

SPT 

230 

SPT 

260 

SPT 

250 

SPT 

260 

SPT 

270 

SPT 

280 

SPT 

290 

SPT 

300 

SPT 

310 

SPT 

320 

SPT 

330 

SP*’ 

340 

SPT 

350 

SPT 

360' 

SPT 

370 

SPT 

380 

SPT 

390 

SPT 

400 

SPT 

410 

SPT 

420^ 

SPT 

430 

SPT 

44 0 

SPT 

450 

SPT 

460 

SPT 

470 

SPT 

480 

SPT 

490 

SPT 

500 

Z‘ 

J'LO 

SPT 

520 

SPT 

530 

5 * 

■VO 


DIMENSION xnOOOO) 


c 

C-T^E FOLLOWING DIMENSIONS MUST 86 EQUAL TO OR GREATER THAN NBAN0S*2^1 

t'~^' ■ ' 

DIMENSION XBARI 100) ,STO< LOO) «6( 100) ,0(100) ,T1LQ0) 


C 

C-THE FOLLOidNG DIMENSION MUST BE 6QUAL„_T0_0R GREATER THAN _ 

C- M*M, WHE^E M=2=»N8ANDS>1 

C _ ■ _ 

OThENSIQN RX(IOOOO) 

_c 

C-THE FOLLOWING DIMENSION MUST BE EQUAL TO OR GREATER THAN (M^l)*M/2 
C 


DIMENSION R(5050) 

COMMON KREADR,KPRI NT _ 

DATA AHNAME/ *AMPL* t* ITUO' f *E ES • , • I MAT* , • E. 

I* •• */ 

DATA 


!• 


TNAME/»T-VA» ,«LUE*» ,• 
•/ 


DATA PHNAME/* PHAS* ,• € ES • , • T I MA • , • T E. 
l»*» •/ 

DATA BANAME/ »FREO* .'UENC* , *V ,-r 
DATA XNAME/* PERI* ,• OD *,* *,• 

KRE"A'DR = 5 

KPRINT=6 


( . t 


•t* 

t 


t 


•/ 

*/ 


READ IN CONTROL DIGITS AS FOLLOWS AND PRINT THEM OUT. 

JOOMAN-DOMAIN OF ANALYSIS. 0 OR BLANK* PER lOD, 1=FRE0UENCV. 

JBAN 0S-H3W VECTOR OF FREQUENCY OR PERI 00 BANOS WI LL BE OBTAINED. 

o' OR eLANK=GENER ATE ACCOROING' TO L'lMITS SUBSEQUENTLY READ’ IN. 
l*REAO IN A VtCTflR OF BANOS TO BE USED FROM CARDS. 

NBANOS-HOW MANY FRFUUENCY OR PERIOD BANOS WILL BE USED. 

NOBS-THE number of OBSERVATIONS IN THE DATA VECTOR 
JANALY-HQW THE DATA SHALL RE ANALYZED. 

0 OR BLANK=DO UNIVARIATE SPECTRUM, SEARCH FOR BEST K-VARIATE 

MODEL 

l=SIMPLY DO UNIVARIATE SPECTRUM 

2=00 K-VARIATE SPECTRUM ACCORDING TO MODEL PARAMETER READ IN 
JDATA-TELLS SUBROUTINE OFPVAR WHERE TO GET DATA. 

0 OR BLANK=R£AO FROM CAROS ACCORDING TO FORMAT <10X*7F10.0) 
1=SENFRATE DATA FROM MODELING PROGRAM CALLED BY DATA, THIS 


MODELING PROGRAM WILL ASK FOR MORE CAROS. 


ill 


jCRiTR-seiecTS criterion for optimi zation of k-vari&te spectrum^ 

0 OR BL A NK = MUL T IP L E™ R- S QUA RED 

l*T-VALu£S 

JPRINT-PRINTOUT CONTROL DIGIT ' 

0 = D0NT PRINT ANY RESULTS _ 

IxpRINT ONLY FINAL RESULTS 

2= PR I NT FINAL AND INTFRMEOIATE RE SULTS 

WRrNT” aT fOR 2 PLUS’ MArpICES, CORRELAT I ONS * ETC FOR EACH 

JPLQT-PLOT CONTROL DIGIT. ~ VALUES AS ABOVE EXCEPT FDR 3. 


^PT 

$PT 

SPT 

SPT 

SPT 

SPT 

SPT 


_ READ ( KREAOR . I ) 
i ’jPRlNTt JPLOT 
I FORMAT I16IS) 


STEPSPT 
SPT 
SPT 

JOOMANt J BANOS •NSANOSf NOBS tJ ANAL YtJDAT At JCRtTR* SPT 

SPT 

. iPT, 

SPT 
SPT 
SPT 

SPT 

IMPROVE- SPT 


WRITE IKPRlNTtl03) 

103 FORMAT (M *♦*♦♦* ♦ PROGRAM PARAMETERS 

WRITE ( KPRINTtl ) JOOMANt JBANDStNBANDStNOBStJANALYt JDATAt JCRlTRt 
1 JPRINTtJPLOT 


KTOP *MAX NUMBER OF TIMES THJ PROGR A M H I LL BE ALLOWED TO TRY 
MENTS ON The WHD^^ MULTIVARIATE SPECTRUM. KQUNT IS THE NUMBER OF 


C 

^-SUCH TRIES. 
C 

KTOP=^ 

kbUNT=0 


SPT 

_SPT 

SPT 

SPT 

SPT 

SPT 


SET SfGLEV* THE T-VALUE FOR THE CRITERION OF ACCEPTANCE OF 
THEN DECIDE HOW TO GET VECTOR OF BANDS. ; 


A PEAK. 


SIGLEV=2.0 
IF (JBANDS) 


100. 100, ilO 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


'READ IN LIMITS OF BANDS 
•OTHERWISE GO TD 1 10 AND 
BANDLQ-LOWeST FREQ. DR 
BANDEL-OELTAIFREU. I OR 
‘*BANOS-VECTOR OF BANDS. 


VECTOR IF JBANOS WAS 0 - GENERATE BANOS. 
READ IN BANOS FROM CARDS. 

LONGEST PERIOD, AS APPROPRIATE TO DOMAIN. 
OELTAIPERIOOI, AS APPROPRIATE TO DOMAIN. 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


100 

2 


120 


110 


READ (KREA0R,2I B ANDLO, 8 ANOEL 
FORMAT (8FID.0) 

BANDS! 1 I =BAMOLO 
00 120 LBANO= 2.N8ANDS 
BANOS! L band ) = BANOS( L8AN0-1 I 
GO TO 130 _ 

R E AD U REAOR , 2) I BANDS! L BAND) - 


BANOEL 


'LBANO'! .NSANOS) 


SPT 

SPT 

SPT 

„ SPT 

SPT 

SPT 

SPT 
SPT 

CONVERT BANOS TO PERIOD DOMAIN IF ORIGINATED IN FREQ DOMAIN. SPT 

•REVERSE the ORDER OF BANDS SO AS TO BE IN ASCENDING WAVELENGTH ORDER. SPT 

■ALSO GENERATE BAN, THE FREQ DOMAIN VECTOR OF BANOS. SPT 

SPT 

SPT 
SPT 
SPT 
SPT 


I 30 
I AO 


( J D n M A N V 150,150,140 


141 


IF 

00 141 IBAND= I , NPANDS 
BAN! LB AND! = 3 ANDSI LtlANOl 
BANDS! L band ) = I . 0/ BANDS ! L BAND I 


CALL REVI BANOS, NSANOS) 


C ____ 

X'-GENERATe' S ! N/COS" PREDICTOR MATRIX AS PART OF DATA MATRIX, X. THEN 


SPT 

SPT 

SPT 


550 
’560' 
570 
580 
590 
600 
_610 
6 ' 20 ’ 
630 
640 
650 
66 0 
_6J0 
630 
690 
700 
710 
720 
_730_ 
740 
_ 750 
76 0 
770 
780 
_T.RO- 
aoo 
810 
820 
830 
840 
_850_ 
860 
870 
880 
890 
900 
„910_ 
92 0 
930 
94 0 
950 
960 

930 
990 
1000 
1 01 0 
1020 
1030 
104 0 
1050 
I 06 0 
1070 
1030 
1090 
'l 10 0 


“Iv 

OJlKcfNAL' PAGE IS 
OP POOR QUAUTSJ 



C-ft eAO IN oe PENDENT VARIABLE PART O F X. THESE D A TA W!LL TH EN BE ENT- S PT 

C-6RED'TNTlJ"'SU8ROUfr^^^ CORRE TQ COMPUTE THE CORRELATION MATRIX OP ALL SPT 

C-VARIABLES. SPT 

c‘ SPT 

150 MsNBANOS^Z I . SPT 

CALL PREGENIX.flAN0S*N08S,N8ANDS ) SPT 

CALL OEPVAR t X^NOBS, JDATA, JDOMAN»JPRINT| _ SPT 



C-GENERATE VECTOR OF VARIABLE NUMBERS FOR LATER USE* SPT 

t~ ' ‘ ' SPT 

00 155 LV=l.M :_SRT 

i55LVAR(L V)=LV ' SPT 

CALL CORRE(NOftS«M« I, X ♦ XB A R > STP , RX,R < 6 tO * T I SPT 

"C SPT 

C-SET TOPP, THE MAXIMUM ALLOWABLE TOTAL AMPLITUDE IN THE FINAL SPECTRUM SPT 
C'-TOPP IS THE total AMPLITUDE OF THE FlTI^aS PERCENT* SPT 

C SPT 

T6»*ST0( l)V. 25*STD( 1 1 SPT 

C _ - SPT 

C^PRINT OUT MATRICES ETC.* IF DESIRED SPT 

C SPT_ 

rF(JPRINT-3) 151*152^152 SPT 

152 WRITE IKPRINT,5) _ ^ ; SPT 

~ “5 FORMAT (• IM A T R IX OF C 0 R R E L A t I 0 N S • . / » * VARIA8SPT 
ILE I IS The dependent variable* ALL OTHERS ARE COSINE AND SINE WAVSPT 


2ES OF VARIOUS PERIODS •,/,T88** 
CALL MATPRTIR.M.LVAR) 




THEN FIND 


C-OeCIDE ON PROPER DATA ANALYSIS PATH. 

~C~ 

151 JANALY = JANALY*1 

GO to I 200, VO'0,6bd)* JANALY 

c _ _ 

C 

_c 

c 

C-STATEHENTS 200-^00 INVOLVE COMPUTATION OF OPTIMIZED K-BANOS SPECTRA 
C 

C-FIRST COMPUTE UNIVARIATE SPECTRUM USING SUBROUTINE USPECT. THEN FIND 
C-ALL SIGNIFICANT SPECTRUM PEAKS IN UNIVARIATE SPECTRUM. 

C _ _ __ 

200 CAL L USPEC T ( R',T i>40BS. M*6 * NBA NDS , S TO, J PR I NT* JPLOT * PHAS E* 8 AN* BA NOS 

1 *AMNAMEt TNAME • PHNAME , XN AM E* JDOMAN * B ANAME ) 

CALL PPIKR (B,T, ISAVE,NPEAK*NBAN0S*SIGLEV » 

IF (NPEAK) 201,201*202 

201 write (KPRINT,66I 

66 FORMAT CO** *_* ♦ NO SIG NIFICANT PEAKS FOUND ♦ ♦♦♦»«) 

- GO ”T3' 2000 
C 

C-NOW LOOP THRU AN OPTIMIZING PROCESS TRYING TO FIND LARGEST VALUE OF 
C-EITHER T or R-SUUARE. THE BEGINNING OF THIS LOOP IS STATEMENT 210 
C 

C-SET A VALJE OF CRIT = 0. THIS WILL LATER BE USED TO STORE_ THE JUST 

■^PRECeo'lNG'VALUE OF THE' DPTImTZATION CRITERION. ALSO’ INITIALIZE THE 


V 


C-variaBLE =-1. mm is the DIR ECnON OF PEAK SLIDING AS WELL AS THE 
C“C0uFtER“'F3R THE NUMBER 6”F STE>S. 

C _ ■ - 

202 LPEAKsl ■ 

CRlTaO.O 

MM=0 

210 CALL KSP6CT (R ,T , NOfl$ »M^NPEAK« !SAVE< RK,Bt JPR INT# XBARtSTOtPtRSO_> _ 
i^hTseV BAV0S”tBA 

C-TEST FOR PRINTOUT ~~~ 

C : 

iF(jPR I NT-3 >208,209, 209 

209 WRITE <KPRlNT,fe) RSQ ' ~ 

(TTQ'RM'Al^i ■OMUL flP^^ SQUARED CORRELATION** ,F8.4) 

WRITE (KPRINT.7) 

FORMAT! • PER lOO/FREO* • ,T20, • AMPLITUDE • ,t40,*T~VALUE* I 

00 214 LO=l,NPEAK ^ 

“ kK»ISAVE(LO> 

IF (JOOMANJ 212, 212,213 

212 BB-BANOSCkk) 

GO TO 214 ... 

213 BB=BAN(KK> 

214 WRITE (KPRINT,8> RB, B (KK » , T I KK I , 

8 FORMAT I 2F12.4,F18.4> 


SPT 1 670 
SPT 1680 
SPT 1690 
SPT 1700 
SPT 1710 
SPT 1720 
SPT 1730 
SPT 1740 
SPT 1750 
SPT 1760 
SPT 1770 
SPT 1780 
SPT 1790 
SPT 1 80 o' 

SPT laio 

SPT 1020 
SPT 1830 
SPT 1840 
SPT 1 850 
SPT 1860 
SPT 1870 
SPT 1880 
SPT 1890 
SPT 1900 


C SPT 1910 

C-SELECTCRTfERION F^R GOODNESS OF FIT - CHECK FOR rMPROVEME!^ ~SPT 1920 

C . ^ SPT 1930 

208 IF (CRlTR(RSQ,T,JCRITR,LPEAK»-CRltT 226,220,230 SPT 1940 

C _ SPT I960 

C-A BRANCH TO 220 IMPLIES NON-IMPROVEMENT* EITHER RESORE OLD VALUE AND SPT I960 

C-GO TO NEXT PEAK 1221 OR 223), OR TRY MOVING UP C222U .SPJ_J9T0 

^ 

220 IF < MM 4^1) 221*222,223 SPT 1990 

SPT 2000 

C-A BRANCH TO 221 IMPLIES THAT WE HAVE BEEN MOVING DOWN FOR 2 OR MORE SPT 2010 

C-STEPS. THIS IS THE FIRST NON-IMPROVEMENT* RESTORE JUST FORMER VALUE SPT 2020 

C-AND GO TO NEXT PEAK. SPT_2030 

— SPT 2040 

221 ISAVEI LPEAK)»ISAVEaPEAKI^l SPT 2050 

MM=-1 - - . - - ^ 2060 

LPEAK=LPEAK«-1 ... ; -SPT 2070 

ISAVEI LPEAK)=1SAVE(LPEAK»-1 SPT 2080 

GO TO 300 ' - S PT 2090 

— 2100 

C-A BRANCH TO 222 IMPLiES FIRST MOVE DOWN AND NON-IMPROVEMENT. TRY SPT 2U0 

C-MOVING SAME PEAK UP SPT 2120 

C „ SPT 2130 

^ 222 MH*l " " ■ "" SPT 2140 

ISAVEI LPEAK 1= ISAVEILPEAK ) 4-2 ___SPT 2150 

' i'E"' II SAVEILPL AK )-NBANOS ) 300*225*225 SPT 2160 

225 WRITE (KPRINT,)! SPT 2170 

3 FORMAT (<0** * ♦WARNING* - JUST TRIED TO MOVE A PEAK I NSPT 2180 

ITO THE LONGEST PERIOD (LOWEST F REQU ENCY ) BAND* TRY WIDER SPECTRUMSPT 2190 

2 LIMITS. * ) SPT 2200 

GO TO 2000 SPT 2210 

- - SPY 2220 


OP page IS 

OP POOR QD Ar.rp y 


Vi 



on no r> on o!n 


C-A BRANCH TO 223 IMPUES THAT_W|^ 
"C^ReSTO'RiE JUST PRIDR^VALUE AND GO 

C 

223 [SAVEC LPEAK ) = ISAVE(LPEAK>-1 
MM=:-l _ 

lpeak*lpeak+i 

ISAVetLPGAK » = ISAVeaPEAKI-l 


HA VE BEE N MOVING UP AND IT DIDNT 
TO NEXT PEAK. 


GO TO 300 


-A BRANCH TO 230 IMPLIES IMPROVEMENT. CONTINUE IN SAME DIRECTION. 


230 IE (MM) 231. 231. 232 


HELP. SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

S^T_ 

SPT 

SPT 

SPT 
__SPT 
SPT 
SPT 


-A BRANCH TO 231 IMPLIES WE WERE MOVING DOWN. “ CQNT-I NUE_SAM„^^ 


^31 MM»MH-l 

CRiT*CRITR(RSO,t,JCR ITR.LPEA kT 
ISAVE(LPEAK) =ISAVEILPEAK)-1 


SPT 

OIAECTIONSPT 

SPT 

SPT 

SPT 

SPT 


IF (ISAVEILPEAK )-l ) 23^,234.300 

tKPRlNT.4) 

4 F0RMAT(*0*« * ♦ ♦ W A 

10 THE SHORTEST PERIOD 
2M LI MI TS* ) 

GO TO 2000 


R N I N 
(HIGHEST 


SPT 

_ _ . SPT 

TRIED to MOVE A PEAK INTSPT 
bano_*„ try wider SPECTRUSPT 

SPT 

SPT 


G ♦ - JUST 
FREQUENCY! 


-A BRANCH TO 232 IMPLIES WE WERE MOVING UP. - CONTI NUE SAM E DIRECTION. 


232 MMsMM^-l 

~ CRlT=CRITR(RSQ,T,JCRITR,LPEAk> 
I SAVE I L peak )=I save (L PEAK) 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


IF 1 1 SAV6(LP6 AK)-NBANDS ) 300,225.225 


STATEMENT 300 CHECKS TO SEE THAT WE ARE NOT OFF THE END OF THE ISAVE 
VECTOR. IF NOT. DO ANOTHER K-BANOS SPECTRUM WITH NEW ISAVE VECTOR 
■VALUES. IF WE ARE OFF THE END, COMPUTE A NEW FULL SP EC TRUM USING A_ 
A-VAR i ATE ' MODEL. ALSO CHECK FOR" VALUES TO EACH OTHER. 

■IF WE HAVE JUST MOVED INTO A ANOTHER PEAK, REDUCE NPEAK BY ONE AND 
PACK THE VECTOR DOWN AND START PROCESS OVER AT LP£AKs l. 

300 IF (LPE AK-NPEAK ) 320,320,310 " 

320 IF (LPEAK-l) 322,322.321 

32r"lF nSAVEILPEAK )^1SAV6(LPEAK«U > 

322 IF (LPEAK-NPE AK ) 324,210.324 

324 IF ( I SAVEILPEAKI-ISAVEILPEAK-i-in 

323 00 330 LL = LPEAK .NPEAK 
330 ISAVE(LL-l)^lSAVE(LL) 

NPFA<*NPEAK-l 

■ L P E AK » I 
MM = 0 

CRIT=0.0 
GO TO 210 

325 NPEAK=NPEAK- I 

00 331 LL = LPf.AK, NPEAK 

“331“ I S 4ve < L L ) = 1 S AV E ( LL *1) 


322,323,322 

210,325,210 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPJ 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

spt 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


2230 
2240 
2250 
2260 
2270 
2280 
^2*5 0 
2300 
2310 
2320 
2330 
2340 
2.350 
2360 
2370 
2380 
2390 
2400 

2420 

2430 

244 0 

245 0 
2460 
247^ 
24B0 
2490 
2500 
2510 
2520 
_25JL0 
2540 
2550 
2560 
2570 
2580 
2590 
2600’ 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2 720' 
2730 
2740 
2750 
2 760 
2770 

“2 78 0 



non • n o o 


C 

C- 

C- 

"C^ 

C- 

”C 


LPEAK=1 . — 

Tim-0 

CRITsO.O 

GO TO 210 

A BRANCH TO 310 OCCURS WHEN OPTIMIZATION OF THE K-BANOS . 

COMPLETE. NOW COMPUTE A FULL SPECTRUM USING_jriHEj3PT^lZ^0 — 

NEW i;-r sr6CTRUM->orPEAKr^-r^^^ IF ANY NEW ONES lUR 
NED UP- IF SO. GO THRU WHOLE OPTIMIZATION ROUTINE AGAIN. 

310 IF ( JPRINT-3I360.36l.36l ~ 

361 WRITE (KPRINT.qi 

9 FORMAT! * ICQMPUTING SPECTRUM 


SPT 2790 
SPT 2800 
SPT 2810 
SPT 2820 
SPT 2830 
SPT 28^0 
_SPT_ 2850 
-SPT '286 0 


USING THE FOLLOWING BA NDS !N MODEL,* I 


T F ( JOOMAN ) 362. 362.363 

362 00 366 LO=l*NPEAK 

KK^ISAVEILOI 

366 WRITE (KPR1NT.8) BANDS! KK I 
GO TO 360 

363 00 366 LO=l.NPEAK 


KKMSAVEILOI 

366 WRITE (KPRINT.8I BAN(KK) . 

"360 CALL MULTI! I S AVE.NPEAK.N8 ANDS.R.T.NOBSt M» RX# B. XBAR. STD. OtRSQ# 

_ I PHASE f SANDS. BAN) — — 

kOUNTsKOUNT+l 

^ *L I 0 U N T- K T OP) 38 4 . 384.385 — 

386 FORMAT!* ITHE PROGRAM IS PROBABLY CAUGHT 
IINPUT PARAMETERS AND SUBMIT AGAIN*) 

GO TO 2000 . . ^ — 

CALL PPIKR!8*T»ISAVE. NPEEK. N BANDS .S IGLHV ) 

IF (NPEAK-NPEEK ) 311 1 350 . 3 il 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


384 


311 


NPE“AK = NPEEK 
LPEAK-l 
MM = 0 

CR IT=0.0 
GO TO 210 


rS^tTBrr'STA'TE IN the" OPTIMIZING ROUTINE RESUTS IN A BRANCH TO 350. 
‘OUTPUT RESULTS.- 

350 WRITE (KPRINT.671 > c r t d h Mi L 

67 FORMAT! *0 MULTIPLE BAND SPECTRUM) 

CALL OUTPUTCPANnS.&AN.B.T . P HAS E .NBANDS . AMNAME . TNAM E .^HNjAMEj, 


IS A NAME . XNAME . J PH I NT. J PL 6 1 * JDOMAN) 
WRITE IKPRINT.6 ) RSO 


-CHECK FINAL SPECTRUM FOR EXCESSIVE TOTAL POWER IN SIGNIFICANT BANOS. 


CALL TOTAM J X . NOBS . NPE AK . I SAVE . M. PHASE . B » SUMS ) 

IF (^UMB-TrjPP ) 3 81.382. 382 

302 WRITE (KPRINr.383) 

383 FORMAT! • IPRJGRAM FAILED BY FINDING TOO MUCH POWER 
SUGGEST ALTERING PARAMETERS OF THE 
2* USE WIDER SPECTRUM LIMITS. FILTER 
381 GO TO 2000 


IN THE SPECTRUM 
ANALYSIS. SUCH AS* 

DATA, ALTER OELTAP, ETC.*) 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

JPX 

SPT 

SPT 

SPT 

SPT 

SPT 

_SP_T 

Spt' 

SPT 

SPT 

SPT 

SPT 

_S_PT 

SPT 

SPT 

.SPT 

SPT 

SPT 

SPT 


2970 
2080 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
_2970 
298 0 
2990 
3000 
3010 
3020 
_3030 
3040 
3050 
3060 
3070 
3080 
_309 0 
3100 
3110 
3120 
3130 
3140 
_315 0 
3160 
3170 
318 0 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
'323 0 
320C 
3300 
3310 
3320 
3330 
3360 


vlii 


page jb 

OF POOR QGALIT 3 J 
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-STATEMENTS AOO-600 INVOLVE COMPUTATION OF THE SIMPLE* SINGLE BAND 
-SPECTRUM O NLY* 

400 CALL USPECT I R ,T ♦ NOBS , M, B • NBANDS *ST0* JPR I NT , JPLOT* PHASE# BAN »BANDS* 

I AMNAME*TNAME.PHNAM6*XNAHE* JOOMAN,8ANAMEI 
GO TO 2000 

**«*«*♦ *V*V*'**<t ♦♦*♦♦♦★*♦*♦♦♦♦♦♦♦♦♦*♦♦♦*♦♦♦****♦♦♦♦*♦♦***** *****^ 


-STATEMENTS 600-800 INVOLVE COMPUTATION OF A SPECIFIC K-BANO MODEL 
-WITHOUT OPT IMIZAT ION, 

-FIRST READ IN MODEL PARAMETE RS . BUI LO yECTOR_ OF BANDS I N MODEL * THEN 
-COMPUTE SPECTRUM. 

600 READ (KREAORfU NPEAK 

READ (KREADR.2) ( PEA KS 1 1 PE AK I »LPEAK«l*NPEAK I 

DO 610 L>EAK=i* NPEAK 

00 620 L6AND= UN8AN0 _ _ 

IF ' (PE AKSILPFAKI-pANOSILBANOn 620*630,620 

630 1 SAVE! L PEAK > =L BAND . 

GO TO 610 

62^_C0NTINUE 

6 id continue" 

CALL multi ( IS AVE, NPEAK, NBANDS, R*T »NOBS * M* RX* B* XBAR* STP, 0,RSQ, 

I PHASE .BANOS. BANI ■ 

WRITE (KPRINT.5) _ 

CALL OUTPUT ( BANDS .BAN, B.T* PHASE , NBANDS* ANAME, TNAME,PHNAME# 

I B ANAME. XNAME , J PR I NT , J PLOT, JOOMAN I 

W”RltE VkpRINT.6) R'SQ 

2000 STOP „ _ _ __ — 

END ~ _ 



r»:r> 


SUBROUTrNe RHV(X,NI 
OlMENSinM X(50) 


REV 

^ " ■"■rev 

REV 

'■■■ REV. 

C „ _ _ ^ REV 

"C REV 

c * _ _______ * 

C * SUBROUTINE REV ’ ' ' ♦ REV 

C ♦ ^ * REV 

C * REV TAKES A VECTOR AND REVERSES ITS OROER. '• REV 

C • ♦ R EV 

~C *1>** #♦**♦♦*♦★'•»*♦♦**•♦♦♦*♦♦♦*♦♦♦♦**♦♦**♦♦**♦♦***♦*•♦****♦*♦♦♦♦*♦♦♦♦♦♦ REV 


J-N»l 

MIO*N /2 

bo to L=l.MID 

J»J-l 

SA'VE=»X(L) 

X(tl-XU) 

lO X(JI=SAVE 

RETURN 

END 



REV 

REV 

_REV 

REV 

REV 

REV 

REV 

REV 

REV 

REV 


10 
20 
30 
40 
_5 0 
60 
70 
flO 
90 
100 
HO 

lid 

130 

140 

150 

160 

_l_tP 

IBO 

190 

200 

210 

220 



SUBROUT INE PREGEN i X* B AMDS * NOBS* NB AMDS J 
DIMENSION XU )» BANDS m 




SUBROUTINE P R E G E N 

PREGEN GENERATES THE SINE AND COSINE PREDICTOR WAVES AND STORES 
THEM IN MATRIXJ^. MATRI X X IS_THE J)_^A HATRt X HAV ING NOB S ROWS 
AND H COCUMNSr HERE NOBS* THE NUMBER OF O0SERV AT IONS AND 
BANDS IS A VECTOR OF PERIOD VALUES IN THE SPECTRUM, _ 


♦ 

■ * 

— ■— ■■ ■ : 

MATRIX X IS STORED IN VECTOR MODE. SEE IBM SYSTEM/360 SUB- ♦ 

routine package I 360A-CM-03X) VERSION lilt PROGRAMMER'S MANUAL* ♦ 
PU 8LI T IO N NUMB E R 2 0 5-3* PAG E S 3-6 , ESSENT IA LLY EACH CQt^ ♦ 

_ ... ENO-TO-ENO into one long VECTOR* ♦ 


UHN OF MATRIX X 


* 

♦ ... _ ._. - 

♦THE FIRST NOBS VALUES OF VECTOR (MATRIXI X WILL BE THE DEPENDENT ♦ 

♦ VARIABLE (SEE SUBROUTINE OEPVARI, THE NEXT NOBS POINTS (COLUMN ♦ 

♦ YWO) WILL CONTAIN COS WAVE* BAND l', THE NEXT NOBS 

♦ THREE CONTAINS SIN WAVE, BAND I* ETC* 


POINTS (COLUMN 


♦♦♦♦♦♦ *♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦**♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦ ♦♦♦♦♦♦♦♦♦♦♦*♦♦♦•♦ 


C 
c 
c 

C-CONVERT TO RADIANS/OBSERVATION AND LOOP THRU NBANpS TIMES* 

c ' " ■ 


L»NO0S 

"■^0^ 1 0 “ L B A NO = i' * N 8 A NO~S 
AFREO=6.2a3l8/0ANOS(LBANOI 
EF*-AFREO 
00 20 LOBS=1.N08S 
L*L^l 

EF = EF«-AFREQ 


20 X(L)=COS(EFr 
EF--AFREO 
00 10 LOBS=l.NOBS 

L = L4^l 

EF=EF*AFREO 
10 X(LI=SIN(EF) 


X CONTAINS THE VECTORS OF SIN/COS WAVES STRUNG ENO-TO-EnO* THIS IS 
EQUIVALENT TO A MATRIX WHERE EACH SIN OR COS WAVE FORMS A COLUMN OF 
NOBS LENGTH. 


RETURN 

■END 


PREO 
PREO 
PREO 
PREO 
_PRED 
PRED 
PREO 
PREO 
PRED 
PREO 
^REp_ 
PREO 
PREO 
PREO 
PREO 
PREO 
PRED 
PREO 
PREO 
PREO 
PREO 
PREO 
PREO 
PRED 
PREO 
PREO 
PREO 
PREO 
PRED 
PRED 
PREO 
PREO 
PREO 
PREO 
J^RE_D_ 
PREO 
PREO 
PREO 
PRED 
PREO 
_PREO 
p R e b 
PREO 
PREO 
PREO 
PRED 
_PREO 
pR^b' 


10 

20 

30 

40 
50 
60 
70 

ao 

90 

100 

uo 

120 

130 

140 

150 

160 

ITP_ 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

2B0 

290^ 

300 

310 

320 

330 

340 

350 

360 

370 

300 

390 

400 

41 0 
420 
430 
440 
450 
460 
470 
480 





xi 



o o 


1 


SUBROUTINE OEPVARI Xf NOBS* JOATAtJOOMAN* JPRINTJ 

OIMENS ION X ( n 

COMMnN KRFAOR 

'♦ 

* SUBROUTINE D E P VAR 

* D6PVAR EITHER RE&OS IN 0R_GEN ERATES IVIA S I MUL A I «_ A_DE PEj£p6J 

“i ^‘RTABCE vector nhF'~LE~NCTH DEPENDENT VARIABLE IS 

* INTO THE FIRST NOBS CELLS OF X» A MATRIX OF VARIABLES STOREI 

* “Vector of columns, see subroutine pregen for rest of x. 

-IF J0ATA'“TS 0 - READ FROM^C/TROS* (10X,7F10.0» 

-IF JDATA IS I - CALL SIMULA WHICH GENERATES DATA USING A MONTE 
-SYSTEM ALONG WITH DETERMINISTIC DATA. SIMULA READS CAROS. 

— IP (jqATA) 10,10,50 

10 READ IKREADR,!) C X( L ) >L~l , NOBS ) 

I format { I OX , 7F I 0 . b I 

GO TO 20 

50 CALL S IMULAI X.NOBS* jDATA,JDOHAN#JPRlNT) 

20 RETURN - 

END 




NT_ ^ 

stored* 

0 AS A* 


*«#«♦** 


CARLO 






o o 


SUBROUT INE S I MUL A ( X, N , J , JODM^N* JPRt NT » 
OIME^4SION X( n ♦ AU30 > ,BAND( 1 I 
CoHMC'ii 


NrvC«L>iV«r\>'i\iln| 


SUBROUTINE S I HU L A 

SIMULA IS A PRIHATIVE MONTf CAR LO ANO SIGNA L GENER ATOR FO R CREAJ:^ 



SIML 

10 


SIML 

20 


SIML 

30 


SIML 

40 


SIML 

50 


SIML 

"60 

** 

SIML 

70 

♦ 

SIML 

80 

♦ 

SIML 

90 


SIML 

100 


ARtTFiCIAL bATA- GAUSSIAN DATA ARE READ IN FROM CARDS AND 
OF SINE WAVES OF VARIABLE FREQUENCY AND AMPLITUDE ARE 


« ING 

*_ A number -- - - 

* ADO ED. * 

* . ... *- 


T 

C 
C 
C 
c 

c 

d-ZERO VECTOR X AND PRINT HEADING 



bo 

100 

X( 


IF 

120 

WR 


100 L=l ,N 
)= 0.0 

(JPRINTJ Il0,ll0tl20 

TE (KPRINT.12M 

- 


121 FORMAT ( * IM 0 N T e 


READ IN MODEL SPECIFICATIONS - GENERATE 
NSIN-NUM8ER OF SINE WAVES IN DATA 
‘ BAM0-FR60. OR PERIOD OF THE LSINTH SINE WAVE 
SO-ST ANOARO OEFI^TIONOF^T HE_ J.S I NT H SINE WAVE. 
SbN-S T and ARD OE V . OF GAUSS I AN NO I S E 


CARLO SIMULATED OATA*,/.*0») 
SINES ANOHivob ‘ NO I SE 


(AMPLITUDE) 


ilO READ (KREADR.U NSIN 

1 FORMAT (151 

~ 00 10 LSIN=1,NSIN 

READ (KREA0R,2) BA NDU) 

2 FORMAf (2FlO.br 

IF UDOMAN 1 20* 20, 30 
30 BAN=BAND(l) 

BANOI l)=l,0/BAND( 1) 

' GO TO 2 1 

20 BAN=l . C/BANDf I ) 

^ , ^ ^ PRINT) 2 5 1 2 5 * 2 2 

22 WRITE (KPRINT,3) LSIN, 

3 FORMAT ( • BAND* » 1 3 
ITUDE= ' ,F 10.4) 

25 AFREO=6.2y3l 8/BAN0( I ) 

EF=-AFRtO 

L ^ 1 


SO 


BAND! Ut BANf SO 
- PERIOO=* tFlO.4* *t 


FREQUENCY^* f Fi0.4f •# AMPL 


DO 26 L^l.N 
£F=EF^AFRF0 

26 AIL ) = SIN( CF ) 1.41421 

00 10 L=l.N 
10 X( L) =X( U »A{ L MSO 

C 

n:-R'EAO“^N"Naise~sb and noise Caros 


ORIGINAL PAOB Et 
OF POOR QUALnU 





xiii 


READ (KREADR,2» SDN 
IF<SDN> 1000*1000*40 

40 IF IJf'RINT} 50,50,41 ' 

41 WRITE { KPR |NT,5» SDN 

5 FORMAT ( ’CRANDOM GAUSSIAN NOISE AHPLITU0E»*,F10.4I 
50 READ (KREAORffel (A(L»,L=l*N» 

6 FORMAT! 4X* 10F6*4 I 
00 60 L=1,N 

60 X(LI=A( L»*SDN^X(Ll 

100 RETURN 

END 


SUBROUTINE MATPRT«R,M,LVARI 
DIMENSION Ril)*X(l5) .LVARI l» 
COMMON KREAOR, KPRI NT 


- - ■ — ♦ 

♦ * 

• SUBROUTINE M A T P R JT ^ 

MATPRT PRINTS OUT TH E LOWER TRIANGULAR MATRIX X WHE gg X IS STOR-._^ 


HATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 


■* to IN MODE 1. SEE REFERENCE L ISTED IN SUBROUTINE PREGEN* ^ 

"tv** 

’LOOP THRU AS MANY TIME S AVNEEDEO^^T^^ "WHOLE MATRIX, EACH TIME 

-PRINTING 15 COLUMNS BY NRO ROWS* 


MATP 

MATP 

MATP 

MATP 

MATP 


N2*15 

110 IF IN2-M) 

16 N2*M 

100 WRITE ( KPRINT, U 


100 , 100,10 

ILVAR(LV),LV*N1,N2I 


MATP 

MATP 

MATP 

MATP 

MATP 

MATP 


1 FORMAT (• 0* , 2X, 151 7! 

WRITE I KPRINT, 2 I 

2 FORMAT! »0» ) 

PRINT ONE ROW AT A TIME.” MATR^ MODE I, UPPER TRIANGULAR 

THEREFORE REVERSE SUBSCRIPTS TO MAKE LOWER TRIANGULAR^ 


MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 


00 20 LRQW=Nl,M 
" IF (LR0W-N2I 21,21,22 

21 limit=lrow _ 

"GO TO 2 3 

22 LlMir=N2 


23 


30 

20 


L = 0 

DO 30 LCQL=N1, LIMIT 

L=L + l 

CALL LOC! LC0L,LR0W,K,M,M,1 1 

X(L)=R(K I 

WRITE IKPRINT.3I (LRQWt! X(K),K^ 1>L I I 


T0R"m A T I U , i 3 , 2 X , 1 5F 7 . 3 ) 
Nl=Nl*15 
N2*Nl^-l^ 

IF (M-NIJ IQOC, 110,110 
1000 RETURN 

END 


MATP 

MATP 

MATP 

MATP 

JIATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 

MATP 


to 

20 

30 

AO 

_50__ 

60 

70 

80 

90 

100 

U0_ 

120 

130 

lAO 

150 

160 

JJQ_ 

idb 

190 

200 

210 

220 

2^0_ 

2^0 

250 

260 

2T0 

230 

300 
310 
320 
330 
340 
350 
360 
370 
38 0 
390 
400 
410 
420 
430 
440 
450 
460 
470 


XV 
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USPT 

SUBROUTINE USPECT(R.T,N , M, B t NBANOS f S TDt J PR INT t JPtOTf PHASE# BAN» USPT 

I BANOS, AMNAMt.TNAME.PHNAME* XNAMEfJOQMANtBANAME) USPT 

0 I ME NS I .ON R 1 1 J , T ; I I , ^ U t S TDI IKPHASE I U , BANU >» BANDS ( U # AHNAME < 

TTr NAME ( ITTP HNAM E r UV'X NAM'ETI ) , BA NA^ME TI ) 

COMMON KREAOR.KPR INT _ . - 


V ■■ " * 


SUBROUTINE U S P E C T 


USPECT COMPUTES THE UNIVARIATE SPECTRUM BY MULTIPLYING THE CORR^ 
ELATION OF THE TIME SERIES WITH A SIN OR COS PREDICTOR 8Y THE 
RATIO OF STANOARO DEVIATIONS OF THE TWO. THIS IS DONE FOR EACH 
PREDICTOR. LET R-THE CORRELATION OF F(TI WITH A SIN OR COS WAVE 
OF SOMF. ARBITRARY WAVE LENGT H, LET S >TH E STANO_ARO OEXIALIQ.^- 
lfF~tHE ‘T lME“TERlEr AN^^ .TOT^tHE STANDARD OEVIATTON OF THE PRED- 
ICTOR WAVE. THEN A(FI = R^(S/.707) .. .. 

WHERE A(FI=THE AMPLITUDE OF THE COMPLEX SPECTRUM AT F. THE SI N 
AND COS COMPONENTS ARE COMBINED AT EACH WAVE LENGTH AND PHASE 
ANGLE IS COMPUTED, T-VALUES ARE ALSO COMPUTED FOR EACH BAND. 

R»UPPER TR I ANGUL AR CORRJL A TI ON M ATR IX , IN STO R AGE MODE I, _S£E 

T^^pObL rCATION 'RE IN SUBROUTINE PREGEN. (INPUT) 

T^VECTOR OF T-VALUES. (OUTPUT) . 

NOBS=NUMBER OF OBSERVATIONS, UNPUT) 

M= TOTAL NUMBER OF DATA VECTORS. M=NBAN0S*2 1. ONE VECTOR 
THE DEPENDENT VARIABLE AND TWO FOR THE SIN/COS WAVES AT EACH 
LENGTH. (INPUT). B* V EC T OR OF S PE CT RU MJES T I MA TES. (OUTPUT) 
'NB'ANDS=LENGTH‘ OF vectors f ,8 AND >HA$E, (INPUT) 

STD-STANOARO DEVIATIONS OF ALL V AR T A9L ES . ( INPlTT ) 

JPRINT=PRINT CONTROL - SEE MAINLINE. (INPUT) 

JPLOT = PLOT CONTROL - SEE MAINLINE. (INPUT) ___ 

PHASE=VECTOR OF PHASE ANGLES. (OUTPUT) 


* 

.. 

* 

FOR * 
WAVE* 

♦ 


USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 

USPT 


■***V*1rJiri*<riV**T* *♦♦♦*♦* *♦♦♦♦♦♦*♦♦♦♦♦♦♦**♦*♦**♦***♦*♦***♦♦***♦*♦♦♦♦♦♦♦ 

L B A NO s O , . .......... — 

extract VECTOR OF CORRELATIONS BETWEEN DEPENDENT VARIABLE AND EACH 
•PREDICTOR VARIABLE. 


10 


DO 10 L=2,M 

CALL LOC( 1,L.K00L,M,M,1) 
8(U-R(K00L ) 


-LOOP THRU THE SPECTRUM COMPUTATION NBANO TIMES. 


FACT=SORT( FLOATlN-2) I 

DO 20 L=2.M,2 

K*L>1 

LB AND = LB AND U 


CO=COS 
•SI=S IN 


COMPONENT 

COMPONENT 


OF 

0 F“ 


SPECTRUM 

SPECTRUM 


ESTIMATE 

eStlMATE 


* USPT 
♦_^SPT 

♦ USPT 

♦ USPT 

* USPT 
» USPT 

USPT 
USPT 
USPT 
USPT 
USPT 
iJS^P^T 

us> f 
USPT 
USPT 
USPT 
USPT 
JJSPJ 
USPT 
USPT 
USPT 
USPT 
USPT 
JUSPT 
USPT 


10 

20 

30 

40 

_50_ 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
_4I0 
42 0 
430 
440 
450 
460 
_470 
480 
4.9 0 
500 
510 
520 
530 
540 


xvi 


* 


OMGiNAL PAGE IS 


b O ;u o o lo 


C-tC*T-VALU6 for cos componen t 

■C^r$=T- VALUE FOR S In' COMPONENT 

c 

CO»B<L ) ♦STO( I )/.707 

SI = B(K »^STO( 1 ) /.707 

TC=B(L )‘FACT/SQRT( l.0“BCLI**2J 

TS=B<K ) <^FJV C T/SO RT < l.Q-B( K )**2 > ■ 

COMPUTE SPECTRUM AMPL ITUOE* T-'VALUE AND PHASE ANGLE 

B(LBANO> = SQRT (CO**2 SI**2) 

T(L8AND)=S0RT (TC^*2 ♦ TS*>2I 
20 PH ASE C L SAND) =AT AN CCO/S! ) 


PRINT AND/OR PLOT IF DESIRED 

JPR=JPRlNT-l 

'JPL*JPL0T-1 

IF < JPR >_50t 50,^ — 

50 I F ( J P'L) ICO Of 1000,60 

60 WRI TE ! KPRI NT, 1 ) 

r FORMAT! *0 SINGLE BAND SPECTRU M*) 

CALL OUTPUT ( BANOS , BAN, B,T * PHASE ,NBANOS* AMNAME» TNAME *_PHNAM^_ 
I/BANAME ,XNAHE, JPRINT, JPLOT, JDOMANI 

1000 RETURN ^ 

END 


JJSPT,_550„ 
OSPT 560 
USPT 570 
USPT 580 
USPT 590 
USPT 600 
_USPJ1610_ 
USPT 620 
USPT 630 
USPT 640 
USPT 650 
USPT 660 
USPT 6 70 
USPT 680 
USPT 690 
USPT 700 
JJSPT 740 
USPT 750 
USPT 7 60 
USPT 770 
USPT 780 
USPT 790 
USPT 800 
USPT 810 
USPT 820 
USPT 830 


u o op u,o o'm up u.u o'u u u uju o o 


SUBROUTINE PPIKR(S,T, ISAVE,NPEAK,NBANOStS IGLEV) 
dlHENSlON B( I > .T( n, I SAVE( II 


* « * « « «« « * 

SUBROUTINE P P I K R 

PPIKR RETURNS THE VECT 
B- INPUT VECTOR OF AMPL 
T- INPUT VECTOR OF f-VA 
ISAVE-OUTPUT VECTOR OF 
NPEAK-OUTPUT SCALAR - 
NBANOS- INPUT SCALAR - 
SIGLEV- INPUT SCALAR - 
EXCEED IN ORDER 




OR INDEX OF ALL SIGNIFICANT SPECTRUM PEAKS** 

ITUOE SPECTRUM EST I MAT E S * 

LUE"S~FdR EACH SPECTRUM ESTIMATE ♦ 

INDEX numbers of SIGNIFICANT PEAKS 

NUMBER OF PEAKS FOUND • 

NUMBER OF BANOS IN SPECTRUM „ 

CRITERION WHICH T-VALUE FOR A PEAK MUST * 

TO BE RETAINED 


*********** ****’********************** 


-SEARCH FOR PEAKS* RETAIN 'SIGNIFICANT* ONES *„ BU ItD IS AVE*^ 


NPEAK^O 

DO 10 L=2*NBANOS 
IF (B(L-U-fl(L)> 
20 IF OIL^l )'3(L) ) 
AO IF I T( L l-SIGLEVI 
50 NPEAK-NPEAK+1 
ISAV6INPEAKI-L 

To contTnue 

RETURN __ 

'END 


20 * 10.10 
AOf 10, 10 
10.50,50 



xviii 


CHIT 

CRtT 

FUNCTION CRITR(RS0*T,JCRITR,U > CRIT 

DIMENSION Tt I ) 

V* ♦ ♦ ♦ ♦ ♦ **~i * * ♦ ♦•♦♦♦♦ ~C R 1 1 

^ _ _ * CRIT 

♦ FUNCTION C R I T R • 

^ ♦„ CR I T 

♦ critr selects a criterion for the spectrum OPTIHIZATIQN by peak • CRIT 

♦ SHIFTING. SELECTS EITHER R-SOUARE OR THE APP RO PRIATE T VA LUE* 

______ * CRIT 

♦ . CRI T 

- ■ ;■■ " ■ ^ CRIT 

'SELECT APPROPRIATE VALUE FOR CRITERION# 5^!! 

'rSQ-INPUT SCALAR - VALUE OF R^SOUARE r^tr 

T- INPUT VECTOR - T-VALU ES 

JCRITR-CONTROL OKTiT for selection of either RSQ or TIL) AS CRITERIONCRIT 

0= R SQ AS C« I TER I ON . . J 

"UT-VALUE AS CRITERION CRIT 

L- INDEX VALUE OF V ECTOR _T CRIT 


CRIT 
CRI T 
CRIT 
CRIT 

CRIT 


IF IJCRITR) 
10 CRITR = RS“0 
GO TO 30 
20 CR1TR=T (L) 
30 RETURN 
END 


10 « 10«20 




n no!r» ! ooo i ? 


KSPJ 

■£- - KSPT 

SUBRDUT INE ASPECT (R, T,N,M*NPEAK* I SAVtRX.B*JPRINT ,XBAR, STO,OtRSO* KSPT 
ftANDS*84N> 

OIMENSmN iSAVEdOOt fR(l)»TCl>tRX(l)tBtH«XBA RCUt STPX1J_.*QJ A 1*— 

IPHASEt IJ , BANOS I I ) . BANrU ,ANS( 10) * ISAVt U ,RY(100) 

COMMON K RE ADR ♦ K PR I NT 

-CONVERT ISAV TO ISAVE* FROM FREQ BAND VECTOR TO VECTOR OF SIN/COS 

J521 FORMAT (•0JKL=*»I5) 



OOH^O L«l.NPEAK 
[SAVE(K)«ISAV(L)*2 

K*K^l . llll 

100 ISAVE<K)-ISAVE(K^l)»l 

-SELECT SUS-'MATR IX OF PREOICTORS ETC FROM R# ACCORDING TO IS AVE- 


call OROERIM.R tl ,K 


, ISAVE, RX*RY) 


-INVERT K-ORDER MATRIX OF PRED I CTOR" INTERCORRELAtf RX 

CALL MlWfRx7ir;^Tf78 ,T) 

IF IJPRINT-3) 300.A00,400 ^ 

“Abo WRITE (KPRINT,7) OET 

7 FORMAT! ‘O*,/, '0 DETERMINANT a*,E3Q. 15) 

-COMPUTE REGRES SICNS(SP ECTRUM) FOR NPEAK SIZED MODEL# ^ ^ 

300 CALL MULTR ( N, K , X8AR * S TO, D,RX, RV# J S A V JE , 8 B ,JT , 

RSO=ANS( 2)**2 

-CONVERT SPECTRUM TO AMPl I TUDES BY COMBI NING S 0, 
-ALSO CONVERT T« 

NF RE 0» 0 ' 

DO 200 L^* I* K, 2 

NFREQ=NFREQ+l — 

PHASE! NFREQI =ATANIB( U /BIL^in 

B!NFREO»=SORT!B(U ** 2_. L»U»*2) ■ 

fob" t( NFREQIsSQRT! t (L)**2 ♦ td+l)**2l 

10 RETURN . ^ 

END 


OBIG^ 



;o o n o o n r» o n n o rJo olo o o o o r> n r«o n o n o o n n o o! | ri;f» 


__ MULT 10 

■ ■ "■ ' mult 20 

SUBRaUTINE MULTK IS4VE,NPeAK*N8ANOS,R,T«N*M,RX,B*XBAR,STO#0*RSO# MULT 30 
I PHASE* BANDS » BAN) 

DIMENSION ISAVtU)*R(UfT< XBAkU )* ST0(U *D(i )t PHASE M ULT 5 

I ( I I,8An 6V( i rrSANri) ,KSAVf< 20) tSI 201 ,TSPEC(20> *TPHI26) MULT 50 

COMMON KR6A0R*KPRINT ^ ...MULT 70 

MULT 8 0 

♦ »♦♦»»»»»»»♦»♦»»»»»»♦♦»»♦»»»»*♦♦♦♦»»»»»* »♦»»»♦♦»♦»» ♦»♦»»♦♦»»♦♦♦*♦*»♦*_■ mult 90 

- - - - -- --- ^ hult'ioo 

♦ SUBROUTINE MULTI ■ ♦ MULT 1 10 

*MULT 120 

♦ MULTI COMPUTES THE FULL SPECTRUM OF A TIME SERIES USING A K^l ♦,MULT J30 

♦ BANOS ' MODEL WHERE K=THE NUMBER OF’ PEAKS IN THE CPTQMIZED ^ MULT lAO 

♦ MODEL. THE EXTRA BAND IS THE WAVELENGTH IN QUESTION* THE OTHERS * MULT 150_ 

♦ ARE THE SIGNIFICANT PEAKS. THE VARIANCE ACCOUNTED FOR BY ♦ MULT 160 

♦ THE SIGNIFICANT PEAKS IS THUS ^ACCOUNTE D FOR* IN EA C H ESTIMATE. ♦ MULT 1 70 

♦ WHEN the current BAND IS ONE OF THE SIGNIFICANT PEAKS* A SKIP ♦ MULT 180 

♦ OCCURS. VARIABLES DEFINED IN MAJNL I NE._ 190 

"♦ " • MULT 200 

♦»»»»»♦♦»♦»♦♦♦♦♦♦♦»♦♦♦♦♦♦»»♦»♦♦*»♦♦»♦»»»♦»♦♦♦♦»»»»♦*♦» *♦♦♦ »♦»»»♦»♦**♦ MULT 210 

- - mult 220 

MULT 230 

-LOOP THRU ONCE FOR EACH FREQ SANO^ NFMOO IS THE NUMBER OF SIGNiriC^ MUlT^AO 
“CANT PEAKS IN MODEL PLUS A CURRENT BAND. ' MULT 250 

■' " “ ; MULT 260 

NFMOO»NPEAK<> 1 , ' MULT 270 , 

00 310 LF=l.N8AN0S MULT 280 

. mult 2 90 ^ 

-CONSTRUCT NEW I SAVE VECTOR ^ L T 30 0 

M ULT 310 

CALL BUlLblkSAVEVrSAVEiLF rNFMOOflSKI^^^ MULT 320 

MULT 330 

-SKIP I F THIS BAND IS IN THE MODEL ALREADY ' “ MUL T 340 

MULT 350 

1 F “(TSlTiP J 32"0, 320* 330 ^ULT 36 O' 

' MULT 370 


BRANCH TO 330 IMPLIES THA T THE CURRENT BAND' fS ONE OF THE' SI GNIF ICANT MULT 380 
PEAKS. 00 AN NPEAK SPECTRUM AND PICK OUT CURRENT BAND. _ _ MULT 390 

MULT 400 

330 CALLKSPECTIR *S*N*M ,NPEAK* 1 SAVE * RX, TSPEC * 0* XBAR, ST0*0,RSQ*T PH, MUL T 4 1 0 

IBANOS.B'AN') ' MULT 42 0 

MULT 430 

TSPEC CONTAINS ALL MODEL EST I MATES . ' P ICK THE ONE' FOR LF. MULT ^40 

MULT 450 

00 331 LLZ=1,NPFAK ' ' ^ " - - -- - 

1_F (LF-ISAVE(LLZ )) 331*333, 331 MULT 47 0 

333 6( L'FT = TSPEC*( LLZV MULT 480 

PHASE ( LF J = TPH(LLZ) J _ .... MULT 490 

T(LF)=5(LLZ) “ “ " " ' MULT 500 

331 CONTINUE „ . MULT 510 

GO TO 310 MULT 520 

C _ _MULT 530 

"C^BR ANCH"'TO"326“'fMPLrES~ THIS IS A FU'LL NFMOD MODEL. DO SPECTRUM MUl't 540 




xxi 



C-ANO 


PICK OUT CURRENT BAND, 


320 CALLKSPECT{R*S.N#H* N F MOO f KS A V E t ® • J * 

I BANOS* SAN » 

00 321 LLL*l,NFMOO 
IF(LF-KSAVE(LLL J I 
TILF)=StLLL ) 


32 I *322 •321 


322 


STun =T S PEC ( L L t ) 
PHASE (LF)*TPH(LLL» 
"321 CONTINUE 

JIO CONTINUE 

RETURN 
END 


MUUT 

MULT 

MULT 

MULT 

MULT 

MULT 


MULT 

MULT 

MULT 

_MULT 

MULT 

Mill T 


550 

5S0 

570 

5B0 

590 

600 

610 


620 

630 

640 

650 

660 



xxii 


SUBROUTINE BU IL D(KSAVE» I SAVE *1. FtNFMOO # I SKIP* NPE AK > 

DIMENSION KSAVEd I »1SAVE( 1) ,ASAVE(30) * RANG I 30 1 

"*■*♦*♦♦♦•♦♦♦♦♦*♦♦♦♦♦♦♦*♦*♦♦♦♦♦♦*♦♦♦♦♦*****♦ *****’****^******^^***** 

♦ BUILD CONSTRUCTS VECTOR OF RANKED FREQ 8ANp_NUMB6RS .G^V^ 

“i particular model and a current frequency* 

* KSAVE - OUTPUT VECT OR TO BE B UILT 




BILD 

BILO 

aiLO 

BILO 

BILD 

BILO 

BILO 

BILO 

_BILD 

BILO 

BILO 


C • I SAVE - iNPlif VECTOR CONTAINING NFREO SIGNIF PEAKS 
C ♦ LF - INPUT SCALAR ~ NFW FREQ BAND TO BE ADDED I CURRENT ONE I, _ 
___ ^fjpMOO - INPUT ~ NUMBER OF FREQ BANOS IN NEW MODEL 

C ♦ ISKIP - OUTPUT CONTROL DIGIT TO KEEP FROM COMPUTING ESTIMATE 

OF ONE OF THE OLD PEAKS AGAIN 


NP6AK ~ INPUT - NUMBER OF FREQ BANOS IN OLD MODEL, 


BILO 

BILO 

BILO 

BILO 

BILO 

BILO 


10 

20 

30 

40 

_50_ 

60 

70 

80 

90 

100 

Up__ 

120 

130 

140 

150 

160 

iTO_ 

iao 

190 

200 

210 

220 

230_ 

240 

250 

260 

270 

280 

300 

310 

320 

330 

340 

360 

370 

330 

390 

400 

420 

430 

440 






SEE EAOY IN MO OEL 

00 10 L-lrNPEAK 


IF(LF - ISAVEILU 10*99*10 
10 CONTINUE 

■STORE VECTOR OF ISAVE PLU S LF I NTO A SA VE AND RANK IT* RAW STOREO__X?* 

VECTOR RANG 


DO 20 L=l#NPEAK 

20 ASAVEI L1 = ISAVE< L» 

AS AVE ( NF MOO I L F 

CALL RANK(ASAVe*RANC * NF M 0 0 > 

c 

C-BUILD KSAVE OF RANKED FREQ BANDS 


00 30 L=l»NFMOO 
“KOOL=RANC(D 
30 KSAVE! KOQL) =ASAVE( LI 
I SK I P=0 

GO TO 100 


99 ISKIP=1 
100 RETURN 
END 


BILO 

BILO 

BILO 

BUO 

BILO 

BILD 

BILD 

_BILO 

BILD 

BILD 

BILD 

BILD 


BILO 

BILO 

SILO 

BILD 

'aiLO 

B!LO 

BILO 

BILO 

BILD 

BILO 

BILD 

BILO 

'BILD 

eiLo 

BILD 


xxlii 



SUBftOUTI NE OUTPUT < BANOS *8 AN I StT t PHASE i NBANOS^ PHNAME* 

1BANAME,XNAME, JPRINT, JPL0T,JDQMANI - - - 

DIMENSION BANDS(l) *B^N(l » « B (l ) f T ( U «PHASE U I ^ AMNAHE < H f TNAMEI H t 
Tl>'H N AM p < I > ♦ B A N A H E ( 1 T» X NA ME“( I'l 

COMMON KREA0R,KPR1NT ^ . 


c 
c 
c 
c 

C ♦ OUTPUT PRINTS ANO/OR PLOTS SPECTRA WITH OOCUHENTAT |ON« 


• •#***« ^* «««*«««***««*««« «♦»♦»»»♦» « «««*♦* 4i« ••♦♦♦♦♦ • 

■ ~ "■ ■ ■" '■ ~ ■" ■ 

♦ SUBROUTINE OUTPUT ■ ' • 


^^**mt:^^*m^**m**^m*******m***************** ************************** 


C 

C-PRINT OUT DOCUMENTATION AND VE CTOR OROER IN APPROPRIATE DQMAINI 
t-FINO”^I^LE values ANO PLOT 



OUTP 

OUT? 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 


TOPICS. 28318 

BTOPI=-TOPI 

"BMAXsBin 

DO 64 LBAN0=2f N6AN0S 


_68 

64 


IF IBMAX-BILBANOn 68t64t64 

BMAX*B(LBANO » 

CONTINUE 

TMAX=T( I ) 

DO 66 LBANO*2t N 8 A NO S 
IF (TMAX-TILBANDI) 67»66«66 


67 TMAX= T ( L bXn 0 1 
66 CONTINUE . 

IF ( JOOMANI 60 t 60 f 50 


C-REV6RSE 

C 


VECTORS IF FREQ DOMAIN IS USED - PRINT - REVERSE AGAIN. 


OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

-OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 


50 CALL REVIB.NeANOS) 

CALL REVIT.NBANOS) - . - ^ ' '■ 

“CALL REVIPHASE.NBANDS) 

IF (JPRINT) 52.52.51 „ _ . , 

51 WRITE IKPRINT.2I 

2 FORMAT! »OFREQUENCY * .JX.* AMPLIT UDE*. 5X.*T-VALUE« . 8X. ‘PHASE* ./.* 0* 
00 55' LDAND=l . NBANOS 

55 WRITE (KPRINT.3I BAN( LBANOI , B (LBANO) .T I LBANOTtPHASEILBANO) _ , _ 

3 F0RMATIF10.3.5X,F9.3*5X,F7.3.5X.F8.3» 

52 IF IJPLOT) 54,54,53 

'53 CALL PLDT2(tt,8MAX,O.0.T,tMAX,0.0,N0ANOS.BAN,AMNAHE,TNAME.BANAME I 
CALL P 10 T?( PHASE, TOP I.eTnPl.T.TMAX.Q.O.NBANOS.BANfPHNAME.TNAME. 


1 BANAME I 

54 CALL REV! B,NBANOSI 
"■ CALL REV(T ,NBANOS) 

CALL REV(PHASE,NBANDS) 
GO TO 100 

60 IF I JPRI NTL,;,62 ,62 .61 

61 WRTTE <KPRINT,4J 


OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

iOUJP_ 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

_J)U.TP_ 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

_pyj^ 

OUTP 


100 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
>70 
480 
49 0 
500 
510 
520 
530 
540 


ORIGINAllBAC^Ji 
OF POOR jQUAUnO 

rv^PRECk)n}G PAGE BLANK NOT FILMED 

t f ^ 
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% 



FQRM4T( 'OPER inO 
65 


• ♦ 5X • * A MPL I TUDE * »5X, *T-VALU£ * » 8X « « PHAS E » f / > * P *_ 

"00 65 1 ft AN 0=1* N B ANOS ^ o ^ . o . k. 

65 WRITE (KPR!NTt3J BAND S < LB AND I t 8 (L 8AN0 1 » T ( L 8 AMO ) » P HASE ( L BAND > 

63 CALL PL0T2 lB*8MAXf0.0iT*THAX.0«0*N8AMOS» BANOS* AMNAME, TNAME * XNAHE I 
CALL PLOT 2 (PHASE* TOP I ,B TOP I . T * TMAX , 0. 0 , NBANOS * BANDS tPHNAME , T NAME* 

__ I XNAME) — — — — 

1 00 RE' turn 


lOUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OJTt* 

_OUTP 

OUTP 

OUTP 


550 

560 

570 

580 

590 

600 

610 

620 

630 










XXX 


u u; 


SUBROUTINE PIOT2 I Yl , YMAXl ,YMlNlfY2* YM&X2 , YH IN2 tNP , XAX »YNAH It 
I YNAM2,XN4H) 

D I ME NS [ ON Y I ( I I # Y 2 ( I I . L I N EJ .1 \ 1 ) < l it XNAMf 5 >t YNAM ir9)#YNAM2C9> 

COMMOY' Kqf APR,KPRNT 

DATA KBLNK/* */ .KBOR 0/ •.*/• KAL /*♦*/* KPLOT /*♦' / 

PL0T2 GENERATES TWO DVER/UNDER GRAPHS PER PAGE ANO CALIBRATES 

THEM. Yl^VERTlCAL VARIABLE l» VHAXl, YM INl= VERT ICAL PLOT LIMITS. 
YMAX2. YMIN2=S1MHAR FOR SECON VARIABLE* NP*NUHBER OF POINTS* 


— Y2* 


PLT2 

10 

PLTZ 

20 

PLT2 

30 

PLT2 

40 

PLT2 

50 

PLT 2 

60 

PLTZ 

TO 

PLT2 

80 

PLTZ 

90 

PLTZ 

100 

PLTZ 

110 

PLT 2 

12 0 

PLT2 

130 


C — 


XAX = X' AX I s’ {COMMON TO BOTH VARIABLES). 

WRITE (KPR^4T,2) YNAM 1 * YNAM2 # XNAM 

2 FORMATIM* / T5.**V T5 . • * * • T25 , • * • *TBO, ' * • / T5 *»♦*, T 25 *•♦**' PL T2 

l*T80t***^*/ T 3 *'♦♦♦♦*** T I 8 ** ♦♦**♦*♦’***** * f9A4* T73* * ♦♦♦*♦♦*♦♦*♦* *PLT2 
2 * 9AA/TA *•♦♦♦** T 25 * *♦♦♦** T80 * * *♦* * /T5* • ♦ * # T25 * • ♦ • *T80* * ♦ * / IX * 5A4 ) PL T 2 

PLT_2 

■ ' PL T 2 

PLT2 

■ ■ ' PLT2 


C-- 


IT. 




GENERATE SCALE FACTORS 

SKALl= SO. 5/C YMAXl-YMI Nl ) 

SKAL2=50.5/ (YMAX2-YM IN2) 

GENERATE CALIBRATION AND OUTPUT 
CALOl = YMINl 

CAL02 = YMIN2 

CALI ='YmAXI 

CAL2 - YMAX2 

CAL 51= YMAXl- I . 5* ( YMAXl-YM INI ) ) 

CAL52=YMAX2- (. 5*( VMA X2-YMIN2I ) 

WRITE (KPRNT*ll CALO I »CAL 5UCAL l*CAL02*CAL52*CAL2 
1 FORMAT <♦ « .F14.4.F?4.4.F23.4«F12. 4»F23.4,F24.4) 

generate VERTICAL AXIS AND OUTPUT it 

00 10 L = 1,111 

LINECL ) = KBORO 

LINE(l) = XAL 

LINEC26) = KAL 
LINE(5l) = 

LI NEC 61) = 

UNEC86) = 

L INF( 1 1 n = 

00 20 L=52,60 
LINEIL.) = KBLNK 
WR ITEC KPRNT.4) 


PLT2 

PLT2 

J*tT2- 

PLT2 

PLT2 

PLT2 

PLT2 

PLT2 

PLT2 


10 


X AL 
XAL 
XAL 
XAL 


20 


L INE 

FORMAT! lOX, U lAl) 


DECIDE HOW MANY LINES BETWEEN POINTS 
NWloe = 32/NP 


LOOP FOR 
00 30 


number 
LP = I 


OF Pq I NT Sj^ 
NP 


EACH TIME plotting A POINT AN0_NWI0E_ 


-CLEAR 

nn 


LINE 
^0 t 


ANO 

- I 


GENERATE 
\ I 1 


GRID POINTS 


PLT2 
PLT2 
PLT2 
PLT2 
PLT2 
__ „^PLT2 
PLT2 
PLT2 
PLT2 
PLT2 
PLT2 

PLT2 

PLT2 
PLT2 
PLT2 
PLT2 
PLT2 
SPACEPL T2 
PLT2 
PLT2 
PLT2 
PLT? 


140 
150 
160 
170. 
180 
190 
200 
210 
220 
_230_ 
240 
250 
260 
270 
280 
.290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
39 0 
400 
410 
420 
430 
440 
450 
460 
470 
48 0 
490 
500 
510 


40 LlNt(L) = Kal.NK 
__ LINE Cn=KBORn 
’’ L I NEC 26 )=KBnRO 


PLT 2 ^ 
PLT2 530 
PLT2 540 


xxxi 


LI NEf51 )=KBORO — 

U NP( 61 ) =KBORO 

HNE{ 86 >=KBQRD ^ 

LIN6U 1 l)=KBnRD 

C . . - . 

C — TEST FOR OVER Gft UNDER RANGE AND TRUNCATE OUT-OF-BOUNDS VALUES 

IF (YMLPI-YMAXll 2000t 2000> 2100 ' 

2160 YA =" YMAXl 

GO TO 5100 .. ^ , 

2000 tF4YULP)-YMINl J5000* 5200* 5200 

5000 YA = YHINl . - ^ 

GO TO 5100 

5200 YA = Yl(LPt -- ^ 

5100 IF (Y2(LPr-YMAX2» 200. 200* 210 

210 Y0 = YMAX2 __ , 

“ GO TD 510 

_ 200 IF (Y2 I LP)-YMIN2) 500, 520, 520 _ 

500 Y8 = YMIN2 

GO TO 510 - 

520 y”B = Y2<"lpT 

~ C — NOW GENERATE POINT INDEXES 

510 Kl = IVA -YM INU^SKAL I ♦! ,0 

K2 = I YD -YM IN2)*SKAL2«-61,0 

UNF(KII= KPLOT ■ 

LINE(X2)=" KPLOT 

C— OUTPUT HORIZONTAL SCALE VALUE AND LINE 

WRITE (KPRNT ,31 XAX(LPI,UNE 

3 FORMAT (IX, F 8. A. IX, 1 1 lA 1 ) 

C „ _ 

C--CrEAR LINE and GENERATE BLANK LINES FOR SPACE 

IF (NWIOE) 30,30, 110 . 1 : „ .. . 

no 00 95 L = l,ni 
95 L INE( L »=KBLNK 
LINE (1 l=KBnRD 

LI NE( 51 > =K80Rn 

L INF( 26 )=KHnH'J 

LI NE (6l» =KBGRO 

LI NE< 86) 'KBORO 

LI NEC I I I ) =Kiu:)Rr» 

c 

C — OUTPUT PLAIN LINE 

DO 90 LOnP=l,\WIDE 

90 WRITE (KPRNT.A) LINE 
30 CONTINUE 

' ■ C 

C — GENERATE RIGHT BORDER AND OUTPUT 


50'LINF(L) ^ KrtClRD 
L INE( 1 ) = KAL 
L : " f ' ■ 

LI Nht Ml = -<AL 
LINF(61J =^,AL 
LINE (86) ="k'aL 
" LI NE(1 I 1) = KAL 


PLT2 550 
PLT2 560 
PLT2 570 
PLT2 580 
PLT2 590 
PLT2 600 
y_T2_61 0 
PLT2 620 
PLT2 630 
PLT2 660 
PLT2 650 
PLT2 660 
Ai-X2 67 0 
PLT2 680 
PLT2 690 
PLT2 700 
PLT2 710 
PLT2 720 
J»Lt2__730 
PLT2 740 
PLT2 750 
PLT2 760 
PLT2 770 
PLT2 780 
_PLT2__790„ 
PLT2 800 
PLT2 810 
PUT2 820 
PLT2 830 
PLT2 840 
PLT2 650 
''pl ’T2 '86 0 
PLT2 870 
PLT2 880 
PLT2 890 
PLT2 900 
PLT2 910 
PLT2 920 
PLT2 930 
PLT2 940 
PLT2 950 
PLT2 960 
PLT2 970 
PLT2 9H0 
PLT2 990 
PLT21000 
PLT210I 0 
PLT21020 

i I. / L 

"PLT21040 

PLT21050 

O • ^ t / r\ 

PLT2I0B0 
PL T21090 
PLT21100 


xxxii 

ORIGINAL PAGB5 IS 
OE* POOR gUAU^. 




00 80 L= 52,60 

WHNE TO = K8LNK ' 

WRITE (KPRNT.4I ONE 

RETUR'4 

END 


PLT21U0 

>021120 

PLT2H30 

PLT2U^0 

PLT2U50 




XXXllX 


non on no noon n o o n o o o n o o o o o o o o o no 


SUSROUTTNE T OT A M ( X , N , NPE AK ♦ IS AVE*M* PHASE • Bt AMP I 
DIMENSION XI 1) ,1 SAVE I I J »P HA SE II » , B ! U 


« ♦'*♦**♦♦♦**♦♦ ♦♦♦♦♦♦♦ 

* _ • 

♦ SUBROUTINE T 0 T A H ♦ 

♦ * 

♦ TOfAM COMPUTES THE RESULTANT AMPLITUDE OF THE SUM OF NPEAK SIN- ♦ 

♦ soios WHICH have their amplitudes and phases oetreminep by the * 

“♦ CALLING PROGRAM.’. ~ TOTAH simply SELECTS TH^ APP'rOPR VaTE SINE ♦ 

♦ AND COSINE WAVES FROM X* WEIGHTS THEM, SUMS THEM AND THEN COM- ♦ 

♦ Pures THE amplitude of the sum. * 

♦ WARNING THIS WILL DESTROY ROW I OF X WHICH ♦ 

♦ CONTAINS THE ORIGINAL DEPENDENT VARIABLE. 


♦ ' * 

-LOOP’THRU THE FOLLOWING NPEAK TIMES TO CREATE THE RESULTANT WAVE INX 
DO 100 LP=1. NPEAK 

K=tSAVE(LPI ^ 


KS=KCf I 


■THE FOLLOWING COMPUTES THE PROPER AMPLITUDES FOR THE COSINE AND SINE 
■COMPONENTS, BASED ON PHASE AND AMPLITUDE ESTIMATES. 


PC-l.O 

PC=PC/ITANIPHASEIK>U PCI 
PS=PC*TAN(PHASe (Kl ) 


SELECT SINE AND COSINO VALUES FROM X A 

SELECT SINE AND COSINE VALUES FROM X, SUM AND PUT INTO X< 


ROW 1. 


100 


DO 100 L=l,N 
CALL LOCIKC ,L , J ,M,N,0 I 
CALL LOCIKS.L , I ,M,N.O» 
X(L) = X(JI*PC ♦■ XI II*PS 


-COMPUTE APLL I TODE OF RESULTANT 


SX=0. 0 
00 10 L=1.N 
10 SX = SX«^XIU 

SX=SX/FinATI N ) 

'SXS=0.0 

00 20 L=lfN 

20 SXS^SXSM XU )-SX>**2 
AHP = SORT { SXS/FL. ' i . 
RETURN 
ENO 


TOTM 

TOTM 

TOTM 

TOTM 

TOTH. 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

'totm' 

TOTM 

TOTM 

TOTM 

TOTM 

_TOTH 

fof M 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

jrpTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 


10 

20 

30 

40 

^0 

hO 

70 

80 

90 

too 

no 

120 

130 

140 

ISO 

160 

1.70 

lao 

L90 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
36 0 
370 
330 
39 0 
400 
410 
420 
430 
44 0 
450 
460 
47 0 
43 0 
4«)0 
500 

520 

530 
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USE 



OF S P E C T 


THE FOILOWlN'; IS ^ LIST OF INPUT PAR &METHRS AND OAT A CAROS AND THEIR 
INTFRPKtTAr I JNS. 

JDO N-0 A I N 1F ANALYSIS* 0 DR ttl ANK= PER I OOt 1 = FREQUENCY« 

jqANHS-Hig UF FRF'iUFNCY OP PFR lOD 3ANns_W.Tt.l TAI NEQ, 

""O [K \aLANK = G£NERATF ACrOPOlNG TO I I^^US SUBSEQUENTLY READ IN. 

l = PEAD IN A VICTOR OF BANDS Tu BE USED FROM CAROS. 

NBANDS-H.1^ many FREOJENCY or PERKU) BANOS WILL BF IJSEO. 

NOBS-THr N IFD^ER OF nASERVATIUNS IN THE DATA VECTOR 
JANALY-hI'-i TUt TATA SHALL BE ANALYZED. 

0 Q;^ _BLANK = nO JJN n/ARIArp SP ECTRUM 1 SEARCH F OR BEST J^VARJA 

' MODEL 

I = S IMPLY DO I.INIVARIATF SPFCTRUM 

^=00 K.-VAPIATF SPECTRUM ACCOROINO TO MODEL PARAMETER READ IN 
JPATA-TFILS SUrROUTlNF ()E PVAR WHE^E TO GET DATA. 

0 OR BLAN<=READ FROM CARDS ACCORDING TO FORMAT <10X.7F10.0I 

1 =GF N t R E data JFR g.*^*!*"'*^ L C Al L E 0 BY DATA,. T_H J_S 

MOOELINS PROGRAM WILL ASK FOR MORE CARDS. 

(j:^-c;pLtCrS CPITERION FOP OPTIMIZATION OF K-VARIATE SPECTRUM. 

0 OR JL ANK =MUL T 1 PL F R-SQUARED 

1 = T- VALUE S - — - - 

JPR I NT -PR I NTOUr fONTROL DIGIT 

0 =_0 0 ^ I r _P R I . N T_A, N Y_R F ,S U I T S 

i=PRlNT ONLY FINAL RESULTS 

2=PRINT = I NAL AND INTERMEDIATE RESULTS 

3=PRINT AS FUR 2 PLUS MATRICES. C ORR ELAT 1 ONS ♦ ETC FOR EACH STEP 
JPLOT-PLOT control digit* - values as above EXCEPT FOR 3. 

HAN OL i ) - L 0 C S T F R F 0 . 0 R L O N G F S T^ _P F P.l C D ,„A $_..AP .P PAIA T_E_J10„^^^^^^ I N . 

PA^ljfLLijfLT A( FRF0.l'"'j3^'lji“LTA(PFRldDI t AS APPROPRIATE TO DOMAIN. 

B ANPS-VfC y,)R of hands. 

NS I N-Nll Y.^r R '1F SINE WAVES IN SlMUt.ATFO DATA. 

!\g p _ r , } , t M F K 1 0 U OF Y H F (. S I N T H SI MU L ATE D . .S I N E W A V F 

Sd-S"^ AN'iA^v) DFA/IATION OF THE L S 1 NT H S I MUL AT ED SINE WAVE. (AMPLITUDE) 

SDN^S T A N ) A Rl) OF V . OF GAUSSIAN NOISE 





THF F(!t L'JW T N" IS A LIST OF INPUT CARDS AND THEIR CONTENTS 


CARO 1 ' - - 
rni UMN 

MANDATORY 
VAN r A lU F 

COMMENTS - - - 






l - ^ 

6-10 

JOOMAN 

JBANOS 

RIGHT JUSTIFY 
RIGHT justify 






11 - IS 

15 - 20 

NMANDS 
- IPS 

RIGHT justify 
RIGHT JUSTIFY 






21 - 2S 
26 - TO 
31 - 3S 
36 - 4 0 

JANALY 
MATA 
JCR T’tq 
JPRI.NF 

RIGHT JUSfIFY 
RIGHT JUSTIFY 
RIGHT justify 
RIGHT JUSTIFY 



■ ' 



- 4S 

- |r*1.0T 







CARO 2 - - - 
COLUMN 

PC JJ I RED ONLY IF JL3AN0S=C, 

vauiahle comments - - - 

OTHERW ISF 

OMIT. 



1 - 10 
LI ” 2"* 

tandlo 
- BANOEl 

PUNCH DECIMAL 
PUNCH DECIMAL 







CARDS 3 - - 

- liLl RED 

ONLY IF JBANDS=l. 






COLUMN 
1 - 10 
11 - 20 

VAR I A.MLH 
RANnS( 1 J 
- RAN0S(2) 

CU«MENTS - - - 
PUNCH DECIMAL 
PUNCH DECIMAL 

- 1ST 

- 2ND 

BAND 

BAND 

IN SPECTRUM. 
IN SPECTRUM. 




♦ _ _ . 

* 

CONT t iN!U^..„F_i P AS MANY CARDS AS NEEDED TO READ IN NBANDS BANDS. 

FORMAT ( 10* 0 I 

CAROS 4 - - - FTOIMRED TNUY IF J0ATA=0* 

COLUMN VA'MAfM. T COMMENTS - - - 

1 - 10 - l:) r(inE may CONTAIN ANYTHING - NOT READ. 

n - T O - •' T A { 1 ) 1S T d ata POI NT, PUNCH DECIMAL OR RIGHT JU $.T I F Y,., 

21 “ 3J - TATA(2> 2ND DATA POINT. PUNCH DECIMAL OR RIGHI JUSTIFY. 

♦ ■ , . 

♦ 

CONTfNUF <^A'-'E format ON CAROS AF NEEDED FOR NOBS DATA POINTS. 

FORMA! { LO K. 7F 1 0.0 I 

CAPOF S - - - R>-iUlRF:0 ONLY If JDATA^L. 

5A - - MA Jf/AT PY . _ ... 

COLUMN VAM A HL COMMENTS - - - 
I - :> - NS IN RIGHT JUSTIFY. 

>0 ;o! ■> ) ‘ONLY^.I r_..NS IN ..I S NfiT,...ji:,.F.Rn. ' 

COl JM'i - COMMFNrS - - - 

I - I ) - MANO(l) PUNCH DFCIMAL. 

11 - ^ - S)(l) PUNCH OFCIMAL. 

CONrp-.i,; ^.\m: FORMAr THRU NSIN CAPOS (BANDS). FORMAT ( 2F 10. 0 1 

- M \ -J!' A "^OR Y ^ .. 

Cf'L.JMN VA-'fABl.: COMMFNTS - - - 


1 - 

M - 

ON 

t’UMCJJ I'ECIF'AL 




\ - _ 

i >U ! 

.1 iiNl Y 

I SON IS Min 2 FRO 




COLUMN VA' 

I A''i r 

COMMENTS - - - 




1 - 

. 1 

■ 1 C li ^ F 

MAY C lUnAfW anything 

7 NOT PF AO. 


S - 

1 2 '1 

r ( 1 

) PUL; CM UK. 1 .'CAl . 

1ST 

GAUSSIAN RANDOM 

NUMjFR 

1 1 - 
* 

1 ;■ N 

) I S[ ( 2 

) PI.INLM nLC imal- 

2ND 

GAUSSIAN RANDOM 

rjOM )ER 

cc\ r ! 


J ^S M 

ANY CARDS AS NEEDED 

FOR 

NOBS RANDOM NUMBERS. 

f 'k* • ’ A 

. ( ♦ A t 1 

‘ 1“ . ) 

xxxvii 






C*^***’*^*****^*****^<‘ <*+^# '^4c^ ** *4i 4i^ ^#%^4i4ii^ #% 4(4 1^*« 44 

c 

C ■ ■■ — ^ 


DIMENSION BANDS! 501, BAN! 50 I . PEAKS I 50 ) ,PHASE I 50 1 ,AMNAME ( 9 ) , TNAME 
1(9! ,PH NAMEIQI ,BANAMF I 5),X NAME 15) »LVA«1 100), I$AVE( 100) 

C 

C-THE FOLLOWING DIMENSION MUST BE EQUAL TO OR GREATER THAN 
C- NDBS«'( N^IANDS*?*-! ) 

C ^ 


DIMENSION XI30000) 

„. C-THE FOLLOWING DIMENSIONS MUST BE EQUAL TO OR GREATER THAN NBAN0S»2«-l 
C 

DIMENSION X3 AR( 100) ,STO( 100) ,B( 100) tDC 100) ,TtlOO) 

C ■'■■' ■ ■ 

C-THE FOLLOWING DIMENSION MUST BE EQUAL TO OR GREATER THAN 

C- M*M, WHE^E M=2->N8ANDS«-1 
C 

DIMENSION RX(IOOOO) ' ‘ ' ‘ 


C 

C-THE FOLLOWING DIMENSION MUST BE EQUAL TO OR GREATER THAN (M«-l)*H/2 

DIMENSION R(5050) 

COMMON KR6 ADR, KPRI NT 

DATA AMNAME/ • AMPL • , * I TUD * , • E ES*, • I MAT* « •E, •"V '* • • , • 

l»*» •/ 

DATA T NAME7'*'T - VA • , • L UE. «T* *T* ^T* *T^ 

•/ 

DATA PHNAME/* PHAS* ,• E E S *,• T I HA •,* T E. • ' *,• i,i 

I* , • •/ 

DATA 8ANAME/ »FREO* , »UENC* , * Y »,* t,* •/ 

DATA XNAME/ • PERI » *• OD •/ 

KREA0R = 5 — — 

KPRINT=6 


c ■- ■ 

C-REAO IN CONTROL DIGITS AS FOLLOWS AND PRINT THEM OUT, 

C JOOMAN-DOMAIN OF ANALYSIS. 0 OR BL ANK = PER I OD, 1 = FRE0UENCY. 

VECTOR OF FREQUENCY OR PFftIQO BANDS WILL BE OBTAINEO. 

C 0 OR PL AN K = GENE R ATE ACCORDING TO "L I M I TS SUBSEOUENTL Y READ IN. 

C 1=READ IN A VECTOR OF BANDS TO BE USED FROM CAROS. 

C N8ANDS-H3W M.\NY FREQUENCY OR PERIOD BANOS WILL BE USED. 

C NOBS-THE N!ivr.ER OF OBSERVATIONS INI THE DATA VECTOR 
C JANALY-MJW THE DATA SHALL DF ANALYZED. 

c 0 OR BLANK=On UNIVARIATE SPECTRUM, SEARCH FOR BEST K-VARIATE 

C MODEL ■ 

C l=SIMPLY DO UNIVARIATE SPECTRUM 

C 2=00 K-VARIATE SPECTRUM ACCORDING TO MODEL PARAMETER READ IN 

C JDATA-TELl 5 SUHRuuTINE DEPYAP. WHtRE TO GET DATA, 

C 0 OR PLANK-RlAD from CARDS ACCORDING TO FORMAT (10X,7FI0.0J 

C ^i = GfNrRArE :>\ra erom modeling program called by data, this 

'C NDDELING PRilGRAM WILL ASK FOR MORE CAROS. 


SPT 

SPT 

SPT 

SPT 

_SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

JPT^ 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

spr 

SPT 

SPT 

SPT 

SPT 




10 
20 
30 
50 
_ 50 
60 
70 
80 
90 
100 
_110 
120 
130 
150 
150 
160 
_IT0 
180 
190 
200 
210 
220 
_230 
250 
250 
260 
270 
280 
290 
300 
310 
320 
330 
35 0 
350 
360 
370 
38 0 
390 
500 
510 
52 0 
530 
550 
550 
56 0 
5 70 
580 
59 0 
500 
510 
520 
5 30 
550 


xxxviii 
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JCRITR-SEL tCTS CRITERION FDR OPTIMIZ ATION OF K-VARIATE SPECTRUH, 

3 r^L ANK = MULT IPLE R-'SdUAREO 

I = f - V « L ' ^ c _ 

JPRINT-PR INTOUT CONTROL DIGIT 

0=nnNT PRINT ANY RESULTS .. 

l=PR|NT ONLY FINAL RESULTS 

2=PRINT FINAL AND INTERNED IATE_RE SULTS 

^ 3=PR[,\ir AS for'?' plus' MATRi^C^^^^ COR'ReL ATI ONS t ETC FOR EACH STEP 

JPLOT-PLDT CGrTROL DIGIT. - VALUES AS ABOVE EXCEPT FOR 3. 

READ I KREAOR 1 1 » JDOMAN, J BANOS * NBANDSt NOBS. J ANAL Y . JOATA. JCR ITR, 

I JPRINT.JPLOT 

I FORMAT 

WRITE {KPRINT. 103) 

103 FORMAT Cl program PARAMETERS ♦♦♦♦*♦ 

WR I TE ( KPR I NT , I ) JOOMAN. J 6 ANDS.NBANOS.no BS.J ANAL Y f JOATA, JCR ITR. 

I JPRINT.JPLDT 

KTOP =MAX NUMBER OF T ! ME S THE PROGRAM WILL BE A LLOWE D TO TR Y IMPRO_VE^ 
■MENTS ON THE WHOLE MULflVARIATE SPECTRUM. KOUNt IS THE NUMBER OF 
■SUCHTRIES. _ 

KTOP=^ , 

KOUNT=0 

“SET irrGLE^^^^ FOR THE CRITERION OF ACCEPTANCE OF A PEAK. 

-THEN DECIDE HOW TO GET VECTOR OF BANDS. . . 


SIGLEV=2.0 
IF (JBANOS) 


100, 100, llO 


-REA D "I N ' L I H I T S 0 F " B A N 0 S "V EC f OR IF J B AND S WAS 0 - GENERATE BANOS. 
'OTHERWISE GO T3 UO AND READ IN BANDS FROM CARDS. 

BANDLO-LOWtST FRFQ. UR LONGEST PFR lOD, AS APPROPRIATE TO DOMAIN. 
BANDEL-D6L TA{ FREU. ) OR DELT A I PE R 100) , AS APPROPRIATE TO DOMAIN. 
BANDS^VECTHR OF BANDS. 

^ '1 00' RE AO ( K R F A 0 R 7 2 ) ■ ‘ BAN 0 Cd78 ANDEL 
2 FORMAT (3F13.0) 

BANDS! I ) ^BANOLU 
DO 120 L rTAN0 = 2,NBANDS 
120 BANOS! L band )-BANOS< LBAND' l ) ♦ 0ANOEL 

GO TtJ 1 in 

no REAL) !KPtA0R,2) ■!BANDS(LBAND),L8VND = 1,N8ANDS) 

'CUNYERT bands to period domain if ORIGINATED IN FREO DOMAIN, 

-REYFRSE Tflf PPOFR i)F BANDS SO AS TO BE IN ASCENDING WAVELENGTH ORDER. 
-ALSO GFNERArt BAN, THE FREO DOMAIN VECTOR OF BANOS. 

I 30' I F ! J 0 OMAN) I S 0 , 1 50 ♦ 1 7o ” 

I AO DO lAl I B ^ I , NBANDS 
0 AN ! L !i A N .1 ) = 3 A NO S ! L 0 A T4D ) 
lAl BANOS! ! " AND > = I . O/BANOS (L BAND) 

CALL RFV( BANOS, N9ANDS) 

C-GENERATESIN/COS pREOICTO R' M a t R Tx ■ A S ”P A~R T“ 0 F~ 0 AT A “M AT R I Y, ' X . THEN 


XXX ix 


C-Re4D IN DHPFNOENT VARI4BL E_P A_R T 01^ X_* THE S E DATA W I L L THE N BE ENT- SPT 

C*ER£D INTj SUBROUTINE COftftE TO 'COMPUTE THE C’ORRELATI ON HATR I X OF ALL SPT 

C-VARIABLES. _ _ _ _ _ 

C •" SPT 

150 M^NDANDS*? 4- 1 _ SPT 

CALL PREGENIX .RANDS, NOBS^NBANOSI ' SPT 

CALL r»EPVAR( X.NURS, JOATA, JDOMAN.JPRINTI SPT_ 

C ■ ■■ SPT 

C“GENERATE VECTOR OF VARIABLE NUMBERS FOR LATER USE. $PT 

C . - 

00 155 LV = l,M , SPT 

155 LVARI L V J =L V ‘ ' ' ' SPT 

CALL CORRE (NOBS ,H, I, X,XBAR,ST0*RX,R,B,0,T» SPT 

C " ■ ■ SPT 

C-SET TOPP, THE MAXIMUM ALLOWABLE TOTAL AMPLITUDE IN THE FINAL SPECTRUM SPT 
C-TOPP IS THE TOTAL AMPLITUDE OF THE FIT)4-25 PERCENT. SPT 

C , _ _ _ _ SPT 

TOPP = STO( l)^.25*ST0( II " ' ‘ SPT 

C 

C-PRINT out MATRICES ETC., IF DESIRED SPT 

C - SPT 

IF{JPRINT-3> 151,152,152 ’ " SPT 

152 WRITE (KPRINT,51 _ SPT 

' 5 FORMAT (MM A T R I X 0 F “C 0 R R E L A T I ON S*,/t* VARIABSPT 

ILE I IS THE DEPENDENT VARIABLE, ALL OTHERS ARE COSINE AND SINE HAVSPT 

2E S OF V A R I 3 U S PERI fit) S • , / , T 8 8 , • • , / 0 • I " S Pf 

CALL MATPRT (R,M,LVAR » SPT 

C" " SPT 

C-OeCIDE ON PROPER DATA ANALYSIS PATH. _ SPT 

C ’ ’ SPT 

151 JANALY= JANALY4-1 SPT 

"GO TO" ( 200 , 400 , 600T, J AN AL V ' SPT 

C SPT 

C " ■ SPT 

«*««*«[* $PT 

C SPT 

C ' " ' '' ' SPT 

C'STATEMENTS 200-400 INVOLVE COMPUTATION OF OPTIMIZED K-BANDS SPECTRA SPT 

C SPT 

C'FIRST COMPUTE UNIVARIATE SPECTRUM USING SUBROUTINE USPECT. THEN FIND SPT 

C-ALL significant SPECTRUM PEAKS IN UNIVARIATE SPECTRUM. SPT 

C SPT 

“ 2 00 C ALL U SPEC T ( R , T , NOBS', M,B, NBA NOS, STD,! I NT , JPLOT , PHASE, BAN , BA NOS SPT 
I , AMNAME f TNA ME , PHNAME , XNAME, JDOMAN,BANAME I SPT 

CALL PPIKR (B, T, ISAVE,NPEAK,NBAN0S,SIGLEV I SPt 

IF <NPFAK) 201.201,202 SPT 

201 WRITE (KPKINT,66) SPT 

66 FORMAT! *0** * * ♦NO SIGNIFICANT PEAKS FOUND ♦ ♦ ♦ SPT 

“go to 2000 “ SPT 

C SPT 

C-NOW LOOP THRU AN OPTIMIZING PROCESS TRYING TO FIND LARGEST VALUE OF SPT 
C-EITHER T OR R-SQUARE. THE BEGINNING OF THIS LOOP IS STATEMENT 210 SPT 
C SPT 

C-S6T A VALUE nt^ CRIT = 0. THIS WILL LATER BE USED TO STORE THE JUST SPT 

C-PRECEUING VALUE OF T HE OP T I Ml 2 A T ION CRITERION. ALSO INITIALIZE THE SPT 


UlO 
1120 
1130 
1 140 
U50 
1160 
1170 
I 180 
1 190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
_129^ 
1300 
1310 
1320 
1330 
1340 
I 350 
1 360 
1370 
138 0 
1390 
1400 
1410 
142 0 
1430 
1440 
1450 
146 0 
1470 
14B0 
149 0 
1 50 0 
1510 
1520 
1530 
1 54 0 
155 0 
1560 
1570 
15?0 
1590 
160 0 
161 0 
162 0 

163 0 

164 0 

165 0 

166 0 


Xl 


nonni oioono niooo io r«r» I 'O o o 


C-VARUBLE MM TS THE OIRECTtON OF PEAK SLIDING AS WELL AS THE SPT 


C-CnUNTEft F3R THE NUMBER OF STEPS. 

C 

202 LPEAK=l 

CRIT^O.O _ _ 

MM^O 

^10 CALL KSPECT ( R . T , NOBS > M«NPEAK» I S AV £♦ RX » S ♦ JPR INT« XS Aft ,STO, D , RSQ » 

ip-R A se • B A .N 0 S , 8 A N > ' 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


“TEST FOR PRINTOUT 

IFIJPRINT-*3 f ?08t20<?, 20« 
209 WRITE (KPR!NT,6) RSO 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


6 FORMAr*! •OMUL TIPLE SQUARED CORRELATION* • ,F3. 41 
WRITE (KPRINT.7) 

7 FORMAT! • PER I OD /FREQ.* fT20* » AMPLITUDE •,“T40t »T-VALUE* I 

DO 214 L0= UNPEAK __ 

KK=lSAVe!LOI ’ 

IF (JDOMANI 212. 212.213 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


212 BB=BAVOS(KK» 

GO TO 214 

213 8B=BAN(KKJ 

214 WRITE IKPRINT,8I BB. BIKK I «T (KK ) 
8 FORMAT ( 2F 12.4.F 18.4» 


-SELECT CRITERION FOR GOODNESS OF FIT - CHECK FOR IMPROVEMENT 

208 IF JCRITR (RSO.T, JCRITR,LPEAK)-C«IT) 220,220.230 ~ 

“A BRANCH TO 220 IMPLIES NON-IMPROVEMENT. EITHER RESORE OLD VALUE AND 
-GO TO NEXT PEAK (221 OR 223). OR TRY MOVING UP (2221. 


220 IF (MM^l) 221.222.223 


-A BRANCH TO 221 IMPLIES THAT WE HAVE BEEN MOVING DOWN FOR 2 OR MORE 
-STEPS. THIS IS the first NON- I MPROVEMENT , RESTORE JUST FORMER VALUE 
-AND GO TO NEXT PEAK, 


221 ISAVGI LPEAK) =ISAVE(IPEAK) 

MM = -l 

LPEAK=LPEA< fl 

I SAVE! LP6AK) = ISAVE{LPEAK»-1 
GO TO 300 


-A BRANCH TP 222 IMPLIES FIRST HOVE DOWN AND NON- I MPROVEMENT . TRY 
-MOVING SAME peak UP 


222 mm=1 

ISAVE( LPEA< ISAVFdPFAK) fr2 
IF I I S A VE < 1. 1’» AK ) -NRANOS ) 300,225.'225 
2 25 WR ITh ( K PR I NT, 3) 


3 FURMAT (•o**^* ^WARNING*- JUST TRIED TO MOVE A PEAK I 

1 TO Tr«E LlTNGFSI PERIOD (LOWEST FREQUENCY) BAND. TRY WIDER SPFCTRU 

2 LIMITS.') 

GO TO 2000 


SPT 

SPT 

SPT 

SPT 

SPT 

„SPJ 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

_SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 
SPT 
SPT 
SPT 
SPT 
SPT 
' SPT 
SPT 
NSPT 
MSPT 
SPT 
SPT 
SPT 


16T0 
~1680 
1690 
1700 
1710 
1720 
_173 0 
1740 
1750 
1760 
1770 
178 0 
J790 

Tabo 
1810 
1820 
1830 
1340 
J850 
1860 
1370 
1880 
1890 
1900 
|910. 
1920 
1930 
1940 
1950 
I960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21 10 
2 120 
2130 
2140 
2150 
2160 
2170 
213 0 
2190 
2200 
2210 
2220 
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‘A BRANCH TO 223 IMPLIES THAT WE HAVE SEEN MOVING UP AND IT OIDNT HELP«SPT 
-RESTORE'JUST prior' VALUE AND “GO TO NEXT PEAK. SPT 

S PT 

"223 ISAVE(LPEAK)=ISAVEaPEAkr-l“ SPT 

MM--1 . SPT. 

LP£AK*lPEAK4>l SPt 

l$ AVE UPEAKIsISAVEgPEAKI»l SPT 

GO to 300 SPT 

SPT 

^BRANCH TO 230 IMPLIES IMPR OVEMENT. CONTINU E IN S AME DIRECTION. . SPT 

'■ . ‘ ~ SPT 

230 IF (MMI 23l»23It232 ■ SPT 

SPT 

■A 6 R ANCH T 0 23 1 IM PL E R E _ MO V I NG_ DOWN * j-_ CON T I NUE S AME_ D I RE_CJ I ONS PT 

'■ 

231 MM»MM-1 SPT 

CRIT=CRITR(RSQ,tVJCRiTR,LPEAkT Spf 

ISAVE(LPEAK) =l SAVE( LPEAKJ-l __ _ SPT_ 

IF ( ISAV£(LPEAki-n 234,23^*300 SPT 

234 W RITE (KPRINT.4> _ SPT 

4 FORMAT ( *0** * ♦ ^ ~ N G * - JUST TRIEO ‘TO MOVE A PEAK INTSPT 

10 T HE SHORTEST PER 100 j, H I GHEST FREQUENCY) JANO. TRY WIDER SPECTRUSPT 
2M LIMITS' ) “ SPT 

GO TO 2000 SPT 

SPT 

•A BRA NCH TO 232 IMPLIES WE WERE MOVING UP . - C ONTINUE SAME OIRECTION. SPT 

SPT 

232 MM^MM»l _ _SPT 

CRIT*CRitRVRSO,T,JCRlTR,LPE SPT 

ISAV EI L P Ej^K » = I S A V e ( A KJMJI SPT 

IF n SAVE ( lpeak) -Wands j 300,225,225 spt 

SPT 

SPT 

statement 300 CHECKS TO SEE THAT WE ARE NOT OFF THE END OF THE ISAVE SPT 
•VECTOR. IF NOT, DO ANOTHER K-BANOS SPECTRUM WITH NEW ISAVE VECTOR SPT 

•VAL UES . _jF WE _ARJ OFF THE ENO, COMPUTE A NE W FULL SPECTRUM USING A SPT^ 

•K-VAftrA‘fE MOOEL. ALSO CWck“ F^^^ OF ("SAVE EQUAL TO EACH OtHER. SPT 

•IF WE HAVE JUST MOVED INTO A ANOTHER PEAK, REDUCE NPEAK BY ONE AND SPT 
>ACK THE VECTOR DOWN AND START PROCESS OVER AT LP£AK«1. SPT 

SPT 

300 IF (LPEAK-NPEAK ) 320,320,310 "" “ SPT 

320 IF {LPEAK-U 322,322,321 SPT 

321 IF n SAvEILPEAK )-ISAVE(LP'EAK-U) 322,323,322 SPT 

322 IF tLPEAK-NPEAK J 324,210,324 SPT 

324 IF (I SAVEILPEAKI-ISAVEILPEAK^U) 210,325,210 “ ‘ SPT 

323 00 330 LL=LPEAK,NPEAK ^ _ SPT 

330 ISAVE(LL-1 > = ISAVE(LL» ‘ ‘ SPT 

NPEAK=NPEAK-l SPT 

LPEAK=i SPT 

MM=0 _ _ SPT 

■ CRIT=0.0 ■' " ■ SPT 

GO TO 210 _ SPT 

325 NPEAK=NPEAK- 1 *“ " ‘ ' SPT 

DO 331 LL=LPE AK ,NPEAK SPT 

TTTTS A VEl L'O = I S A vl?I L L +i ) spt' 






non on ojo 


lpeak= I 

"MH = 0 
CR IT = 0 .fl 
GO Tfl 210 


C 

C- 
c- 
c- 

C-NEO OP. 

c 


K^BANOS MODEL IS 


A BRANCH TO ^\0 OCCURS WHEN OPTIMIZATION OF THE 

complete. now compute a full spectrum using the 

models SEARCH NEW K-8 SPECTRUM FOR PEAKS AND SEE IF ANY NEW ONES TUR 

IF SOt GO THRU WHOLE OPTIMIZATION ROUTINE AGAIN. 


USING THE FOLLOWING BA NOS IN MODEL* I 


36^ 


BANDS! KK I 


363 


310 IF (JPRINT- n360. 361*361 

361 WRITF <KPRINT,9> 

9 FORMAT! • ICOMPUTING SP ECTRUM 

TF^~( JDOMAN ) 36 2 ,362* 363 

362 DO 36^ LO-l.NPEAK 
KK=ISAVeUOI 
WRITE (KPRINT*8> 

GO TO 360 

DO 3 66 L 0 = I , N PEA K — 

KK=ISAVE(Lni 
366 WRITE (KPRINT*8I BAN(KK) 

360 CALL MULTI ( I S AV E , NPE AK, NB ANOS * R* T * NOBS* M* RX, B * X BAR , STD# 0 t R SO* 

1 PHASE , BANDS, BAN) _ 

KQUNT^KOUNT »• I 

_IF IKOUNT-KTOP) 384 *_3 8 4 3 8 5_ 

'385 wr ite' ( KPrt I NT , 3861 

386 FORMATIMTHE PROGRAM IS PROBABLY CAUGHT 
lINPUT PARAMETERS AND SUBMIT AGAIN*) 

GO TO 2000 

384 CALL PPrKRIB.T, ISA\/E,NPeEK.N0ANOS,S IGLEV) 

IF CNPEAK-NPE EK ) 31 1*35 0* 311 

NPEAK = NPEFK 

LPEAK=1 : 

MM = 0 

CR1T = 0.0 

GO TO 210 


IN A LOOP. TRY DIFFERENT 


311 


'A STABLE"''5TATE' r f HE ' OPT I H I Z ING ROUTINE RESUTS IN A BRANCH TQ 350. 
>nUTPUT RE SUL TS.- . 


LE band SPECTRUM*) 
,B,T, PHASE. NBANOS, A M N A M J * J N A ^ 

oi riT , inn MA Ml 


350 WRI TE (KPRINT,67) 

67 FORMAT! » 0 M U L, T I P 

CALL OUTPUT ( PAN OS. BAN..., 

IBANAME , XNAME , JPrt INT, JPLOT » JOOMAN) 

WRITE (KPKINr.6) RSO 

-CHECK FINAL SPECTRUM FOR EXCESSIVE TOTAL POWER IN SIGNIFICANT BANOS. 

CALL TOTAM! X , NOBS ,NPE AK , I S AV E * M . PHA^E , B * SUMB ) 

IF ( SUMB-TfU’P ) 381,382,382 

302 WRITE ( KP'^ I NT , 3 83 ) 

383 F DKMAT! * IP- JGKAM FAILED BY FINDING TOO MUCH POWER IN THE SPECTRUM 
I*./.* SU-GfST alter ING PARAMETERS OF THE ANALYSIS, SUCH AS**/. 

2* USE WIDER SPECTRUM LIMITS, FILTER DATA, ALTER DELTAP, ETC.*) 

381 GO TO 2000 


"SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

^PT 

SPT 

SPT 

SPT 

SPT 

SPT 

JSPT_ 

SPT 

SPT 

SPT 

SPT 

SPT 

_SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

spt' 

SPT 

SPT 

SPT 

SPT 

„spt 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

.SPT 

SPT' 

SPT 

SPT 

SPT 


2790 
2800 
2810 
2820 
2830 
2840 
2050 
2860 
2370 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 _ 
”298 0 
2990 
3000 
3010 
3020 
3030 
'3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
M^O 
3150 
3160 
3170 
318 0 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3200 
3300 
3 ■*. I 0 
3 320 
3330 
3340 


xllii 
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-STATEMENTS ^OO-SOO INVOLVE COMPUTATION OF THE SIMPLE# SINGLE BAND 
SPECTRUM ONLY. 


4 00 CALL USPECT(RtT »NO0S , M# 8 * NBANOS # STO# JPR INT# JPLOT » PHASE# BAN# BANOS# 
1 AMNAME.TNAME rPHNAME # XNAME* JDQMAN#BANAME) 

GO TO 2000 


**♦******♦♦*♦♦**♦♦*♦*♦***♦*»*♦♦*♦♦*♦**♦*♦♦♦* ♦♦*♦♦♦♦*♦♦ ♦♦♦**♦♦♦*♦*♦♦■•♦♦ 


-STATEMENTS 600-800 INVOLVE COMPUTATION OF A SPECIFIC K-BANO MODEL 
-without OPT! M I Z A T IO N. - 

-FIRST READ IN MODEL PARAMETERS# BUILD VECTOR OF BANDS IN MODEL# THEN 
-COMPUTE SPECTRUM. 

600 READ (KREAOR.IJ NPEAK 

READ (KREADR#?) ( PEAK S I L PE AK ) , LPEAK^ I , NPE AK ) 

~60 6 1 0 \ P E A K = I , N P E AK ' ' 

00 620 LBAN!>= I .NBAND 

IF (PE AKS( LPFAKJ -BANDS! LBANOn 620,630*620 

630 ISAVE( LPEA< > =LBAND _ 

GO TO 610 

620 continue : 

610 CONTINUE 

CALL MULTI ( IS AVE,NPEAK,NBANDS#R#T,NOBS#M#RX#B,XBARf STO»0#RSQ* 

1 PHASE , BANOS. BAN! 

WRITE (KPPIMT.5) 

CALL OUTPUT ( BANOS. BAN#B»T # PHAS E , NB ANOS # ANAME# TNAM E# PHNAME # 

I BANAME.XNAME , JPRINT, JPLOT^JDDMANI 

'WRITE (KPKINT.6) RSO ' " 

2000 STOP 
END 


SPT 3350 
SPT 3350 
SPT 3370 
SPT 3380 
SPT 3390 
SPT 3400 
_S:PT_3410 
SPT 3420 
SPT 3430 
SPT 3440 
SPT 3450 
SPT 3460 
SPT_3470 
SPT 3480 
SPT 3490 
SPT 3500 
SPT 3510 
SPT 3520 
SPT 3530 
SPT 3540 
SPT 3550 
SPT 3560 
SPT 3570 
SPT 3580 
J^PT_3590 
SPT 3600 
SPT 3610 
SPT 3620 
SPT 3630 
SPT 3640 
SPT 3650 
SPT 3660 
SPT 3670 
SPT 3680 
SPT 3690 
SPT 3700 
SPT 3710 
SPT 3720 
SPT 3730 
SPT 3740 




xllv 


n o 


subroutine revix.nj 

DIMENSION XI50I 


"C ♦ * ♦ iT ^ 4c « « ^ * *ftV4'«V4ricV«*4>**«*«4t'***4i*4t4i«4i««4c4t««*«4[4ii^4i*«4t4t4E«4t4t« 

C * «i 

C ♦ SUBROUTINE REV ' ' ' • 

C* 

C * REV TAXES A VECTOR AND REVERSES ITS ORDER. ♦ 

C * ♦ 

_ ^ «* ♦ * * 4tV*ir*4c 4t * 4i«V*«* ** «'^*4( 4>4t4E« ♦ ♦ * * 4t^«* 

c 

j«N*i ; 

MI0=N/2 

DO 10 l=l,MlO 

J*J-l 

SAVE=X(Lr 

X<LI=XIJ) ______ 

10 X(J) = SAVE ■' ■ “ 

RETURN 

END 




xlv 



SVJBROUT INf PREGENCX* BANOS, NOBS* NBANOS ) 
01 HENS TON XUtfBANOSd) 


« 

SUBROUTINE P R E C E N 

PREGEN GENERATES THE SINE AND COSINE PREDICTOR WAVES AND STORES ♦ 
_T HE_M J^N MAT P T X_ X • MATR IX X I S THE D ATA MA TRIX HAVING NOB S ROWS * 
AND M COLUMNS,' HERE' Nb8S= THE NUMBIER OF OBSERVATIONS AND ♦ 

BANOS IS A VECTOR OF PERIOD VALUES IN THE SPECTRUM* * 

♦ 

MATRIX X IS STORED IN VECTOR MODE. SEE IBM SYSTEM/360 SUB- * 

ROUTINE PAC<AGE ( 360A-CM- 0 3X I VERSION III* PROGRAMMER'S MANUAL* ♦ 
PUBLICATION NUMBER H20-0205-3, PAGES 3-6. ESSENTIALLY EACH COL- * 


UMN OF MATRIX X IS STRUNG END-TO-ENO INTO ONE LONG VECTOR. 


THE FIRST NOBS VALUES OF VECTOR (HATRIXI X WILL BE THE DEPENDENT ♦ 
VARIABLE I SEE SUBROUTINE OEPVARJ. THE NEXT NOBS POINTS ICOLUMN * 
TWO) WILL CONTAIN COS WAVE* BAND I. THE NEXT NOBS POINTS (COLUMN ♦ 
THREE CONTAINS SIN WAVE, BAND I, ETC. 1 


CONVERT TO RAO I AN S/OBSER V AT 1 ON AND LOOP THRU NBANOS TIMES. 


20 


10 


L = N03S ^ 

60 lb L BAND ^ ri N BANDS 

AFRE0 = 6. 2d3l8/8ANOS( LBANO) 

EF=-AFREO 

00 20 LneS=l,NOBS 

L = L<^1 

EF=EF»Af REQ 
xi O = C0S< EF) 

EF=-AFREQ 

00 10 LOBS=l .NOBS 

L = L4-l 

eF = FF*-AFREO 
X( LdSINIEF) 


X CONTAINS THE VECTORS 3F SIN/COS WAVES STRUNG END-TO-ENO. THIS IS 
EQUIVALENT TO A MATRIX WHERE EACH SIN OR COS WAVE FORMS A COLUMN OF 
NOBS LENG’^H. 

RETURN 

END 


PREO 

PREO 

PRED 

PREO 

PREO 

'pREO 

PREO 

PREO 

PRED 

PREO 

PRED 

PREO 

PREO 

PRED 

PRED 

^REO 

PRED 

PREO 

PREO 

PREO 

PREO 

PREO 

PRED 

PREO 

PRED 

PREO 

PREO 

_PREO 

PREO 

PREO 

PRED 

PRED 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 

PRED 


10 
20 
30 
40 
_ ^0 
60 
TO 
80 
90 
100 
no 
120 
130 
140 
ISO 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
29 0 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
47 0 
480 
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SUBRnuriNE nEPVAMX,NOBS. JDATA, JOOMAN, JPRINTI 

oi'mens ion X u ) 

rOMMOM KHFanR 

♦ • 

* SUBROUirNE DEPVAR _ 

* "■ ■' 

*_ DEPVAR EITHER READSJN OR GENFR^ES (VtA SIMULA), ADEPENOENT _♦ 

♦ ” V A R [ A B L E VECTOR OF" L E N GT'H' N 0 B S . Y H E " D E P E N 0 E N T V AR I A Bt E ' I S STORE D * 

♦ INTO THE FIRST NOBS CELLS OF X, A MATRIX OF VARIABLES STORED AS A* 

♦ VECTOR OF COLUMNS. SEE SUBROUTINE PREGEN FOR REST OF X. • 

« A««« ««««««*«« AAA 

l‘F“ J DA Tin 5 0 “'“rTadTrOM CAROS, 1 1 OX * 7 F lO^O") 

IF JDATA IS I - CALL SIMULA WHICH GENERATES DATA USING A MONTE CARLO 
SYSTEM ALONG WITH DETERMINISTIC DATA, SIMULA READS CARDS. 

IF (JDATA) 10 , 10,50 

10 READ (KREALTR.l) ( X ( L ) # L = I , NOBS ) 

I > U R M a' T" Y ■ 1 0 X , 7 fY 0 0 ) 

GO TO ZO 

50 CALL SIMULAIX, NOBS, JDATA, JDOMAN,JPRINTI 
20 RETURN 

END y - - . --- 


OEPV 

OEPV 

OEPV 

D6PV 

_DEPV 

oepV 

OEPV 
OEPV 
OEPV 
OEPV 
OEPV 
‘OEP V 
OEPV 
OEPV 
OEPV 
DFPV 
J)EPV 
OEPV 
OEPV 
OEPV 
OEPV 
OEPV 
_DEP_V 
OEPV 
OEPV 
OEPV 
OEPV 
OEPV 
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SUBROUT INE S IMUL 4 ( X, N , J , JOOMAN , JPR I NT I 

dimension X in , a( nm ♦bano( I ) 

CGM'TuN X RE A P K » K r R I NT 


♦ * 

* SUBROUTINE SIMULA _ _ ♦ 

♦ * 

* SIMULA IS A PRIMATIve MONTE^CARLO AND SIGNAL GENE RA TOR FOR GRE A T-* 

V“lN&”ARnFICrAL'OATA, GAUSSIAN OAf A ARE READ IN FROM CAROS AND * 

* A number of sine waves of variable frequency and amplitude are * 

* ADDED* ♦ 

* ♦ 


'ZERO VECTOR X AND PRINT HEADING 


DO 100 L=l.N 
100 X(L)=0.0 

IF (JPRINT) 110*110,120 
_120_ WRITE I XPRl NT 
i'2 I F 0 RM AT i • IM 0 N ' T E 


CARLO 


SIMULATED OATA*,/»*0*l 


READ IN MODEL SPECIFICATIONS - GENERATE SINES AND ADD NOISE 
NSIN-NUMBER OF SINE WAVES IN DATA 
BAMD-FREO. OR PERIOD OF THE ISINTH SINE WAVE 
_S0-STAN0AR0 OEFIATION of THE LSINTH SINE W AVE* (AMPLITUDE) 
SON-STANOAPO OS V . " OF GAUSS I AN NO’l SE 

no READ (KREAOR,U NSIN 

I FORMAT (15) , 

DO ID LSI N= I, NS IN 

READ (KP.EA0R,2) BAND(l),SD 


2 FORMAT (2F10*0) 

IF (JDOVAN) 20,20,30 

30 BAN = BANnU ) 

BAND( I ) = 1. O/BANDI n 
GO TO 21 

20 8AN = I.C/?AN[)( n 

21 IF (JPRINT) 25,25,22 

22 write (KPKINT,3) I.SIN, BAND(l), BAN, SO 

3 FORMAT!* HAND', 13,* - PE R 100^ • , F 10* <», • , F REQUENCV= • , F 10 *4, • , AMPL 
XTUDE= • , F 10. A ) 

25 AFRE0=b,2 J3l rt/OAND( I ) 

___ EF=-AFRFO ^ 

DO 2S L= U N 
EF=EF+AFRFO 

26 A(L ) = SIN(rF)*< I.A1421 
DO 10 1=1, N 

10 X(L)=X(U*A(L)*SO 


-READ' IN’ NOISE SD AND' NOISE' CAROS" 


SI ML 
SIHL 
SIML 
SIML 
AIML 

SI ML 
SIML 
SIML 
SIML 
SIML 
^IML 
SIML 
SIML 
SIML 
SIML 
SIHL 
51 ML_ 
SIML 
SIML 
SIML 
SIML 
SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

_SIHL 

SIML 

SIML 

SIML 

SIML 

SIML 

^SIML 

SIML 

SIML 

SIML 

SIML 

SIHL 

SIHL 

SIML 

SIML 

ISIML 

SIHL 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIHL 

_SIML 

SIML 


10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
no 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240' 
250 
260 
270 
2B0 
290 
300 
310 
320 
330 
34 0 
350 
360 
370 
380 
39 0 
400 
41 0 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
54 0 
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READ C K R e A"D R 7 2T ~ S D H 
IF (SOM I lOOOt 1000,<*0 


SIML 


IF 


« I • 1 ^ I 


I j 0 *4 1 


41 write (KPR IMT,5) SON 

5 FORMAT ( *CP ANDO‘S GAUSSIAN NOISE ANPLITUOE*»tFlO*4) 
50 RFAO (KREA0R.6J (ACL}»L=1»N> 

6 FORMAT (4X, 10F6.4 »■ ““ ~~ 

DO 60 I =l.N 

60 X(U=A(LI*SON<^X(LI 
1000 RETURN 

END ■ "■ ■■ 


SIML 

SIML 

SIML 

SIML 

SIML 

_S^!ML_ 

SIML 

SIML 

SIML 

SIML 

SIML 


550 

560 

570 

580 

590 

600 

610 

620 

630 

64 0 

65 0 

66 0 
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SUBRQUTlNe MATPRT(R,H*LVARI 

DIMENSION R n } ,X( 151 # LVAft( II 

KRF AOR , KPRJ NT 

♦ ***♦**»*****♦•**!»♦**♦*♦**♦♦♦♦»♦♦♦♦♦*♦*♦♦♦★♦♦*♦♦♦*♦♦♦•♦♦♦♦♦*♦♦♦***♦** 
♦ ' ♦ 

♦ SUBROUTINE H A T P R T • 


J*_ MATPRT PRTNTSOUT JHE LOWER TR I A NGUL ARJM AT R I X X WHER E X IS STDR- * 
*’ to’ IN'^^One iV "see reference list ED IN SUBROilTINE PREGEN. ♦ 

m ♦ 

4i ***«**«[i4i*i^««* ************************** 

-LOOP THRU as many TIMES A$ NEEDED TO PRINT WHOLE MATRIX, EACH TIME 
-PRINTING 13 C O L'UMN S BY NRO ROWS« ' 

Nl»l _ „ 

N2=15 

110 IF IN2-MI 100, 100,10 ^ 

10 N2=M 

10 0 WRITE IKPRlNT,n ( LV AR I L V ) , LV»N1 , N2 ) ________ 

1 FORMAT! * 0 ' , 2X ,'l 517) 

WRITE I KPR [NT, 2 1 _ _ _ _ 

2 FORMAT! »0' ) 

-PRINT ONE ROW AT A TIME, MATRIX R IS MODE I, UPPER TRIANGULAR", 

-T HERE FORE REVERSE S UBSCRIPTS TO MAKE LOWER TRIANGULAR, 

00 20 LR0W = N1,M _ 

IF ILR0W-N2I 21,21,22 

21 LIMIT^LROW 

GO TO 2 3 ' 

22 L IMIT = N2 . 

_ '23 L = o 

00 30 LCOL=Nl, LIMIT 
L = L + l 

CALL LOC! LCOL,LROW,K,M,M, U 
30 X! L >=R (K ) 

20 WRITE !KPRINT,3) (LRQW , ! X ! K I , K= I, L ) I 

' 3 " FORMAT! IX, I ^ ,2X,15F7. 31 

Nl=Nl*lS 

N2=Nl+14 

IF IM-NU 100C,U0,110 
1000 RETURN 
END 




1 
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SUSR3UTINE USPECT(R,T,N , M. B. NB AMDS t S TO, J PR tNT , JPLOTt PHASE* BANt 

I BANOS, AMMA^t , TNAMF, PHNAME , XNAME » JOOHAN , 8ANAM6 ) 

0 I .IE NS I ON R { 1 S , T { 1 I , B a J , STD { I I , P H A 5 E i I >* BA N{ U, BAN DS 111 MN A ME I I ) 

ItTNAME ( IKPHNAMEJ I » , XNAME ( I r,"BANAMF{ lT 
COMMON KftEAOR.KPR INT 

««[«««« «<[««*« 1^4^ « If 4c 4. 4c«» *«* «*4>4i 

4i ♦ 

SUBROUTINE U S P E C T 

♦ ” . * 

♦ USPECT COMPUTES THE UNIVARIATE SPECTRUM BY MULTIPLYING THE CORR- * 

♦ ELATIO.N OF THE TIME SERIES WITH A SIN OR COS PREDICTOR BY THE * 

♦ RATIO OF STANOARO DEVIATIONS OF THE TWO, THIS IS DONE FOR EACH * 

♦ PREOICroR. LET R=THE CORRELATION OF F<T» WITH A SIN OR COS WAVE ♦ 

♦ OF SOME ARBITRARY WAVE LENGTH. LET S =THE STANDARD DEVIATION _* 

l^QF'THE TIME SERIES” AND .707= THE STANOARO OEVTATION'oE ♦ 

♦ ICTOR WAVE. THEN AIFI = R+IS/.707) * 

♦ WHERE A(FI=TH6 AMPLITUDE OF THE COMPLEX SPECTRUM AT F. THE SIN ♦ 

♦ AND COS COMPONENTS ARE COMBINED AT EACH HAVE LENGTH AND PHASE • 

♦ ANGLE IS COMPUTED. T-'VALUES ARE ALSO COMPUTED FOR EACH BAND. ♦ 
__R = UPPER TRIANGULAR CORR ELAT MATRI X, IN STORAGE MODE 1. SEE * 

♦ rSM PUBLICATION REFERENCED IN SUBROutiNE PREGEN. (INPUT) ♦ 

♦ T=VECTOR OF T-VALUE5. (OUTPUT) ♦ 

♦ NOBS=NUMBER OF OBSERVATIONS. (INPUT) ♦ 

♦ M= TOTAL NUMBER OF DATA VECTORS. M=NBAN0S*2 «■ 1. ONE VECTOR FOR ♦ 

♦ the DEPENDENT VARIABLE AND TWO FOR THE SiN/COS WAVES AT EACH WAVE* 

♦ LENGTH. (INPUTI. 6=VECT0R OF SPECTRUM ESTIMATES. (O UTPUT) * 

♦ N8ANl)S = LENGTH UF VECTORS T, 8 AND phase, riNPUtr’ ♦ 

♦ STO=STANDARD DEVIATIONS OF ALL VARIABLES, (INPUT) • 

♦ JPRINT^PRINT CONTROL - SEE MAINLINE. (INPUT) • 

♦ JPL0T=PL0T control - SEE MAINLINE. (INPUT) ♦ 

♦ PHASE=VECTOR OF PHASE ANGLES. (OUTPUT) * 

♦ ^ ; 

"if 4( 4t f f 4i4>*'f * * * * * « ♦ * ♦ if V« * 4t ♦ If **if 4c**4i**4c*4i*)it4c4i «4iif4E#4(#«4i4i4ciff«4c«4t 

lband=o 

EXTRACT VECTOR OF CORRELATIONS BETWEEN DEPENDENT VARIABLE AND EACH 
■PREDICTOR VARIABLE. 


00 10 L=e,M 

CALL LOCU. L.KOOL ,M.M,l) 
10 B(L)=R(KOOL) 


-LOOP THRU THE SPECTRUM COMPUTATION N8AND TIMES. 

F AC T = SO R T ( F 10 A 1 1 N -2 ) Y”“ ”” 

DO 20 L=2,M,2 
K = L + l 

LflAND = L BAND*- I 


C-CO=CnS COMPONENT OF SPECTRUM ESTIMATE 
C-SI = SIN component i>E SPECTRUM ESTIMATE 


PMCHNAL' PAGE Ig 
m POOR ^UAUKf 



C-TC»T-VALl>H FOR COS COMPOM€NT_ ____ 

C^fTif-VALUE FOR SIN COMPONENf ' 

C „ 

CO*B(L I*ST0( II/.707 

SI’BCU ♦STO( n /- 707 . 

TC=8(L ><^F/\CT/SORT( UO-BI L»**2I 

TS=fi(K)*F^CT/SORTI 

C 

C-COHPUTE SPECTRUM AMPLITUDE, T-VALUE AND PHASE ANGLE 

C 

BILBAN0) = SQRT (C0«*2 SI*’«'2) 

■' t(L8AND J = SQRT (TC**2 TS**2» 

20 PHASE! LBANDI =ATAN(C0 7SI ) ^ 

C 

C-PRINT AND/OR PLOT IF DESIRED ’ 

^ c 

__JPR=JPRlNT-l ■ 

JPL»JPLOT-l 

IF (JPR> 50,50,60 

50~I F ( J PL ) I CO 0,1000,60 

60 WRI TE ( KPR I NT, 1) . .. . _ 

I FORMAT! »0SINGLE BAND SPECTRUM*) 

CALL OUTPUT! BANOS , BAN, B,T, PHASE, NBANOS, AMNAME, TNAME .PHNAME 
■ 1,8ANAME,XNAME, JPRINT, JPtOT, JDOMAN) 

1000 RETURN 

END 


JUSPT 550 
USPt 560 
USPT 570 
USPT 580 
USPT 590 
USPT 600 
_U SPJT 610 
USPT 620' 
USPT 630 
USPT 6<»0 
USPT 650 
USPT 660 
US PT 6 70 
USPT 680 
USPT 690 
USPT 700 
USPT 7^0 
USPT 750 
U SPT 760 
USPT 770 
USPT 780 
USPT 790 
USPT 800 
USPT 810 
US PT 8 20 
USPT 830 


I^GINAi; PAGE J3 
POOB .QUALmr 


lii 
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SUBROUTINE PPI KR ( 6 , T , I SAVE * NPEAK, NB ANDS t S IGL EV ) 

DIMENSION BfU,Tri)*ISAVE(l» 

'«r***ir*n*i('* ♦♦ ♦♦ « « if ** **«**««# *V*~* 

♦ _ _ • 

♦ SUBROUTINE P P I K R ♦ 

♦ « 
* PPIKR RETURNS THE VECTOR INDEX OF ALL SIGNIFICANT SPECTRUH PEAKS.* 

♦ BMNP'JT_ VECTOR OF AMIPLIJUDE SPECTR UM ESTIMATES - ♦_ 

-INPUT VeCfOR ~0F T-vaLUES FO EACH SP EC TRUM ESTIMATE * 

* ISAVE-OUTPUT VECTOR OF INDEX NUMBERS OF SIGNIFICANT PEAKS _♦ 

♦ ' NPEAK-OUTPUT SCALAR - NUMBER OF PEAKS FOUND V 

• NSANOS- INPUT SCALAR - NUMBER OF BANDS IN SPECTRUM 

* SIGLEV-INPUT SCALAR - CRITERION WHICH T-VALUE FOR A PEAK MUST * 

_* EXCEED IN ORDER TO BE RETAINED » 

* * 

*^*#****^4141 ****** ****** 41 * ************* ****** 4 t 


4*^^ y. P E 4 KS . RET A I N • S J GN I F I C ANn ON ES * BU I L 0„ I S AVE, 


NPEA K^O _ 

00 10 l=2«n’ban6s 

IF (B{L-U-8(D) 
IF )-8(U ) 

IF IT(L I-SIGLEVI 
NP£AK=NPEAK*l 
ISAVE«NPEAK»=L 
ZONTnUE 
RETURN 


20*10. 10 
40 . 10, 10 
10 * 50,50 


no 


FUNCTION CRITRf RSOfT*JCftlTR,L» CRIT 

OIHENSION Tin ' ' CRIT 

_C C RIT 

C ************** i^*****tim'*****^************************it*****^******m*** CRIT 

C ♦ * CRIT 

C * FUNCTION C R I T R " “ '♦ CRI^ 

C * ■* CRIT 

C > ■ CRITR SELECTS A CRITER ION FOR' tnF SPEC TRUM' OPT ♦ CRIT 

C » SHIFTING. S EL ECTS EITHER R-SQUARE OR THE APPROPRIATE T VALUE, ♦ C RIT 

C • ' ♦ CRIT 

C **m***************************^************************************** CRIT 

"C - - -- CmT 

C-SELECT APPROPRIATE VALUE FOR CRITERION, CRIT 


RSa- INPUT SCALAR - VALUE OF R-SQUARE 
T-INPUT VECTOR - T-VALUES 


CRIT 

CRIT 

CRIT 


JCRITR-CONTROL digit 


SELECTION OF EITHER RSQ OR TIL) AS CRITERIONCRTT 


0*RSQ AS CRITERION 

laT-VALUE AS CRITERION 
L-INOEX VALUE OF VECTOR T 


IF (JCRITRI 
10 CRITRsRSO ■ 
GO TO 30 
'20 CRITR=T (L) 
30 RETURN 
END 


10 * 10*20 


CRIT 
CRIT 
CRIT 
CRIT 
CRIT 
CRI T 
CRIT 
CR I T 
CRIT 
CRIT 



liv 

itStlGINAi; PAGE IS 
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KSPT 

" ■ ■■ ■" ~ ' KSPT 

SUBROUTINE < SPECT (R , T ,Nt M,NPEAK, I SAV « RX » B t JPR I NT ,XBAR * STD * D* RSO f KSPT 
IPH4SE* BANOS, BAN) KSPT 

DIMENSION I save MOO) ,R (1) . T( I WRX( 1 ) , 6 (1) , X8AR ( I ) ,STD( U.OUI, KSPT 
iPHASe i I ) , BANDS ( L ) , RANTI ) "■ , ANSl'i bl *TSAVf li ,RV| lbb> KSPT 

COMMON KRE AOP ,KPRINT KSPT 

-CONVERT ISAV TO I SAVE ♦ FROM FREQ BAND VECTOR TO VECTOR OF SIN/COS PAIRKSPT 

KSPT 

8521 FORMAT (•0JKL**,I5) KSPT 

K=0 KS PT 

bo lob t*l,NPEAK KSPT 

K » K ♦ I K SP T 

~l SAVeiK )« ISAV IL )*2 KSPT 

K = K*l KSPT 

lOO'iSAVEIK »=ISAVE(K-lUl‘ KSPT 

K L 

-SELECT SUa-MATRIX OF PREDICTORS ETC FROM R, ACCORDING TO ISAVE. KSPT 

_ KSPT 

CALL OROER(M*R,l *K ' ,1 S AVE * RX ,RY ) ' ~ KSPT 

KSPT 


-invert K-ORbER MATRIX OF PREDICTOR INTE RCORRELATI ONS t RX 

CXCl MINV(RX,K /DET^Tb ,TI ~ 

IF IJPRINT-T) 300*400,^00 

'400 WRITE (KPftINT*?) DET 

7 FORHATI •O*,/, *0 DETERMINANT » V*E>0.15) . 

-COMPUTE RE GR E SS T CNS ( SPEC TRUM ) FOR NPEAK SIZED MODEL, 

300 CALL MULTRIN,K,xaAR,STO*0,RX,RY, ISAVE, B.SB^^ANS) 

RSQ=ANS( 2»*»2 

-C ON V E R T SPECTRUM TO AMPLITUDES By“ COMB IN I NG S I N /COS . PUY TnT 0 ” 0 ^ 
-ALSO C ON VERT T« 

NFRE0*0 _ __ _ _ ' 

00 200 L= 1,K, 2 ' 

NFREO=NFRFQ^l 

PHASE! NFR6Q » =ATAN( BC U /B( L> II ) 

B ( NFRE 0 ) =SORT ( B ( L > 2 B < L ~ 

200 T( NFReO)=SORT(Ta»*^2 ♦ t(LYl)*Y2l 

10 RETURN 

END 


KSPT 

J<S^T_ 

KSPT 

_KSPT 

KSPT 

KSPT 

KSPT 

KSPT 

KSPT 

KSPT 

KSPT 

KSPT 

KSPT 

KSPf 

KSPT 

KSPT 

KSPT 

KSPT 

J^SPT 

KSPT 

KSPT 

KSPT 


SUBR3UTINE MUL T U I SA VE*NPE AK • N3AN0S,R , T *NtH ,RX t B « X5AR, STD t 0*AS0« 

I PHA<;e, BANOS ♦8Ar4» 

01 HENS ION i SAVEUi tR<U*TU),RX(Uf B<n,X6AR(l ),ST0(n,0( ILiPHA^ 


I ( IK BANOS ( I ) ,3AN( U «KSAVE( 20US< 201 *TSPEC(20» *TPH(20) 
COMMON KRE AOR.KPRINT 


♦ ’ ’ ~ ' ♦ 

• SUBROUTINE MULTI * 


MULT 

10 

MULT 

20 

MULT 

30 

MULT 

40 

MULT 

50 

MULT 

60 

MULT 

70 

MULT 

80 

MULT 

90 


MULT 

MULT 


♦ ♦ MULT 

♦ MULTI COMPUTES THE FULL SPECTRUM OF A TIME SERIES USING A K + l • MULT 

'* 'BANOS MODEL WHERE K=THF NUMBER OF PEAKS IN THE GPTQMIZEO ' MULT 

♦ MODEL. THE EXTRA RANG IS THE WAVELENGTH IN QUESTION. THE OTHERS ♦ MULT 

♦ ARE THE SIGNIFICANT PEAKS. THE VARIANCE ACCOUNTED FOR BY ♦ MULT 

♦ THE SIGNIFICANT PEAKS IS THUS »ACCOUNTEO FOR* IN EACH ESTIMATE. * MULT 


100 

120 

130 

140 

ISO 

160 

170 


♦ WHEN the current BAND IS ONE OF THE SIGNIFICANT PEAKS. A SKIP 

♦ OCCURS. VARIABLES DEFINED IN MAINLINE. 


* 

___ _ ♦ 

♦ ’ “ ' ' ' * 


-LOOP THRU ONCE FOR EACH FREQ BAND. NFMOO IS THE NUMBER OF SIGNICIC- 
-CANT PEAKS IN MODEL PLUS A CURRENT BAND. 

NFMOO=NPEAK+ I 

'do 310 LF^l.NBANOS 



MULT 
MULT 
MULT 
MULT 
MULT 
jiyLH 
MUL T 
MULT 
mult 
MULT 
MULT 
MULT 


CONSTRUCT NEW ISAVE VECTOR 

CALL BU ILD ( K SAVE . I SA VE,LF7NFHdd , I 
SKIP I F THIS BAND I S IN t HE MODEL ALrEAOY 


IF nSKIP) 320, 320.330 

BRANCH TO 330 IMPLIES THAT THE CURRENT BAND IS ONE OF THE SIGNIFICANT 
•PEAKS. DO AN NPEAK SPECTRUM ANO PICK OUT CURRENT BAND. 

330 CALLKSPECTI R, S.N.M .NPEAK, I SAVE , RX, TSPEC . O.XBAR , STD.D.RSQ. TPH, 


IBANOS.BANI 

•TSPEC CONTAINS ALL MODEL ESTIMATES. PICK THE ONE FOR LF. 

DO 331 LLZ=l.NPFAK . 

IE {LF-ISAVf ( LLZn 331,33 3 , 331 

"333 8( LFMTSPEC ( LLZ » 

PHASE I Lf KTPHILLZ) 

TILFI =S(LLZ) 

331 CONTINUE 

GO TO 310 


:-8RANCh to 320 rMPLiES^THAt THIS IS A FULL NFMOO MODEL. 00 SPECTRUM 


MULT 

MULT 

MULT 

MULT 

mult 

mult 

MULf 

MULT 

MULT 

MULT 

MULT 

MULT 

MUL f 

MULT 

MULT 

MULT 

MULT 

mult 

MULT 

MULT 

MULT 

MULT 

MULT 

MULT 

MULT 


180 
190 
200 
210 
220 
230 _ 
240 
250 
260 
270 
260 
290 
300 
310 
320 
330 
340 
350 
360 
370 
38 0 
390 
400 
410 
420 
430 
44 0 
450 
460 
470 
460 
490 
500 
510 
520 
530 
540 


Ivi 




c-and pick out current band. 5ykl_-”y_ 

^ MUt.T 560 

320 CMLLK3PECT{R»S*r'j#M*NFriGD» KS AV E t RX t T S P . . ^ 570 

1 BANOSfBAN) ’ '■ 580 

00 321 LLL = l,NFMOO 590 

IF( LF-KSAVE( LU. ) J 32l.322»32l ^ULT 600 

322 T(LFJ = SfLLL) 

aiLFi =T5 pec(lCCJ ^20 

PHASE (LFI=TPH(LLL) — WULT 630 

‘321 CONTINUE HULT 640 

310 CONTINUE WpLT 650 

RETURN 660 

FNn MULT 670 


ooo j noino 


SUBROUT INE ftUILO<KSAVE» I SAVE *L FtNFMOD# I SKIP»NPEAKI 
DIMENSION KS AVE( I » * I SAVE ( I I * AS AVE ( 50 1 .RANCt 30) 


♦ 


SUBROUTINE BUILD 

♦ BUILD CONSTRUCTS VECTOR OF RANKED FREQ BAND NUMBERS GIVEN A 

♦ PARTICULAR MODEL AND A CURRENT FREQUENCY. 

♦ KSAVE - OUTPUT VECTOR TO BE BU ILT - 


BUD 
BILO 
BILD 
BILD 
_BIL_0_ 
BI LO 
BILO 
BILO 
BILO 
BILO 
BILO 


C'~* ISAVE -^INPUf” VECTOR CONTAINING NFREO SIGNIF PEAKS 
C ♦ LF - INPUT SCALAR - NEW FREQ BAND TO BE ADDED (CURRENT ONEI _ 
C'V “ NFMOO - INPUT - NUMBER OF FREQ BANDS IN NEW MODEL 
C • ISKIP - OUTPUT CONTROL DIGIT TO KEEP FROM COMPUT ING EST I MATE 

OF ONE OF THE OLD PEAKS AGAIN 

C ♦ NPEAK - INPUT - number OF FREQ BANDS IN OLD MODEL 


BILO 

BILD 

BILO 

BILO 

BttO 

BILO 




C 

c 

c 

C-SEE IF LF IS ALREADY IN MODEL 




00 10 L=l»NPEAK 


IFILF - ISAVE(L)) 
10 CONTINUE 


10# 99. 10 


10 
20 
30 
40 
_50. 
60 
70 
60 
90 
100 
11 0 _ 
120 
130 
140 
150 
160 
170 


BILO 

BILD 

BILO 

BILO 

BILO 

BILO 


BILO 

BILO 

BILD 


180 

190 

200 

210 

220 

23Ql 

240 

250 

260 


STORE VECTOR OF IS AVE PLUS Lf INTO ..A 
VECTOR RANG 


00 20 L=l. NPEAK 

20 ASAVE( L >^I SAVE! L» 

ASAVei NFMOO I=LF 
CALL RANK! ASAVe.RANC.NFMOO) 

-B UILD KSAVE OF RANKED FREQ BANOS 

00 30 L = l. NFMOO 

KOQL-RANC(L ) 

30 KSAVE( KOOL) =ASAVE(L) 

ISKIP^O 

_G0 TO 1^0 

99 ISKIPsl 

100 RETURN 

END 


RANK STORED IN BILO 270 
BILO 280 

6 0 

BiLD 300 
BILO 310 
BILD 320 
BILO 330 
" BILO 340 
BILO 350 
BILO 360 
BILO 370 
BILO 380 
BILO 390 
BILD 400 

BI LD 4 10 

BILD 420 

BILD 430 

BILO 440 
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SUBftOUTI NC OUTPUT (a ANnS* BAN, 8* T ♦PHASE. NB AMDS* AMNAHE,TNAME* PHNAME. 
IBANAMEt )(NAMe * JPRINT, JPLOT, JOOMANJ ’ 

D I ME NS I ON bands { I ) » 8 AN( I »_♦ B (I t,T(I>,PHASE(Il. AHNAHE ( 1 > .TNAHEC U t 
IPHNAME ( II ,0ANAME( I ) VxNAMETU 

COMMON KREAOR.KPRINT . .. 

' * * 

* SUBROUTINE OUTPUT - 

♦ OUTPUT PRINTS ANO/OR PLOTS SPECTRA WITH OOCUMENTAT ION. * 

■'* * 

-PRINT OUT DOCUMENTATION A NO VECTOR ORDER IN APPROPRIATE DOMAIN! 

^rNOT^CALE value's a'NO PLOT 

TOP 1 = 6 *283 18 

8TOPI*-TOPI 

" BMAX=e(l) 

00 64 LflAN0 = 2«NBANDS 

IF ('BMAX-B'rtBANOT) 68,64,64 

68 BMAX = B( LBAND ) _ - 

54 CONTINUE 

TMAX=T(LI 

00 66 L8AN0=2.NBAN0S 

IF ITMAX-TI LBAND M 67,66,66 - 

67“ T M A'X= T ( LBAND 1 

66 CONTINUE — — 

IF IJOOMANI 60,60,50 

•REVERSE VECTORS IF FR EO DOHA I N I S' USE 0 ' - ‘PR 1 NT ~ • REVE R SE “AG A I N 


50 CALL ReVrBfNBANOSI 

CALL REVIT.NBANDSI 

"CALL REV(PHASE,N8AN0S) 

IF (JPRINT) 52, 52,51 - 

51 WRITE IKPRINT.2I 

2 FORMAT! • OFREQUENCY * , 5X , * AMPL 1 TUDE * , 5X, * T-VALUE * , 8X, * PHAS E * , / , * 0 » 
DO 55 LB ANO = I ,NBANOS ' 

55 WRITE !KPRINT,3I B AN ( LB AND ) , B (L BAND ) ,T I LBANCH ,PH ASE ! L BAND ) 

3 F0RMAT(F10- 3, 5X , F9. 3 , 5X, F 7 • 3, 5X,F 8. 3 1 

52 IF IJPLOTI 54.54,53 

53 CALL PLQT2I B , 8M A X . 0. 0 , T , TMA X ,0* 0 , N3 ANOS , B AN, AHNAME ,TNAME , BANAME I 

CALL PL 0T2 ( PHAS E_,Tqi^I , BT QP I , T , T MAX, 0* 0,NBANDS, SAN , PHNAME , TNAM E, 

T^BANAME) 

54 CALL RE VI fl.NBANOSJ _ 

CALL REVIT.NBANDS) 

CALL REV! PHASE, NBANDS) 

GO TO 100 

60 IF U^RINTI 62,6^,61 — — 

lbl~HR'lTe <KPRINT,4> 




O' O' 

\jt t'j 




lx 
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SUBROUTINE PLOT2 ( VI t TMAXl *YMINl,Y2t YHAXZt YMIN2tNP,XAX»YNAHlt 

^0 I MENS lOrvi ( I » ,Y 2 ni ,L INFI 11 i ) ,XAX ( 1 K XNAH <5 K, YN AMljgj » YNA^ 

~T^MU 'i " k R E'A b R . K P R N T 

DATA KBLNK/* •/ *KBORO/ *.*/ ,KAL /•«'*/ t KPLOT/V»*/ 


PLT2 

PLT2 

PLT2 

PtT2 

PLT2 


PLT2 

PLT2 

PLT2 

PLT2 

PLT2 

PLT2 


C — 


PLOT2 GENERATES TWO OVER/UNOER GRAPHS PER PAGE AND CALIBRATES 

THEM* YI=VERTICAL VARIABLE I, YMAXlt YH I N I VERT !C AL PLOT LIMITS. 

Y2. YM4X2» YMINZ^SIMUAR FOR SEC ON VARIABLE^ NP»NUH 6ER OF POINTS,. 

Ta 1^Y AXIS“('CnMMON“TO BdfH' VAR^f^^^ 

WRITE (KPRNT.2) YNAM 1 * YNAM2 • XNAM . 

2 FORMAT! / T5.»**/ T5 • • ♦ • . T25 . * * • .T80, • * • / T5 , • * • , T25, • PLT2 

1 T80. **♦♦•/ T 3, *♦♦***• *T 18 * ,9A4*T73» * ♦♦**♦♦*♦♦♦♦♦ *PLT2 

2 ♦9A4^/T4, *!*♦♦* fT2Sf • « T80t*****/T5t***#T25»*** » TSO» •♦•/IX.SA^l PLT 2 


-GENERAfe SCALE FACTORS 

SKALl=S0.5/( YMAXl-YMINU 
SKAL2-S0.5/ (YMAX2-YMIN2I 
I^GENERATE CALIBRATION AND OUTPUT 
CALOl = YMINl 

CAL02 * YMIN2 


IT. 


PLT2 
PLT2 
PLT2 
PLT2 
PLT 2 
PLT 2 




CALI = YMAXl 

CAL2 - YMAX2 . „ 

CAL5l=YMAXl- (*5«( YHAXl-YMINl ) ) 

CAL52=YMAX2- (- 5*( YMAX2-YMIN2) I 

WRITE (KPRNTtll CALO I *CAL 51 »CAL l»CAL02t CAL52#CAL2 
I FORMAT I* *»FlA,^fF_24. ^tF23.A «F 12.4f F23.AtF2A»A) 
^NERAte VERTI 

00 10 L * It 1 1 1 



io 


20 


LINEILl = KBORO 
LINE! I I = RAL 
LINE(26I = KAL 
LINE(5U = KAL 
^L1NE(61 I 
LINE! 86) 

L INF! 1 I 1 )= KAL 
00 20 L=52.60 
LINE(U = KflLNK 
WR I TE( KPRNT, 4 ) LINE 
FORMAT ( lOX, I U A1 ) 


— OeClOE HOW 

NWine = 


MANY LINES 
32 /NP 


- -LO OP FOR 
DO "30 


NUMBER 
LP'= I 



BETWEEN 


OF POINTS^ 
* NP 


EACH 


—CLEAR LINE ANO GENERATE GRID POINTS 
DO ^ = 1 f I II 

AO LlNc»L*~KtiLNK 
^_LIN 6 ! l»=KBnRD _ 

L l'NFr?6 I^KflORO 


10 

20 

30 

40 

60 
70 
BO 
90 
100 
UO 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230_ 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
47 0 
480 
490 
500 
510 
52 0 
530 
540 



L I N6( 5 I ) =K BOR p 

LI NE~(6l »=KBGRC) 

HNe( 86 »=K80RD ' . — — 

LINEU ll l^KBORO 

c 1 - — - - 

■"C--TEST POR OVER OR UNDER RANGE AND TRUNCATE OUT-'OF-BOUNOS VALUES 

IF I Y I ( L P »_r_Y H A X 1) 2Q0 0« 2000* 2100 

2100 YA = YMAXi 

GOT05100 — T 

~20b0 IF ( YU LP >-YH I Nl ) 5000 * 5200* 5200 

5000 YA = YHINI 

GO TO 5100 

5200 YA-YULP) 

5100 IF CY2('lP>-YMAX2> 200# 200* 210 

210 Y8 = VMAX2 

GO TO 510 

200 IF I Y2 I LP I-YM IN2 ) 500# 520^ 520.^ , 

~ 500 YB - YMIN2 

GO TO 510 — 

520 YB=«Y2(LP» 

C_ _ ; - 

C—’NOW' GENERATE POINT INDEXES 

510 Kl * «VA -YMINl »*SKAL1>1 #0 — 

__ K2 i ( YB -YMIN2I<«SKAL2>61,0 

LINFIK1I= KPLOT__ 

LINEU2 J* KPLOT 

C 

‘(C--OUTPUT HORIZONTAL SCALE VALUE AND LINE 

WRITE (KPRNT ,31 XAX<LP)#IINE 

3 FORMAT (IX,F 8«A#lX#lllAn 

C 

C^COWXrNE”AN0 GENERATE BLANK LINES FOR SPACE 

_ __ (NWIDE) 30,30*110 

110 DO 95 L = l*Ul 

95 L INE(L )=KBLNK _ — 

_ line (I 1=KBriRD 

LI NE( 51 » =KBORD . 

X I N E (26 ) = K Bn R 0 

LI NE(6l I =K80RD _ , 

' ' LINE(86)=KBGRD 

LINE( 1 IU=KBORO , 

c 

C--OUTPUT PLAIN LINE 

' 'DO' 90 LOOP=l f NMIOE 

90 WRITE (KPRNT.^I LINE .. 

30 CONTINUE 

C — GENERATE RIGHT BORDER AND OUTPUT 


OJ 00 L = 1. I 1 1 

5b"L I NE ( L 1 = K fl 0 R 0 
L INE( I ) = KAL 
L ! 126 > = K AL 

LI N6( 51 1 = KAL 
L INF(61 I = KAL 
LU.E(B6I KAL 
LI NE( 1111= KAL 


_PtT2J55 0 . 
PLT2 560 
_PLT2 570 
PLT2 580 
PLT2 590 
PLT2 600 
PLT2 610 
PLT2 620 
_PtT2 63 0 
PLT2 6^0 
_^tT2_650 
PLT2 660 
PLT2 670 
PLT2 68 0 
JPLT2 690 
PLT2 700 
_PLT2 .7 10 
PLT2 T20 
PLT2 7 30 
PLT2 740 
PLT2 750 
PUT 2 760 
PLT2 770 
PLT2 780 
PLT2 79 0 
PLT2 800 
_PLT2 810 
PLT2 820 
PLT2 830 
PLT2 840 
PLT2 850 
PLT2 860 
PLT2 870 
PLT2 880 
PLT2 890 
PLT2 900 
PL T2 910 
PLT2 920 
PLT2 930 
>LT2 940 
PLT2 950 
PLT2 960 
PLT2 970 
PLT2 980 
PLT2 990 
PLT21000 
PLT21010 
PLT21020 
_PLT21030 
PL'T21040 
PLT21050 
PLT21060 
PLT21070 
PLT21080 
_PL T21090 
PL T 21 100 
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SUBROUTINE TOT AH( XtN, NPEAKf IS AVEt M , PHASE* 6* AMP I 
DIMENSION xm flSAVEIU ,PHASE(l)*BI II 


c * 

c ♦ SUBROUTINE T 0 

C ♦ 

C * ' TOTAM COMPUTES 

c ♦ SO I OS WHICH H4V 
TATcfNC PROGRAM 
AND COS INE WAVE 
PUtES THE AMPLt 
WARNING 
CONTAINS THE OR 


TAM ' . . • 

THE RESULTANT AMPLITUDE OF THE SUM OF NPEAK SIN- • 
E THEIR amplitudes and PHA SES PETR EHI_NED __B Y JTH1__*. 
.7 — t OTAir S I mMV "SELECTS THE APPRQPRiATE SINE ♦ 
S FROM X* WEIGHTS THEM* SUMS THEM AND THEN COM-_* 

TU06 OF The sum. * 

— THIS WILL DESTROY ROW L Of _X WHICH 

IGINAL DEPENDENT VARIABLE* * 






-LOOP thru THE FOLLOWING NPEAK TIMES TO CREATE THE RESULTANT WAVE INX 


DO' 100 LP=l*NPEAK 
K= ISAVEILP ) 


KC*K*2 

KS^KC^l 


THE FOLLOWING COMPUTES THE PROPER AMPLITUDES FOR THE COSINE 
COMPONENTS* BASED ON PHASE' AND AMPLITUDE ESTIMATES* 


AND SINE 


PC =*1.0 

PC*PC^( tan ( PHASE! K U *■ 
PS«PC*TAN<PHASE(KI j 


PC) 


■SELECT SINE AND COSINO VALUES FROM X A 

SELECT SINE AND COSINE VALUES FROM X* SUM AND PUT INTO X* ROW 1. 


100 


DO 100 L=l*N 
CALL L0C(KC,L*J.M,N*0) 
CALL LOC(KS*L*I tM,N,0) 
X(L)=K(J)»PC ♦ X(I)*PS 


^■^QM'Puf E APLLlfUOE OF RESULTANT 

C 

SX'0.0 

00 10 L*l.N _ 

10 SX = SX4-X!L J 

SX=SX/FLOAT(N) 

SXS=0*0 
00 20 L=l*N 

20 SXS=SXS*< XtL )-SX)**2 
AMP=SQRT(SXS/FL0AT(N) ) 
RETURN 

END 


TOTM 

fOTM 

TOTM 

TOTM 

TOTM 

to'tm 

TOTM 

TOTM 

TOTM 

TOTM 

TO TH 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTH 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

„totm 

TOTM 

TOTM 

TOTM 

jraTM_ 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

jrojM 

TOTM' 

TOTM 

"TOTM 

TOTM 

TOTM 

_TOTM 

to'tm" 

TOTM 

TOTM 

TOTM 

TOTM 

_TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 


10 
20 
30 
^0 
_^0 
60 
70 
80 
90 
100 
110 
12 0 
130 
160 
150 
160 
170 
180 
190 
200 
210 
220 
230 
260 
250 
260 
270 
280 
290 
300 
310 
320 
330 
360 
350 
"360 
370 
380 
39 0 
400 
410 
420 
430 
440 
450 
46 0 
470 
“480 
490 
500 
510 
520 
530 
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SUBROUTINE DATA 



> .... 


"return 

END 










NOTE: The following IBM Scientific Subroutines must 

be included in the program: 

CORRE, LOG, ORDER, MINV, MULTR, RANK 






Ixvl 




K R e ADR * 5 ^ — 

KPRINT=6 

“READ IN CONTROL DIGITS AS FOLLOWS AND PRINT THEM OUT, 

JDOMAN-OOMAIN OF ANALYSIS. 0 OR BLANK* PER lOD* INFREQUENCY. 
JBANO S-HJW VECTOR OF FREQUENCY OR PER I 00_B ANDS_W I L L BE OBTAINED. 

0 0 R~ B L A N K = G E N E R AT E ACCORD I NG TO L TM 1 1 S S IJ 6 S E 0 UE N TL Y " R EA tTlTT 
l=REAO IN A VECTOR OF BANOS TO BE USED FROM CARDS. 

NflAN0S*H3W MANY FRFOUENCY OR PERIOD BANDS WILL BE USED. 

NORS-THE NUMBER OF OBSERVATIONS IN THE DATA VECTOR 

JANALY-HQW THE DATA SHALL RE ANALYZED. “ 

3 OR BLANK=00 UNIVARIATE SPECTRUHt SEARCH FOR BEST K-VARIATE 

MODEL ■ 

1*SIHPLY DO UNIVARIATE SPECTRUM 

2-00 K-VARIATE SPECTRUM ACCORDING TO MODEL PARAMETER READ IN 
JDATA-T6LLS SUBROUTINE OEPVAR WHERE TO GET DATA. 

0 OR BLANK=RE40 FROM CAROS ACCORDING TO FORMAT (lOX.TFlO.O) 

^GENERATE DATA FROM MODELING PROGRAM CALLED BY DATA. THIS 

MODEL ING'PRUGRAM'WrLL ASK FOR MORE CAROS. 


page is 


OS poo® f 


Ixvil 


up oio u'o oL uju UUUU'U ouujo opooupo i I ouoou 


CRITERION FOR OPTIMIZATION OF K^VARIATE SPECTRUM^ 


JCRITR-SELECTS 

BL ANK = HULrlPlE R-SOiMREO 

i=T- VALUES 

JPRINT-PR INTOUT CONTROL DIGIT ' 

0=0ONT PRINT ANY RESULTS 

l*PRINT ONLY FINAL RESULTS 

2»PRINT FINAL AND INTERMEDIAT E RESULT S 

CORRELATIONS, ETC FOR EACH STEPSPT 

JPLOT-PLOT CONTROL DIGIT, - VALUES AS ABOVE EXCEPT FOR 3. SPT 

SPT 

SPT 
SPT 
SPT 


sy 

SPT" 

SPT 

SPT 

SPT 

SPT 

SPT 


READ (KREAOR.II 

i jpRlNT, JPLOT 
I FORMAT (1615) 


JOOMAN, J BANOS , NBANOS , JDAT A, JCR ITR , 


WRTTE (KPRINT, 103) 

103 FORMAT (*l ♦♦♦♦♦♦♦ PROGRAM PARAMETERS ■♦***•♦♦♦ ) 
~ WR I TE I KPR I NT * I I JOOMAN, j BANOS , NBANDS# NOBS , J ANALY# JOAT A, jCR ITR , 

I JPRINT,JPLOT 


‘KTOP -MAX NUMBER OIL LI ^ jy E_^OGR A^-W ILk 
-MENfS ON THE WHOLE MULTIVARIATE SPECTRUM. 
'SUCH TRIES. _ 


BE A LLOWED TO TRY 
KOUNT IS THE NUMBER 


KT0P=4 

kOUNTxO 



“SET SIGLEV, THE T-VALUE FOR THE CRITERION OF ACCEPTANCE OF A PEAK. 
-THEN DECIDE HOW TO GET VECTOR OF BANOS. 


SIGLEV=2.0 
IF (J8AN0S) 


100, lb 0,1 10 


-READ IN L'INTTS OF BANDS 
-OTHERWISE GO TO UO AND 
■ BANDLO-LOWEST FREQ. OR 
BANOEL-OELTAI FREO. ) OR 
6ANOS-VECTOR OF BANOS, 


100 

2 


120 


no 


Tr A 0 i K RE AD R ♦ 2 » BA NOl 0 , B A N D E L 

FORMAT (3F10.0) 

BANOS! I ) =BANOLO 
DO 120 LflAND=2, NBANDS 

BANDS(LBAMD»=BANOS<LBANO-l » ♦ BANOEL 

GO TO 130 

*REA't)~( KREAORi 2 ) IBANDSI LBAND) 


"vector if JBANOS MAS 0 - GENERATE BANOS. 
READ IN BANOS FROM CAROS. . 

LONGEST PERIOD, AS APPROPRIATE TO DOHA IN. 
DELTA! PERIOD) , AS APPROPRIATE TO DOMAIN. 


SPT 

SPT 

SPT 

SPT 

SPT 

_S?_L 

SPT 

SPT 

SPT 

SPT 

SPT 

SP T 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


LBAND^l , NBANDS) 


SPT 
SPT 

-CONVERT BANDS TO PERIOD DOMAIN IF ORIGINATED IN FREQ DOMAIN, SPT 

—REVERSE THE ORDER OF BANDS SO AS TO BE IN ASCENDING WAVELENGTH ORDER* SPT 

-ALSO GENERATE BAN, THE FREQ DOMAIN VECTOR OF BANOS. SPT 

— $ P T_ 

Tf~! j b om a n ) ~ i 5 b”, Ts 6 ♦ i <►0 ^ J 

00 lAl LBAND= I , NBANDS ' . .. 

8AN<L3AN0)=3AN0S(LBAN0) SPT 

BANDS! LBAND )= I .0/BANDS!LBANO) 

CALL RE V! BANDS ,NBANOS ) 

^ s ll 

"C - C E N E ft A T E” Tl N / C 0 S~P R E 0 1 1 T 0 R" M A f R lT“ A S”"P A Rt^Q^^ M ATR 1 X , X. THEN SPT 


130 
I AO 

lAl 


550 _ 
56 0 
570 
580 
590 
600 

620 
630 
640 
650 
660 
J^JO_ 
680 
690 
700 
710 
720 
_73_0__ 
740 
750 
760 
770 
780 
„T90_ 
800 
810 
820 
830 
840 
_850 
860 
870 
880 
890 
900 
J?10 
920 
930 
940 
950 
960 
970 
"930 
990 
1000 
lOlO 
1020 
1030 
1040 
lOSO 
106 0 
1070 
lOBO 
1090 
1 100 
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C*READ IN DEPENDENT VARIABLE PART OF X- THESE DATA WILL THEN BE ENT' 

INTO SUBROUTINE CORRE TO COMPUTE THE CORRELATION MATRIX OF ALL SPT 

C-VARIABLES. SPT 

C “ ■ “ - - 

150 MsNBAN0S^2 *■ 1 __SPT 

CALL PREGENIX .BANOS. NOBS.NBANOSV" " SPT 

CALL OEPVARI X.NOOS. JOATA, JDOMAN.JPRINTI SPT 

SPT 

C-GENERATE VECTOR OF VARIABLE NUMBERS FOR LATER USE« ' S_PT 

'C ■ ■ SPT 

00 155 LV»l»M SPT 

155 "CVARI LVI'LV " SPT 

CALL CORRE (NOBS »M, l» X,xaARtSTO,RX,R«BiOy T> ^ SPT 

C SPT 

C-SET TOPP, THE MAXIMUM ALLOWABLE TOTAL AMPLITUDE IN T HE FIN AL SPECTRUM SPT 
C-TOPP Is THE TOTAL AMPLITUDE 6^ PE ^ SPT 

C SPT_ 

TOPP«StD( l)».25#SfbTil " SPT 

C ' SPT 

C-PRINT OUT MATRICES ETC*, IF DESIRED SPT 

C __ SPT 

IFIJPRTnT-BI 151,152,152 SPT 

152 WRITE (KPRINT,5) __ SPT 

5 FORMATC IM A T R I X O'p C O R ~R £ L A T I 0~N $*,/•• VARIA8SPT 
ILE I IS THE DEPENDENT VARIABLE, ALL OTHERS ARE COSINE AND. SINE WAVSPT 

2ES of various periods *,/,T88,* SPT 

CALL M AT^P RT ( R ,JltL^ AR ) SP 

"C ■ SPT 

C-OECIDE ON PROPER DATA ANALYSI S PATH* S PT 

t ' ' ‘ ' ’ SPT 

151 janaly=janaly+i _ . SPT 

CO TO (200;^00,600», JANALY SPT 

C SPT 

C " SPT 

~C^4************’^***************t********t***m*¥*i*^*itt**Hf****ttt*m***tit****^** SPT 

C SPT 

C SPT 

C-STATEMENTS 200-AC0 INVOLVE COMPUTATION OF OPTIMIZED K-BANOS SPECTRA SPT 

C SPT 

C-FIRST COMPUTE UNIVARIATE SPECTRUM USING SUBROUTINE USPECT. THEN FIND SPT 

C-ALL SIGNIFICANT SPECTRUM PEAKS IN UNIVARIATE SPECTRUM* SPT 

C ^ SPT_ 

20CrTA L L U S PE C TTR , T , NOBS , M , 8 , NB a nos , S to , J P R I NT , J PL OT , PHAS E ♦ B AN , 8A NOS SPT 
1 , AMNAME,TNAME , PHNAHE , XNAME, J00MAN,8ANAME) _ SPT 

CALL PPIKR<B,T, lSAVe,NPeAK,NBANOS,SIGLEV I ' 

IF INPEAK) 201,201,202 

201 write |KPRINT,66» 

66 FORMAT! »0** • » ♦ NO SIGNIFICA NT PEAKS FOUND ♦ ♦ » »*M 

Gd“ TO '2000 ■ 


C-NOW LOOP THRU AN OPTIMIZING PROCESS TRYING TO FIND LARGEST VALUE OF 
C-EITHER T OR R-SQUARE. THE BEGINNING OF THIS LOOP IS STATEMENT 210 
C 

C^SET A VALJE OF CRIT = 0. THIS WILL LATER BE USED TO STORE THE JUST 

C-'P RE ce 01 NG VALUE OF T H E 3 P 1 1 M T Z A T I ON 'C R 1 1 6 RION^ "AL SO IN I fl AL I ZE ‘THE 


Ixtx 


-VARIABJ.E 
Counter Faft thF 


MM IS THE DIRECTION OF PEAK SLIDING AS WELL AS THE 


'number of steps- 


202 L PEAK* I 

CRITaO.O _ 

MMaO 

210 CALL KSPECT (R »T fNOfiS t H«NPEAK» ISAVEt RXtSt JPR INTt XB AR tSTPt D» RSQ » 
r^ASETSANOSTBANI 


c 

C'-tEST FOR PRINTOUT 
C 


209 


TFI JPRINT-3 ) 208f209, 209 

WRITE <KPftlNT,6) RSQ _ 

■^R^HXr r'OMULTfPL'E squared CORRELATION** #F3.4> 

WRITE (KPRINT.7I . ...... 

FORMAT! • PEftlOO/FREO. » ♦T20**AMPLlfU0E*f T40**T-VALUE*I 
DO 21A LO=l#NPEAK 
kK»ISAVe(LOI 
IF (JOOMAN) 212t2l2*213 


SpT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


212 


06=8ANDS(KK» 

_ _GO_TO 214 ^ 

213 BB*BAN(KK) 

214 WRITE (KPRINT,8I BB, B ( KK » , T f KK » 
'8 FORMAT (2F12.4iFl8.4) 


C 

C - s E L E C fTTfre R I O N for goodness of fit - CHECK FOR IMPROVEMENT 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


208 IF (CRITR (RSO .T , JCRltRtLPEAKI -CRl'T) 220t 220» 230 


A BRANCH TO 220 
•GO TO NEXT PEAK 


IMPLIES 
<221 OR 


NON-lMPROVeHENT." EftHER RESORE OLD VALUE AND 
223 1, OR TRY MOVING UP 1 222), 


220 IF (MM4U 221.222,223 

•A BRANCH TO 221 
•STEPS. THIS IS 
•AND GO TO NEXT PEAK, 


SPT 
SPT 
SPT 
SPT 
SPT 

SPT 
SPT 
SPT 

IMPLIES THAT WE HAVE BEEN MOVING DOWN FOR 2 OR MORE SPT 

THF FIRST NON-IMPROVEMENT. RESTORE JUST FORMER VALUE SPT 

SPT 

ISAVEILPEAKI^ISAVEaPEAK)-*-! 

SPT 

SPT 
SPT 

S PT 

SPT 

BRANCH TO 222 IMPLIES FIRST MOVE DOWN AND NON- I MPR0VEM6NT . TRY SPT 

SAME PEAK UP 

SPT 

222 MM=l 5PT 

I SAVE! L peak ) = ISAVFILPEAKI *2 „>PT 

SPT 


221 


MMa-1 

LPEAK=L PEAK *1 

ISAVE! LPEAK) = ISAVE(LPEAK»-1 
GO TO 300 


A 

■MOVING 


^ IF ( I SAVE( LPE AK )-NBANOS) 

225 WR ITE ( KPR INT, 

3 FORMAT ♦ * * W A R N I 

ITO THE tONGEST PERIOD (LOWEST 
2 L IMUS. • I 
GO TO 2000 


300,225,225 


N G * - JUST TRIED TO 
FREQUENCY I BAND- TRY 


SPT 

MOVE A PEAK INSPT 
WIDER SPECTRUMSPT 
SPT 

SPT 

SPT 


1680r 

1690 

1700 

1710 

1720 

1730_ 

TtVo 

1750 

1760 

1770 

1780 

1800 
laio 
1820 
1830 
1840 
_1850 
I860 
1870 
1880 
1890 
1900 
J.9t0_ 
1920 
1930 
1940 
1950 
I960 
1970 
198 o' 
1990 
2000 
2010 
2020 
_2030 
2040 
2050 
2060 
207 0 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
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C-4 BRANCH TO 223 [MPITES THAT_W^ 
C-RESTORE JU"rf''PRio'R‘ VALUE” AND GO 


223 ISAVE( LPEAK 1= ISAVE(LPEAK I-1 
LPEAK=L PE A K+ I 

I$AV£< LPEAK)= ISAVEaP6AK|-l 


HAVE BEEN MOVING UP AND IT OIDNT HELP«SPT 
t o NE XT 


PEAK. 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 


GO TO 300 


C-A BRANC H TO 230 IMPLIES I M^R 0 ^E ME NT^ C ONT I NU E . IN ^ S A^ E .D IRE C„T I ON . 

c 

230 IF (MM) 23I.231>232 


SPT 

_SPT 

SPT 

SPT 

SPT 

SPT 


C-A BRANCH TO 231 IMPLIES WE WERE MOVING DOWN. CONTI N UE. S 


231 


2<M«MM- I „ 

CRIT«CR ITR(RSO» t , JCRl TR.LPEAk I 
ISAVE(LPEAK) = I SAVE( L PEAK ) -I 


SPT 

01 RECTI ONSPT 
SPT 
SPT 
SPT 
SPT 


IF ( ISAVE( LPF AK1-U 234, 234 t 300 

234 WRITE (KPRINT,4) „ 

4 FDRMATt'O** * * ♦ W A R N I N G ♦ - JUST 

10 THE SHORTEST PERIOD ( H I GHESX FREQUENCY )_ 
2M LIMITS*) 

GO TO 2000 


SPT 

SPT 

TRl EO to MOVE A PEAK “iNT SPT 
BAND. TRY WIDER SPECTRUSPT 

SPT 

SPT 


A BR A NC H T 0 2 3 2 IMPLIES W E W E RE MOV 1 NG UP. - C ONT I NUE S AM E_. OJ R E C TJ ON. ^ 

'c ' 

232_MM=MM4-l ; 

CRIT=CRITR(RSQ,T,JCRITR,LPEAK) 

ISAVEt LPEAK ) = ISAVEtLPEAK) 


SPT 

SPT 

SPT 

SPT 

SPT 


IF USAVE(LPEAk)-NaANOSI 300*225,225 


SPT 

SPT 

C-STATEMENT 300 CHECKS TO SEE THAT WE ARE NOT OFF THE END OF THE ISAVE SPT 
0-VECTOR. IF NOT, 00 ANOTHER K-BANOS SPECTRUM WITH NEW ISAVE VECTOR SPT 

IF WE ARE OFF THE END, COMPUTE A NEW FULL SPECTRJJM USING A SPT 

FOR 'value S”0F‘ ISA TO EACH OTHER. SPT 

SPT 
SPT 
SPT 

ILPE AK-NPE AK ) 320,320.310 SPT 

a PEAK- I I 32 2,3 22,321 SPT 

t I SAVE) L RE AK ) - I SA VE( CPEAK^TI I 322,323,322 SPT 

(LPEAK-NRE AK ) 32^,210,324 „ SPT 

U SA VE(L REAK I - 1 S A V E ( L PE AK^ 1 1 I 210,325,210 SPT 

330 LL = LPEAK,NPEAK SPT 


C-VALUES. 

C-K-VARIATE MODEL. ALSO CHECK 
C-IF WE HAVE JUST MOVED INTO A ANOTHER PEAK, REDUCE NPEAK BY ONE AND 
C-PACK THE VECTOR DOWN AND START PROCESS OVER AT LP£AK»l. 


300 

320 


321 

322 
324 

323 
330 


IF 

IF 

IF 

IF 

IF 

00 


ISAVEILL-D-ISAVEILL) 
NPEA< = NPEAK- 1 


325 


TTl 


LPEAK- I 
MM*0 

CR IT-0.0 
GO TO 210 
NPEAK=NPEAK-1 
DO 33,,l.LL = LPEAK, NPEAK 
■rSA#tt L) = I SAVEILL ^1 )' 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

_5PT 

SPT 


2230 
2240 
2250 
2260 
2270 
2280 
2290_ 
2300 
2310 
2320 
2330 
2340 
_Z350_ 
2360 
2370 
238 0 
2390 
240 0 
_^410_ 
2420 
2430 

244 0 

245 0 
2460 

_247 0_ 
248 0 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2560 
2590 
' 2600 
261 0 
2620 
2630 
264 0 
2650 
2660 
2670 
268 0 
269 0 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
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lpca k-1 

HM*0“ 
CftlTaO^O 
GO TO 210 


SPT 


A BRANCH TO 310 OCCURS WHEN OPTIMIZATION OF THE K-BANOS MODEL IS 
COMPLETE. NOW COMPUTE A FULL SPECTRUM US ING TH E OPT I 

S'E ARCH"New~K-B'SPECfRU'M’>OR''P^^ AND SEE IF ANY NEW ONES TUR 
NEO UP. IF SOt GO THRU WHOLE OPTIMIZATIDNJiqUTINE AGAI^^ 


C- 

t^ooet: 

c- 
c 


310 



361 1 


9 1 


362 1 




IF (JPRINT-3)360,36l,36l 

WRITE IKPRINT.9) 

FORMAT! • 1C3MPUT ING SPECTRUM 


USING THE FOLLOWING BANOS IN MODEL*! 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

-s“pT 

SPT 

“SPT 

SPT 

SPT 

SPT 


363 


IF Ij'obMAN) 362, 362, 363 

00 36A LO^ UNPEAK 

♦CkilSAVE (L0> 

WRITE ! KPR I NT , 9 ) B ANOSI KK I 
GO to 360 

00 366 LO-l,NPEAK 


SPT 

SPT 

SPT 

SPT 

sp't 

SPT 


KK=ISAVE(L0» 

366 WRITE (K PRINT, 8) BANlKKl — --- 

360 CALL MULTI ( IS AVE,NPEAK*N8AN0S,R, t , NOBS tM#RX.B,XBAR* STD* 0»RSQf 

I PHASE , BANOS, SANI ■ 

kOUNT=KOUNT*-l 

IF (KOUNT-KTQP) 3 8A,3BA,385 

YKPRrNf,386l 


IS PROBABLY CAUGHT 
SUBMIT AGAIN*) 


385 WRITE 

386 FORMAT!' ITHE PROGRAM 
ilNPUT PARAMETERS AND 

CO TO 2000 

"3“8V CALL P P I K R ( 9 . T , I S AVE . NPE 6 K , N BANOS , S IGL EV > 

IF (NPEAK-NPE EK ) 311,350,311 

NPEAk*=NPEEK 

LPEAK = l 

MM*0 

C R I T « 0 . 0 ^ 

GO TO 210 


IN A LOOP. TRY DIFFERENT 


311 


X^STaBLE^STATE' in the OPTIMIZING ROUTINE RESUTS IN A BRANCH TO 350. 
■OUTPUT RESULTS.- , 

350 WRITE (KPR1NT,67) ^ . 

67 FORMAT! *0 MULTIPLE BAND SPECTRUM) 

CALL Oli’^PUT(BANDS,8AN,8, T , P HA S E .NBA NO S,AMNAME,TNAME«P H N AM^ 

TaXNAME , XNAME , JPR I NT , JPLOT , JpOHAN) 

WRITE (KPRINT,6) R$0 


C 

C-CHECK 

c 


FINAL SPECTRUM FOR EXCESSIVE TOTAL POWER IN SIGNIFICANT BANOS. 


CALL TOTAM(X 


NO B S ♦ N P E AK , I S A V E , M, PHASE , B ♦ SUMS ) 

i'E( SUMB-TOPP ) 391, 382, 382 

382 WRITE ( KPR I NT , 383 ) 

383 FORMAT! * IPRJGRAM FAILED BY FINDING TOO MUCH POWER IN THE 
I*,/,* SUGGEST altering PARAMETERS OF THE ANALYSIS, SUCH 
2* USE WIDER SPECTRUM LIMITS, FILTER DATA, ALTER 

381 GO TO 2000 


SPECTRUM 
AS* ,/, 
DELTA?, ETC.*) 


SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

SPT 

__SPT 

SPT 

SPT 


.SPT 

SpT 

SPT 

JSPT 

SPT 


2790 _ 
2 800 
2910 
2020 
2830 
28*V0 
2850_ 
'2860 
2870 
2880 
2890 
2900 

2920 
2930 
29^,0 
2950 
2960 
_2970_ 
2980 
2990 
3000 
3010 
3020 
_3030_ 
30A0 
3050 
3060 
3070 
3080 
_3090 . 
3100 
3110 
3120 
3130 
3140 
_3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3200 
3300 
3"^10 
3320 
3330 
3340 




Ixxii 



C************^‘***************m*******************t^***^***************** 

~C~ ■ 

c _ 

C-STATEHeNTS 400-600 JNVOLVa COMPuf AT JON FF THE SIMPLEt SINGLE BAND 
C-SPECTRUH ONLY. 

T ^ 

400 call USPECT (R *T, nobs. M.S tNBANOS, STD. JPRINT,JPLOT.PHASEtBAN.BANOS, 
I AMNAHE.TNAME.PHNAME , XNAME',JOOMAN»BANAM6r 
GO TO 2000 


Q^m*«***********^**if*******************************^m***^*m*************m 
C******^**************‘ *********************** *******nf**^if**********’^**** 

~c 

_C 

C-STATEMENTS 600-800 INVOLi/E COHPUTAtlON OF A SPECIFIC K-'BANO MODEL 
C-WJTHOUT OPTIMIZATION. 

c 

^^FIRST READ IN MODEL PARAMETERS. BUILD VEC T OR OF B ANDS I N MOD EL. THEN 
C-COMPUtE SPECTRUM. " ■ 

_C 

600 READ IKREAOR.I) NPEAK 

READ IKREADR.2) I PEAKS I L PE AK I , LPEAK® 1 * NPEAK I 

DO 610 L PE AK=1. NPEAK 
00 620 LBANO-l.NBANO 

IF (PE AKS(LPeAK)-BANOSTLBANOn ‘620,630,620 
630 ISAVE< LP6AK l=L8AND 
GO to 610 
620 CONTINUE 
610 CONTINUE 

CALL MULTI USAVE, NPEAK, NBANOS.R.T, NOBS. N.RX.B.XBARtSTO.OtRSQt 
1 PHASE, BANOS, BAN) 

WRITE IKPRINT.5) _ ___ _ 

CALL OUTPUT ( BANOS , BAN , B , T , PHASE , NBANbSVANAME,TNAME',PHMAME, 

A N A M E , X N A M E P^l NT , JPLQT , JOOMAN ) 

WRITE IK PR I NT, 6") R5Q 

2000 STOP ' 

END “ " ‘ . 





u u u u u 


9 9 


R6V 

subroutine REVIX,N» rev 

DIMENSION XI 501 " REV 

. REV 

REV 

_* : ♦ REV 

♦ SUBROUTINE REV ' ♦ RPV 


REV" tAKES A VECTOR AND REVE^ "♦ 

* » 


J-N^l 

MIp.N /2 

00 ioT»r 

J«J-l 
SAVE=X( LI 
X<L MX( Jl 
10 XUT*SAVE 

RETURN 

END 



c OCPV 10 

— C ' ' “ ' oepv 20 

SUBROliTINe OEPVAR(X,NOBS* JDATA, JOOMAN* JPRINTI OEPV 30 

DIMENSION Xm ' OEPV 40 

COMMON kRF&or OEPV 50 

OEPV 60 

Q OEPV 70 

CT* ♦ OEPV 80 

C * SUBROUTINE 0 E P V A R _OEPV 90 

m OEPV 100 

C ♦ DEPVAR EITHER READS IN OR GENERATES <VIA SIMU L AI, A D EPEND EN T » OEPV tlO 


C ♦ INTO THE FIRST NOBS CELLS OF X, A MATRIX OF VARIABLES STORED AS A* OEPV 130 

C ♦ VECTOR OF COLUMNS. SEE SUBROUTINE PAEGEN FOR REST OF X. • OEPV 140 

C * ♦. OEPV 150 

*#♦*♦'*♦♦♦♦♦♦♦♦**♦*♦♦***♦*• OEPV 160 

C ^ O EPV I7 Q 

C-IF JOATA IS 0 - READ FROM CAROS. (10X»7F10.0) ^EPV 180 

C-IF JDATA is 1 - CALL SIMULA WHICH GENERATES DATA USING A MONTE CARLO DEPV 190 

C-SYSTEM ALONG WITH DETERMINISTIC DATA. SIMULA READS CARDS. OEPV 200 

C DEPV 210 

IF TjOATAI 10.10.50 “ ~ OEPV 220 

10 READ (KREADR.U ( X ( L I »L = I . NOBS J OEP V 230 

I FORMAt ( lOX.'TF LO.OI OEPV 240 

__ GO TO 2 0 OEPV 250 

~“5CT CALL SIMULA! X.NOBS.JDATA/jOOHAN.JPRlNf r OEPV 260 

20 RETURN OEPV 270 

END ” ~ OEPV 280 


iSSSS'fipl^ 
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nnoorJ ! i r»in oo o n o n no o n o o n o rvo r>o r»o on | o|o 


SUBROUT INE PREGENI X, BANOS. NOBS# NBANDS) 
OIHENSION XlUfRANDSU) 


♦ ♦♦♦♦♦♦♦♦ •♦♦♦'♦**<i***^******»****%*^4i ♦♦♦♦♦*♦*♦♦♦♦**♦♦*♦♦**♦♦♦* ♦♦**** 


SUBROUTINE P R E G E N 


PR E GEN GENERATES THE SINE AND COSINE PREDICTOR WAVES AND STORES 
THEN IN MATRIX X. MATRIX X IS THE DATA MATRIX HAVING NOBS ROWS 
AND M COLUMNS . HERE N06S= THE NUMBER OF GBSERVATfONS A'nO 
BANDS I S A. yec TOR O F P ER 1 00 VALUES I_N_THE SPEC 

MATRIX X IS STORED IN VECTOR MODE. SEE IBM SYSTEM/360 SUB- 
ROUTINE PACKAGE (360A-CM-03X) VERSION III, PROGRAMMER'S MANUAL, 
PUBLICATION NUMBER H20-0205-3. PAGES 3-6* ESSENTIALLY EACH CDL~ 


UHN OF MATRIX X IS STRUNG ENO-TO-ENO INTO ONE LONG VECTOR. 

t ME F I R S T NOBS V A L UE S OF VECTOR ( MA f R I X } X ' W I LL ~ B E T HE ’ D E P E NOE N T 
VARIABLE ISEE SUBROUTINE DEPVAR). THE NEXT NOBS POINTS ICOLUMN 
TWO) WILL CONTAIN COS WAVE, BAND I, THE NEXT NOBS POINTS (COLUMN 
THREE CONTAINS SIN WAVE, BAND I, ETC. 




PREO 

PREO 

PRED 

PREO 

J>REp„ 

PREO 

PREO 

PREO 

PREO 

PREO 

PREO 


-CONVERT TO RADIANS/OBSERVATION AND LOOP THRU NBANOS TIMES._ 


_L»J40BS 

(To to LB“AND=‘t /NBANOS 

AFRE0=6.2B313/6AN0S(LBAN0) 

EF=-AFREO 

00 20 LOBS=l,NOeS _ _ _ 

■L = L + 1 

EF-EF^AFREO 


20 X(L)=C:OS(EF> 
EF=-AFREO 
00 10 LOBS-l,NOBS 
L*L*1 

■ EF=EF^AFREQ 
10 X(L»-SIN(6F» 


•X CONTAINS THE VECTORS 3F SIN/COS WAVES STRUNG ENO-TO-ENO. THIS IS 
■EQUIVALENT TO A MATRIX WHERE EACH SIN OR COS WAVE FORMS A COLUMN OF 
■NOBS LENG'^H. 


RETURN 

^ENO 


PREO 
PREO 
PREO 
PREO 
PREO 
PREO 
PRED 
PREO 
PREO 
PREO 
PRED 
_PREO 
PRED 
PREO 
PREO 
PREO 
PREO 
_PR£P_ 
PREO 
PREO 
PRED 
PREO 
PRED 
PREO 
PREO 
PREO 
PRED 
PREO 
PREO 
^REO 
PRED 
PRED 
PREO 
PREO 
PREO 
^REO 
PRE 0 


10 
20 
30 
40 
_50 
60 
70 
80 
90 
100 
Jto. 
120 
130 
140 
ISO 
160 
no 
180 
190 
200 
210 
220 
_230 
240 
2S0 
260 
2T0 
260 
29 0 
300 
310 
320 
330 
340 
350 
360 
370 
380 
39 0 
400 
410 
420 
430 
440 
450 
46 0 
470 
48 O" 






« 
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n n 


SU6RQUT INE S IHUl A( X* N, J , JOOHAN, JPRINTI 
DIMENSION XU ) ,An30» ,8ANDI I ) 

COMMON KREAOR,J<Pa;NT 


SIHL 

SIML 

SIHL 

SIML 

SIML 


10 
20 
30 
40 
_50 
60 
70 
80 
90 
100 
llO 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
230 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
49 0 
500 
r • Q 

520 

530 

540 


V ■ - • - - - ■ - ^ 

♦ SUBROUTINE S l,_OA 1—- * 

♦ - - - ^ 

♦ SIMULA IS A PRIMATIVE MONTE CARLO AND SIGNAL GENERATOR FOR CR 6AT>* 

♦ ING AATIFICrAL DATA. GAUS'SIAN DATA ARE READ IN FROM cARDS AND ♦ 

♦ A NUMBER OF SINE WAVES OF VARIABLE FRE QUENCY_ AND AMPLITUDE ARE ♦ 

4 ADDED. ‘ "y 4 

♦ 4 


SIML 

SIML 

SIML 

SIML 

SIML 

SIML 


C-iERO VECTOR X AND PR”i¥t” HEAOTf^ 

c 


DO 100 L=l.N 
JOO X(LI=0.0, 

~IF (JPBINT) UOrllO, 
120 WRITE (KPRINT.12U 


120 


SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 


121 FORMAT ( • IM 0 N T E 


CARLO 


SIMULATED OATA»,/**0»J 


READ IN MODEL SPECIFICATIONS - GENERATE SINES AND ADO NOISE 

NSIN-NUM8ER OF SINE WAVES IN DATA 

BAMO-FREO. DR PERIOD OF THE LSI NTH SINE WAVE ‘ 

SD“ST AJ^DARD^EFl^I 01^ Q F T HE LSINT H SINE WAVE. (AMPLITUDE! 
SON-StANOARO OW. OF GAUSS'lAN NOISE 


SIML 

SIML 

SIML 

SIML 

SIML 

SIML 


110 READ (KREAOR.U NSIN 

_ I FORMAT (15) 

DO 10 LSIN=1,NSIN 
READ (KREA0R.2) BANO(l) 


SD 


30 


FORMAT (2F10.0) 

IF (J00MANI20. 20*30 _ 

BAN=BAND(l) 

BAND! II *UO/BANO( U 

GO TO 21 

20 BAN=l.O/BAND( I I 

21 IF (JPRINT) 25/25T22 

22 write (KPR1NT.3) LSIN* BANOdl# BAN*SD 
3 FORMATS BAND*. 13** - PER 100** *F10.4, * , 

1TU0E*»,F10.41 

25 AFR£a=6.233l 8/BANOU ) 

EF=“AFREO 

DO 26 'L^rt N 
EFsEF4-AFRE0 

26 A(LI*SIN(FF)41. 41421 
DO 10 L = UN 

10 X( L)=X(U*A(L I^SO 
C 

■C^REAO"Ifi NLnSE' SO AND NOTSE’ C ARDS 


SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

■ SIML 

SIML 

_ SIML 

" SIML 

SIHL 

SIML 

SIML 

FREQUENCY** »F10. 4, *, AMPLISIML 

SIML 

SIML 

SIML 

SIML 
SIML 
SIML 
S T ■■■. 
SIML 
SIML 
Sf “L 


Ixxvli 


SIML 550 


READ iKREA0R,2l SON 

IF(SONI L000,1000,40 

‘40 iF UPRINTJ 50,50#4l 

^IMRITE (KPRINT.5) SDN 
5 FORM4T<»ORANOQM GAUSSIAN NOISE 
50 READ (KReADR,6l (A(UltL*ltN» 


“AM PL I tub Ei • , F I 6^ 4 ) 


b FORMAT ( 4X , To F”^ 4T” 
_ 00 60 L=l .N 

60 X4U=A( LI*SON^X(L » 
JOW RETURN 
END 


SIHt 

SIML 

SIMt 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 

SIML 


560 
570 
580 
590 
600 
610 
620 
630 
64 0 
650 
660 
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SUBROUTINE MATPRTIR*M»LVAR) 

' "DIMENSION R< I » ,XU5) ,LVARI U 
COMMON KREAOR* KPRiNT 

^ 

~C~* ■■ ■ ♦ 

C_* SUBROUTINE M A T P R T _• 

c • 

C ♦ MATPRT PRINTS OUT THE LOWER TRIANGULAR MATRIX X WHERE X IS STOR- « 


• EO IN MODE 1. SEE REFERENCE LISTED IN SUBROUTINE PREGEN* * 

C • _ * 

C ♦♦♦♦ ♦♦♦♦*♦♦♦♦ ♦ i* * * «« ♦Vi* ♦"'Ml ’♦♦’♦ ♦♦'♦ * ♦ ♦♦•♦♦♦♦♦ ♦ ♦ ♦ ♦ • ♦ ♦ ♦'♦ ♦" 

C 

C-LOOP THRU AS MANY TIMES AS NEEOEO' tO PRINT MATRIX; EACH TIME 

C-PRINTING 15 COLUMNS BY NRO ROWS. 


Nl-l 

N2*15 

UO IF IN2-MI 100,100.10 

iO N2=M 

100 WRITE (KPRINT, II ( LVAR(LV) «LV«N1 • N2) 


1 FORMAT! •0*,2X.15I7) 

WRITE (KPRINT.2) _ 

2 FORMAT I *0* I 


2 F 
C 

C-PRINf 


"one row at 


time, HAtRix 


IS MODE 


UPPER TRIANGULAR. 


C-THEREFORE REVERSE SUBSCRIPTS TO MAKE LOWER TRIANGULAR. 


DO 20 LROW-Nl,M 
IF ILR0W-N2I 21,21.22 

21 LlMITsLROW 
GO TO 2 3 

22 LIMIT=N2 

. 23 L=0 

DO 30 LCOL^Nl, LIMIT 
" L*L+l 

CALL LOCUCOL,LROW,K,M,M, n 
30 X(LI=R(KI 

20 WRITE (KPRINT, 31 (LROW , I X( K ) ,K=i I, L M 
3 F ORHA T U X , I 3,2 X ,T5F7 . 31 

Nl.= Nl*l5 ♦ _ 

N2 = NUl^ " 

IF (M-Nll 1000, 110,110 „ 

"1000 RETURN ~ “ 

END 



c USPT 10 

^ - y 20 

subroutine USPECT(R,T,N ,M,B#NBANOS*STO*JPRINT, JPL0T*PHASE*BAN* USPT 30 

r BANOS* AMNAME, TNAM6,PHN4ME, XNAME,JD0MAN*8ANAMe» USPT 40 

01 MENS ION R{ 1 } ,T( 1 3(1? »STOJ I) *P HASEI II* BAN! 1 ) * BANDS I I ? » AHNAHg I DUSP T 5 0 

FTTn a P HN A met 1 ) , X NA ME U ) , B A NA ME fU USP T 6 0 


jj ■ " USPT 80 

C USPT 90 

~ ' ■ ♦ USPT 100 

C SUBROUTINE U $ P E C T USPT I IQ 

Z~-i • USPT 120 

C * USPECT COMPUTES THE UNIVARIATE SPECTRUM BY MULTIPLYING THE CORR- • USPT 130 

C“*“ ELATION OF THE TIME SERIES WITH A SIN OR CDS PREDICTOR BY THE * USPT 140 

C * RATIO OF STANDARD DEVIATIONS OF THE TWO* THIS IS DONE FOR EACH * USPT ISO 

C * PREDICTOR. LET R=THE CORRELATION OF F(T» WITH A S IN OR COS WAVE ♦ USPT 160 

C * OF SOME ARBITRARY WAVE LENGTH. LET S ^iTH E ST ANDAR D DE VIAT ION * U SPT I7Q 

” OF“T'Hri’lME SERIES AND .TOY^THE StANOARO DEVIATION OF THE PREO- * US>T 180 

C * ICTOR WAVF. THEN A(F) * R4(S/.707) _ • USPT 190 

C *' HHERE MF) = THE AMPLITUDE OF THE COMPLEX SPECTRUM AT F. THE SIN • USPT 200 

C AND COS C0MP1NENTS ARE COMBINED AT EACH WAVE LENGTH AND PHASE * USPT 210 

C * ANGLE IS COMPUTED. T-VALUES ARE ALSO COMPUTED FOR EACH BAND. • USPT 220 

C ♦ R-UPPER triangular CORRELATION MATRI X, IN STORAGE HO PE 1. SEE * U SPT 230 

C ♦ IBM PUBLICATION REFERENCE 6 IN SUBROUTINE PREGEN. UNPUT) * USPT 240 

C • T-VECTOR OF T-VALUES. (OUTPUT! ... * USPT 250 

C * NOBS^NUMBER OF OBSERVATIONS. (INPUT! ♦ USPT 260 

C * M*TOTAL NUMBER OF DATA VECTORS. M=NBAN0S*2 ♦ 1. ONE VECTOR FOR * USPT 270 

C * the dependent variable and two for the SIN/COS WAVES AT EACH WAVE* USPT 200 

C * LENGTH. (INPUT). 0=VECTaR OF SPECTRUM ESTIMATES. (OUTPUT ! * U SPT 290 

C“* N BA N DS = L E N G T H 0 F V EC TO R S t , 0 A NO PH A S E . UN PUT I ♦ U S P T 30 0 

C * $TD-$TANDARO DEVIATIONS OF ALL VARIABLES, (INPUT! * US»T 310 

C * JPRINT=PR*INT CONTROL - SEE MAINLINE. (INPUT! * USPT 320 

C ♦ JPL0T=PLGT control - SEE MAINLINE. (INPUT! ♦ USPT 330 

C“* PHASE-VECTOR OF PHASE ANCLES* (OUTPUT! * USPT 340 

C * * U SPT 350 


USPT 350 
USPT 370 
USPT 3B0 
USPT 300 
USPT 400 
J^SPT 410 


00 10 L=2*M • USPT 42 0 

CALL LOCI l.L.KOOL.M.M.l) USPT 430 

V ^ 10 B(L)-R(KOOL) USPT 440 

C USPT 450 

C-LDOP THRU THE SPECTRUM COMPUTATION NBANO TIMES^ USPT 460 

C USPT 470 

FACY^S 0 R T ( F L Q^T\ N- 2 ! ! USP T 4 8 0 

. ^ 00 20 Ls=2,M,2 _ _ USPT 490 

" K-L+l. ' ’ USPT 500 

LBANO = LBAND*' I USPT 510 

C USPT 520 

C-CO=COS COMPONENT OF^SPECTRUM ESTIMATE USPT 530 

C^ST= SIN COMPONENT OF SPEC T RUM” ESt I MATE USPT 5*0 


******* ********* *ti************************^<t>*^**^**^*^** 

LBANDsO 


C-EXTRACT VECTOR OF CORRELATIONS BETWEEN DEPENDENT VARIABLE AND EACH 
C-PREDICTOR VARIABLE. 

C 



* 


Ixxx 


C-TC«T-VALUE FOR COS COMPONE NT USPT S 50„ 

C*tS»T-VACOE FOR S IN COMPONENT USPT 560 

C OSPT 570 

tO»BCL »*STO( U/.T07 USPT 580 

SI-B(K)»STO< I >/.707 MSPT 590 

TC~B<L M‘FACT/SQRT< 1.0-B(L)»*2T WSPT 600 

TS«B( K ) ♦FACT/ SORT ( l.O-B( K )^*2) USPT 6 10 

^ — USPT 620 

C-COMPUTE SPECTRUM AMPLITUDE, T-VALUE ANO_PHASE_ ANGLJE OSPT 630 

C ■ ' ' ' USPT 6^0 

BILBANO»=SORT(CO**2 ♦ Sl**2) OS^T. 65 0 

T(LBAND»=S0RT(TC^*2 ♦ TS^*2I USPT 660 

20 PHASE(LfiANO>*ATAN(CO/Sn OSPT^6 7_0^ 

_ — — — — USPT 680 

C- PR I NT AND/ OR PLOT IF DESIRED O^PT 69 0 



JPR«JPRINT-1 ' _:_0 >PT.7^0 

JPL*JPLOT-l USPT 750 

IF CJPR» 50,50,60 USPT 76 0 

50 IF < J PL ) 100 0,1000,60 USPT 770 


60 WRITE <KPR1NT,U USPT 780 

I FORMATfO S I N G L E BAND SPECTRUM*) USPT 790 

CALL OUTPUT (BANOS, BAN, B, T , PHASE ,NBANOS, AMNAME, TNAMEtPHNAME 

iVBANAME,XNAME, JPRINT, JPLOT, JDOMANI USPT 810 

1000 RETURN USPT 8 20 

USPT 830 




SUBROUTINE PPl KR < B f T t ISAVE # NPE AKfNBANOS t S IGLEV I 

dimension B( I > ,TUM SAVE! l» " 

* 

• SUBROUTINE P P I K R~ "" 

* 

"P returns THE VECT OR INDEX OF' ALL “ S I GN I F YC A N T SP EC TR UM P E A K S . 

* S - 1 N PiJT VECTOR OF AMPLITUDE SPECTRUM ESTIMATES 

"* r=TNPUf VECrCJR OF T-VALUTS" FOR EAYtH SPECTRUM ESTIMATE 
_♦ ISAVe-OUTPUT VECTOR OF INDEX NUMBERS OF SIGNIFICANT PEAKS 

* NPEAK-OUTPUT SCALAR - NUMBER OF PEAKS FOUND 

.• NSANOS- INPUT SCALAR - NUMBER OF BANOS IN SPECTRUM 

* SIGLEV-INPUT SCALAR - CRITERION WHICH T-VALUE FOR A PEAK MUST 

EXCEED IN ORDER TO BE RETAINED 


* 

♦» *♦♦*»»♦ »* * * * ************* * ************ ♦♦.*♦*♦♦*♦♦♦ ♦♦♦♦♦♦•«♦♦♦♦*♦•♦* 

SEARC H F OR P E A K^. MTA^I, N _^S I GN C ANn ONES t _ BU I L 0 I S AV E . 


NPEAK»0 

DO 10 L=2*N8ANDS 
IF (BCL-D-BILM 
20 IF (BtL>l)-BtLn 
40 IF (TIU-SIGLEVI 
50 NPEAK»NPEAK«-l 
_ ISAVEINPEAK»=L 
To CONfiNUE 
RETURN 
END 


20 , 10,10 
40, 10, 10 
10,50,50 


PPIK 
PPI K 
PPIK 
■pPlK 
PPI K 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 
PPIK 


PPIK 

PPIK 

PPIK 

PPIK 

PPIK 

PPIK 

PPIK 

_PPIK 

PPIK 

PPIK 

PPIK 

PPIK 

ppU 

PPIK 

PPIK 




Ixxxti 


o o n o r»p o r» o n o n OiO r>ln rJr« ! n o 


ctat 

CRIT 

FUNCTION CRITR(RSO«r«JCRlTR«LI CRIT 

01 HENS ION T<U CRIT 

CRIT 


♦ _ • CRIT 

♦ FUNCTION C R I T ft ♦ CRIT 

♦ ♦ CRIT 

♦ CRITft SELECTS A CRITERION FOR THE SPECTRUM OPTIMIZATION V'CRIT' 

» SHIFTING. SELECTS EITHER R-SOUARE OR THE APPROPRIATE T VALUE. » CRIT 


♦ * CRIT 

. _ CRIT 

SELECT APPROPRIATE VALUE FOR CRITERION. CRIT 

RSQ- INPUT SCALAR - VALUE OF R-SOUARE CRIT 

T-INPUT VECTOR - T-VALUES ’ CRIT 


JCRITR-CONTROL DIGIT FOR SELECTION OF EITHER RSQ OR TIL) AS CRITERI ONCRI T 

0«RSQ AS CRITERION ■ CRIT 

l«T-VALUE AS CRITERION CRIT 

L-INDEX VALUE OF VECTOR T CRIT 

CRIT 


IF IJCRITR) 10, 10,20 CRIT 2 30 


10 CRITR=RSO CRIT 240 

GO TO 30 _ CRIT 250 

20 CRITR=*T(L) CRIT 260 

3^ RE TURN CRIT 270 

END CRIT 280 



Ixxxiil 


c 



KSPT 

, 10 

c 



KSPf 

20 

SUBROUTINE < SPECT (R , T * N, M ,NPE AK, I SAV ,RXtB*JPRlNT 

,XBAR,STO, 

D,RSO, KSPT 

30 

IPHASE, BANOS, BAN) 



KSPT 

40 

OIHENSION ISAVEdOOl ,R ( 1 ) ,T( 1) ,RX( U « BI U ♦XBARIll ♦ 

STOI 11 

,0111, KSPT 

50 

TphASEI n, BANOS ( t), BAN 1 U , ANSI 10), 

ISAVI 11 ,RY(100I KSPT 

60 

COMMON KREAOR,KPRINT 



KSPT 

TO 

"C-CONVERT ISAV TO I SAVE, FROM FREQ BAND VECTOR TO VECTOR OF 

SIN/CQS PAIRKSPT 

80 

C _ 



KSPT 

90 .. . 

98521 FORMAT I»0JKL**,I5) 



KSPT 

100 

K»0 



KSPT 

110 

00 100 L»l,NPEAK 



KSPT 

120 




KSPT 

130 

I SAVE I K I* ISAV (L 1*2 



KSPT 

140 

K*K^l 



KSPT 

150 

100 rSAVE<K)s=iSAVE{K-Uf 1 



KSPT 

160 

C 



KSPT 

170 

1- SELECT SUB-MATRIX OF PREDICTORS ETC FROM R, ACCORDING 

TO 

ISAVE. 

KSPT 

180 

C 



KSPT 

190 

CALL 6RbERIM,R,l,K , ISAVE, RX ,RV ) 



KSPT 

200 

C 



KSPT 

210 

C-iNVERTK-ORDER MATRIX OF PREDICTOR iNTERCORRELAtlONS, 

RX 


KSPT 

220 

C 



KSPT 

230 

CALL MlNV(RXtK.DET*B,T) 



KSPT 

240 

IF I.JPRlNT-3) 300,400, AOO 



KSPT 

250 

400 WRITE (KPRINT,7» DET 



KSPT 

260 

7 FORMAT! '0 DETERMINANT »*,E30»15) 



KSPT 

270 , , 

c 



KSPT 

280 

C-COHPUTE REGRESSIONSISPECTRUM) FOR NPEAK SUED MODEL. 



KSPT 

290 

c 



KSPT 

300 

300 CALL MULTRI N,K ,X8AR, STD, 0,RX,RY, ISAVE, S, SB, T, ANSI 



KSPT 

310 

RSQ»ANSt 2)**2 



KSPT 

320 

c 



KSPT 

330 

“C-CONVERT SPECTRUM to AMPLITUDES BY COMBINING SIN /COS. 

PUT 

INTO 

0. KSPT 

340 

C-ALSO CONVERT T. 



KSPT 

350 

c 



KSPT 

360 

NFRE0*0 



KSPT 

370 

00 200 L*1,K,2 



KSPT 

380 

NFR£0“NFRE0^1 



KSPT 

390 

PHASE! NFREO) -AT AN(BI L ) /B (L* IM 



KSPT 

400 

B!NFREO)»SORT!B(L)**2 * BIL+l)**2) 



KSPT 

410 

200 T(NFRE0)-S0RTITtU**2 * TIL*1)**2) 



KSPT 

420 

10 RETURN 



KSPT 

430 

END 



KSPT 

440 
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ft 



SUBR3UT1NE «UL Tl< ISAVE#NPeAK,N8ANOS,R,T,N*M,RX.B»XBAR*STO.O.RSO, 

I PHASE, BANOS* BAN) 

ye ( I J , RCI),T( n_, RXJ UfB(l),XBARn),STOa),0(l).PHASE 

I < n,BANOS( I ),3AN( n *KSAVEC20)tS(20l,fS>6C<20) *TPHl"20) 

COHMON kreaor,kprint __ 

^ * ** ******* * ♦ **. * * * ♦ ♦♦♦♦♦♦** ♦♦ ♦♦♦♦♦♦♦ ♦»♦*»»»♦♦ ♦♦» «««««« •««««« 

C * SUBROUTINE MULTI « 

— 

C * MULTI COMPUTES THE FULL SPECTRUM OF A TIME SERIES USING A K>1 ♦ 

C * BANOS MODEL WHERE K=THE NUMBER OF PEAKS I N THE CPTQMI Z60 • 

T'Hf extra band IS THE WAVELENGTH IN QUESTION, THE OTHERS • 
C * ARE THE SIGNIFICANT PEAKS. THE VARIANCE ACCOUNTED FOR BY ♦ 

C ♦ THE signific ant PEAKS IS THUS *ACCOUNTEO FOR* IN EACH ESTIMATE. • 

WH^EN the current BAND IS ONE OF fHE SIGNIFICANT PEAKS. A SKIP * 

C * . o ccurs . VA R I A at E $ DE F I NED J N MAI NL I NE ^ 4> 

C *♦*♦*♦**♦♦♦ ♦*♦♦* **< » »»*♦»♦»»♦»»»»»»»♦» »»»»»»»!»»» 

C 

t-LOQP THRU ONCE FOR EACH FREQ BAND. NFMOO IS THE NUMBER OF SIGNICIC'^ 
C-CA NT PEAKS IN MODEL PLUS A CURRENT BAND. 

c 

NFMOOaNPEAK*- 1 

DO 310" LFsUNBANOS " 

C 

C-C ON STRUCT NEW I SAVE VECTOR 

C 

^ CALL 3 U I l6 ( K S A VE , I S A VE .LF , NFMOO^ ™ 

pSKlP~ I r THIS BAND I $ IN THE MOO EL ALREADY 
^ • 

I F 1 1 S KTP ) 320 rr2"073^ 

:-BRANCH TO 330 IMPLIES THAT THE CURRENT 8AN0 IS dNlTl)F THE SfCN IF I’CANT 
>PEAKS. 00 AN NPEAK SPECTRUM AND PICK OUT CURRENT BAND. 

330 CALLKSPECT ( R , S ,N,M, NPEAK , ISAVE,RX,TSPEC,0,XBAR, STO.O.RSO.TPH. 

I BANDS, BAN) : 

:-TSPEC CONTAINS ALL MODEL ESTIMATES. PICK THE ONE FOR LF^ 

00 331 LLZ*l, NPEAK ' 

^I£^(LF«ISAVF ( LLZ )) 331,333,331 

333 BI LF) = TSPEC < LLZ ) “ ' — 

PH4SE f LF )=TPHfLLZ) 

TILF)=S(LLZ) ~ " 

331 CONTINUE 

GO TO 310 - 


MULT 

10 

MULT 

20 

MULT 

30 

MULT 

40 

MULT 

50 

MULT 

60 

mul t 

TO 

mult 

80 

mult 

90 


.-BRANCH TO 320 TMPLTE5~tHAt THIS IS A FULL NFMOO MODEL. 00 SPECTRUM 


► MULT 100 
» MULT n o 
MULT 120 
' mult 130 
MULT 140 
MULT 150 
MULT 160 
MULT 170 
MULT 180 
MULT 190 
MULT 200 
MULT 210 
MULT 220 
mult 2 30 
MULT 240 
_MULT 250 
MULT 260 
_MULT 270 
MULT 280 
MULT 2 90 
MULT 300 
_MULT 310 
MULT 320 
„MULT 330 
MUL T 340 
_fiykJ„350„ 
MULt 360 
mult 370 
MULT 300 
MULT 390 
MULT 400 
_MUL T_4 1 0 
MULT 420 
MULT 430 
MULT 44 0 
MULT 450 
MULT 460 
_MULT_470 
MULT 48 0~ 
MULT 490 
MULT 500 
MULT 510 
MULT 520 
JIULT 530 
MULT 5^,0' 


Ixxxv 



C-AND PICK OUT CURRENT BAND. 


320 CALLK3?ECTIR,5»N,H,NFH0D»KSAVE,RX,TSPEC,0»XBAR,ST0*0»RSQ*TPH, 

I BANOS* BAN J “ ■ ' 

00 321 LLL=l,NFMOO 
IF(LF-KSAVE(LLD) 321*322*321' 

TUF»^SILLL) 


322 


"321 

310 


"STTF ) = T5 P E C < L L L j 
PHASE ( LF ) =TPH(LLL ) 
CONTINUE 

continue 

RETURN 


_MULT 

MUct 

MULT 

MULT 

MULT 

MULT 

MULT 


MULT 

MULT 

MULT 

MULT 

MULT 


550 
56 0 
570 
580 
590 
600 
A* 0 
620 ' 
630 
64 0 
650 
660 
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C _ __ _ BllD._lO 

C ' BILO 20 

SUBROUTINE BU IL DIKSAVEf I SA VE,L F* NFMOO , I SK IP .NPEAK > _...BILO 30 

DIMENSION KSAVEU I *ISAVE( U *ASAVEI30) ♦RANCI30) BILD ^0 

C 8 T LD S O 

C «*«««* B I LO 50 

C_* _ » SILO 70 

c SUBROUTINE ft U I L D BltD SO 

* BUILD CONSTRUCTS VECTOR OF RANKED FREQ BANp_NUMBERS_^^^^ •..BILO, 90 

C“i' PARTICULAR MODEL AND A CURRENT FREQUENCY, • BILO 100 

C ♦ KSAVE - OUTPUT VECTOR TO BE BUILT . ♦ BILD I IP 

C ♦ ISAVE - ’input VECTOR CONTAINING NFREQ SIGNIF PEAKS ♦ BILO 1^0 

C » LF - INPUT SCALAR - NFW FREQ BAND TO BE ADDED (CURRENT ONE) • BILD 130 

C* NFMOO - INPUT - NUMBER OF FREQ BANDS IN NEW MODEL ♦ BILO 140 

C ♦ ISKIP - OUTPUT CONTROL DIGIT TO KEEP FROM COMPUTING JEST IHATE„ __*.BILO ISO 

C • ‘ OF ONE OF THE OLD PEAKS AGAIN ♦ BILO 160 

C * NPEAK - INPUT - NUM BER OF FREQ BAN DS IN OLD MODEL » BILO 1 70 

C • ♦ BILD 180 

C RILO 190 

~(T~ ’ ~ " BILO 200 

C-SEE IF LF IS ALREADY (N MODEL BILO 210 

C BILO 220 

00 10 L = 1 _AK BILD 2 30 

IFILF - ISAVeiD) 10,99*10 BILO 240 

10 CONTINUE BILO 250 

C BILO 260 

C-STORE VECTOR OF ISAVE PLUS LF IN TO ASAVE AN D RANK IT, RANK STORED IN BILO 270 

C-VEC TOR RANG BILO 280 

_C BILD 290 

DO 20 L=l*NPEAK BILO 300 

20 ASAVE( U=I SAVEIU ^ BILO 310 

'■'ASAVE( NFMOO )=LF BILO 320 

CALL RANK! ASAVE, RANCtNFMOD) BILD 330 

C ■ “ " BILO 340 

C~6UIL0 KSAVE OF RANKED FREQ BANDS BILD 3 50 _ 

______ 

00 30 L = l, NFMOO 

KOOL“RANC(L) BILD 380 

30 KSAVE( KOOL>=ASAVE(LI BILO 390 

ISKIP=0 ” BILO 400 

GO TO 100 _ B ILO 410 

99 ISKYPsI ' BILO 420 

100 RETURN _ BILO 430 

END ’ ’ ’ ’ ^ 
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n 0:0 


subroutine OUTPUTIBANOS,BAN*BtT,PHASE#NBANOS* AMNAME,TNAME»PHNAHe, 

IBANANE* XNAM 6 , JPR (NT, JPLOT • JOOHANI 

DIMENSION BANOSdUBANHUBU l,T(U,PHASE(l),AMNAMEmtTNA«EI II # 

TPHNAMEI U f BANAMEI 1 ) /XNAMEfi I 

COMHON KREADR.KPRINT 

**♦****♦♦♦♦****♦♦*♦♦♦**♦*♦♦♦♦***♦**♦*♦♦»•♦♦♦**♦♦*♦** I* ♦♦#***#*^***^^^* 

* SUBROUTINE OUTPUT * 


and/ or plots spectra with documentation. 


_?1^^*”**** **•****♦*****♦*♦♦♦♦***♦♦** ♦♦♦♦**♦**♦♦**♦♦*♦♦♦♦•♦*♦♦*•**»♦♦ 


c 

C-PRINT OUT DOCUMENTATION 

T^ino scale 


VALUES AND 


_ANO^ 
PLOT 


VECTOR ORDER IN APPROPRIATE OOMAINI 


OUTP 

OUTP 

OUTP 

OUTP 

_PUJTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

-OUT^ 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 


T0PI*6. 28318 

BTDPI»-TOPI 

BMAX-B(l) 

00 6^ L8AND=2,NBAN0S 


IF ( B M A X - B ( L B AN D ) I fi 8 * 6 A , 6 ^ 
68 BMAX=8« LBANO ) 

64 CONTINUE 
THAX=T( I I 

DO 66 LBAND = 2,N8ANDS 

IF ITMAX-TI LBAND) I 67,66,66 


OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 


67 TMAX='T<LBANOi 

_66_C0NTINUE 

IF (JO OMAN) 60,60,50 

•REVERSE VECTORS IF FREQ OOMA IN' TS' USE 0 • ^ I N t”~RE VE R AG AI N . 


50 


OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 


CALL R E V I B , N BAN 0 S ) ~~ 

CALL REV(T,NBANOSI 

CALL REV(PHASE,NBANDS) : 

IF (JPRINT) 52,52,51 

51 WRITE (KPRINT.2I ' ' ~ 

2 FORMAT ( 'OFREQUENCY * ,5X,<AM PLITUDE' ,5X,* T-VALUE* >8X, * PH ASE * ,/ , * 0* I OUTP 

OUTP' 
OUTP 
OUTP 
OUTP 
OUTP 
OUTP 
OUTP*^ 
OUTP 
OUTP 
OU^’o 
ou • - 

OUTP 
OUT 'p' 


DO 55 LOANO=l ,N6ANnS 
55 WRITE (KPR1NT,3> B AN ( L BA NO ), B ( L BAND ), T ( L 8 ANCTI , PHASE ( L BAND ) 

3 format IF 10, 3,5X,Fq. 3, 5X,F7.3, 5X,F8. 31 

52 IF (JPLOT) 54.54,53 

53 CALL PLDT21 a , BM AX , 0. 0 , T , TMA X ,0. 0 ,NB ANOS , B AN, AMN AME ,TNAME , BANAME I 
CALL PLOT2(PH4SE,T0PT,BTnPI,T,TMAX,0.O, NBANOS,BAN,PHNAME , TNAME . 

I SANAME ) ' ^ ^ ^ 

54 CALL RFV( B, NBANOS I 

CALL RFV(T,NRAN0S) 

CALL REV(PHASE,NBAN0SI 
GO TO 100 

60 IF (JPRINT) 62,62,61 

61 WR 1TF“ (KPRINT,41 " 


10 

20 

30 

40 

50 

60' 

70 

80 

90 

100 

_110 

120 

130 

140 

150 

160 

iso 

190 
200 
210 
220 
_230 
240 
250 
260 
270 
280 
290 
300' 
310 
320 
330 
340 
350 
360 
370 
380 
39 0 
400 
410 
420 
430 
440 
45 0 
460 
470_ 
480' 
490 
500 
5tO 


530 

540 
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4 FORMATI «OPER lOD * , 5X. * AHPL 1 TUOE « >5X > «T-VALUE» ,8X> > PHASE* »/»* 0* lOUTP 550 


00 65 LBANO^l.NBANOS * 

65 WRITE (KPftINT,3» BANOSCLBANOI ♦8RBAN01 fTILBAMD) tPHASEILBAND) 

62 IF (JPLOT) 100*100*63 

63 CALL PtOT2 ( B* BMAX »0*0# T * TMAX *0*0*NBAN0S* BANDS* AMNAME* TNA ME • XNAME ) 

■ CALL PL0T2(PHASE,T0PI *BT0P1 •T,TMAX,0.0*NBAN0S*BANDS*PHNAME*TNAME. 

I XNAME) 

100 RETURN 

END — 


OUTP 

OUTP 

OUTP 

OUTP 

OUTP 

OUTP 


OUTP 

OUTP 


560 
570 
580 
590 
600 
61 0_ 
620 
630 





u u !ou u u 


SUBRDUT I NE Pi. QT2 C Y I . YMAX I *Y MI_N1 tV2» XMAxa , VH I N2 # NP# XAX f YNAM 1 » 

1 YNAH2,XNAM» 

DIMENSION Yim ,Y2( U,LINEmU>XAXaKXNAM<5),YNAMU9 >,YN^21U 
COHHQN KREAORtKPRNt 

DATA KBLNK/* • f *KBORD/* , V,KAL/«^V/fKPLpT/|*»/ 


_PLT2 

PLT2 

PLT2 

PLT2 

PLT2 


10 

20 

30 

40 

50 


PLT2 

PLT2 

PLT2 

PLT2 

PLT2 

PLT2 


C — PU0T2 GENERATES TWO OVEft/UNOER GRAPHS PER PAGE AND CALIBRATES 

“C— them. Yl-VERTICAL VARIABLE 1. YMAXlt YMINl*VERT ICAL PLOT LIMITS. 

C— Y2* YMAX2. YHI N2 = $t MI 1.AR FOR SECO N VARIABLE# NP«NUMBER OF POINTS^ 

C— XAX»x"AXI S (COMMON TO BOTH VAR I A B L E S ) . 

WRITE (KPRNT.2) YNAM 1 ,YNAM2 * XNAM 

2 FORMAT! *1* / T5.»*»/ T5 . , T25, • .TBO, •♦• / 

l*T80f •♦*♦*/ T3t **♦*♦♦• *T 18 • *♦♦♦*♦*♦?♦*** **9A4* T73. •♦♦♦♦♦♦*♦•♦♦* *PLT2 
^9 A 4/ T4.*****»T25. * ♦♦♦• « T80* •♦♦♦• /T5f •♦•tT25****»TBO# •♦•/IXt 5A4 IPLT 2 

PLT2_ 


PLT2 

PLT2 

T5,*** ,T25.»****PLT2 


60 
TO 
80 
90 
100 
110 - 
12 0 
130 
140 
150 
160 
170 


C — GENERATE SCALE FACTORS 

SKALl*50.5/( YMAXl-YMINU 
SKAL2*50.5/ (YMAX2-YMIN2I 
C — GENERATE CALIBRATION AND OUTPUT 
CAL 01 a YMINl 
CAL02 * YMIN2 


IT. 


PLT2 

PLT2 

PLT2 

PLT2 

PLT2 

PLT2 


CALI 

CAL2 


YMAX I 
YMAX2 


CALSlaYMAXl- (.5*( VMAXl-YMINl I \ 

CAL52=YMAX2' (. 5*( YMAX2-YMIN2n . 

WRITE (KPRNTtU CALO I .CAL 51 tCAL ltCAL02f C;AL52f CAL2 
1 FORMAT (• *pFl4 .4»F24. 4tF23 .4 .F 12 .4 t F23 .4»F24.4 t 
C — GENERATE VERT rCAt' AX I S AND OUTPUT IT 

DO 10 L « l.lll 


PLT2 

_PtT2 

PLT2 

PLT2 

PLT2 

PLT2 


10 


LINE(L» » 
LINE(I) a 


KBORD 
K At 

LINE! 26 » = KAL 
LINE(51) a KAL 


20 


LINE(61 ) a kaL 
LINE! 861 a KAL 
LINE( I llJa KAL 
DO 2D L=52*60 
LINE(L> a KBLNK 
WRITE! KPRNT.4 I LINE 


PLT2 

PLT2 

PLT2 

_PIT2 

PLT2 

PLT2 

PLT2 

J>LT2 

PLT2 

PLT2 

PLT2 

PLT2 


4 FORMAT! I OX. I HAD 


— DECIDE HOW 
NWIOE » 


MANY LINES 
32/NP 


BETWEEN POINTS 


— LOOP POR 
do" 30 


NUM BER 

LP a I 


OF P OINTS. 
, NP’ 


EACH TIME PLOTTING A POINT AND NWIDE 


—CLEAR LINE AND GENERATE GRID POINTS 

00 40 L = I . I 1 1 . . 

'40 LINE(L I = KBLNK 

LINE (ll=KftORO 


LINE! 26»=KBnRD 


PLT2 

PLT2 

PLT2 

PLT2 

PLT2 

$PAC EPLT2 

PLT2 

PLT2 

PLT2 

PLT? 

PLT2 

PLJ 2 

PLT2 


180 
190 
200 
210 
220 
219^ 
240 
250 
260 
2T0 
28 0 
29 0_ 
300 
310 
320 
330 
340 
35 0 _ 
360 
370 
380 
390 
400 
>10_ 
420 
430 
440 
450 
460 
470 

48 0 

49 0 
500 
510 
52 i 
530 
54 0 








xc 



LINE<5 IJ^K BOR D 

LI NE<6U=KR0RD 

__LINEI 861 =K BORO ^ 

LINEU I n = KSnRO * 

c __ ■ 

C — TEST FOR OVER OR UNDER' RANGE AND TRUNCATE OUT-OF-BOUNOS VALUES 
IF lYl (LPl -VHAXl J 2000. 2000. 2100 


PL T2 5S0 
PLT2 560 
PLT2 570 
PLT2 580 
PLT2 590 
PLT2 600 
PLT2 610 


2100 YA s YMAXl 

GO TO 5100 



PLT2 

PLT2 

620 
63 0 

2000 IF ( Yll LPI-YMINl )5000. 5200. 5200 
5000 YA = YMINl 



PLT2 

PLT2 

660 

650 

GO TO 5100 
5200 YA*Yl(LP) 



PLT2 

PLT2 

66 0 
670 

5100 IF IY2(LP)-YMAX2) 200, 200. 210 
210 YB = YMAX2 

• 


PLT2 

PLT2 

68 0 
690 

GO TO 510 

200 IF <Y2 ILP)“YMIN2I 500, 520, 520 



PLT2 

PLT2 

700 

710 

500 YB * YHIN2 

GO TO 510 



PLT2 

PLT2 

720 

730 

520 Y8*Y2(LP) 



PLT2 

760 

C 



PLT2 

750 

C — NOW GENERATE POINT INDEXES 



PLT2 

760 

510 Kl = (VA -YMIND^SKALHI .0 



PLT2 

7T0 

K2 = < Y8 -Vf1lN2)*$KAL2^6l.O 



PLT2 

780 

LINEIKl)^ KPLOT 



PLT2 

790 

LINE(K2)» KPLOT 



PLT2 

800 

C 



PLT2 

810 

C — OUTPUT HORIZONTAL SCALE VALUE AND LINE 



PLT2 

820 

WRITE (KPRNT ,3) XAX(LP),LINE 



PLT2 

830 

3 FORMAT (1X,F 8.6, IX, lUAl) 



PLT2 

860 

C 



PLT2 

650 


— CLEAR LINE AND GENERATE BLANK LINES FOR SPACE 
IF (NWIOE) 30,30,110 
HO DO 95 L = l,lll 
95 LINE! L I =KBLNK 

"line (i)=KBnRo 

LI NE( SI I =KBORO 
L I’NEI 26 ) = KttnKU 

LINE! 61) =KBORD _ _ . 

■ LI NE< 86)=K0ORO 

LINE! I 1 1 ) = KBORD 

'C" " 

C — OUTPUT PL A IN L INF 

00 90 LO0P = l .NW lO^E 

90 WRITE (KPRNT,4) LINE _ 

■ 30 CONTINUE 


C — GENERATE RIGHT BORDER AND OUTPUT 
00 50 L * I, 1 U 
50 LINEIL) = KflORO 
L INEI I ) = KAL 
LINE (26) - KAL 
_ LI NE( 5U = KAL 
LINE(6l) = KAL 
LINE(86) * KAL 
LINE( 111) = KAL 


PLT2 860 
PLT2 870 
PLT2 880 
PLT2 890 
PLT2 900 
PLT2 910 
PLT2 920 
PLT2 930 
>LT2 960 
PLT2 950 
PLT2 960 
^LT2 9T0_ 
PL“T2“980 
PLT2 990 
PLT21000 
PLT21010 
PLT21020 
PLT21030_ 
PL T2 1060 
PLT21050 
PLT21060 
PLT21070 
PLT21080 
JPLT2109 0 
PL T2 I 100 




T 

C 

C 

c 

c 

c 

~c 

c 

c‘ 

c 

c" 

c 


SUBROUTINE TOT AM(X,N,NPEAK. ISAVE* M* PHASE# 8. AHP 
DIMENSION X( 1) « I SAVE ID ♦PHAS6Cl)#BI U 


SUBROUTINE 


TO TAM 
SOIOS 


«r 

_ ♦ 

COMPUTES the RESULTANT AMPL ITOoE OF tHE SUM OF NPEAK * 

WHICH HAVE their AMPLITUDES AND PHASES OETREMINED BV THE ♦ 


TOTH 

TOTM 

TOTM 

TOTM 

XOTJl. 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 


TA^LLING PROGRAM., tOTAM SIMPLY SELECTS THE APPROPRIATE SINE ♦ 

AND COSINE WAVES FROM X, WEIGHTS THEM, SUMS THEM ♦ 

PUTES THE AMPLITUDE OF THE SUM. • 

WARNING this HILL OESTROY^ftOH I OF X 

CONTAINS THE OR IGINAr DEPENDENT VARIABLE. ♦ 




C 

^^OOP ThRU THE FOLLOW ING NPEAK tTmES TO CRE Af rtHE VlESULf ANT " H INX 

C 


DO 100 LP=1. NPEAK 
K« ISAVE(LP) 


TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

XQJJl 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 


10 

20 

30 

40 

^.Q_ 

60 

70 

80 

90 

100 

Xi^ 

120 

130 

140 

150 

160 

170 


KC*K*2 

KS»KC«-l 


THE FOLLOWING COMPUTES THE PROPER AMPLITUDES FOR THE COSINE.>?LP„SINjE 
COMPONENTS# BASED ON PHASE AND AMPLITUDE ESTIMATES. 


PC^l.O 

PC*PC/C TANIPHASEIKM^ 
PS»PC*TANIPH4SEIKn 


PC) 


^SELECT SINE AND COSIND VALUES FROM X A 

SELECT SINE AND COSINE VALUES FROM X# SUM AND PUT INTO Xt 1' 


100 


00 100 L*l.N 
CALL LOCIKC.L *J#M#N,0) 
CALL LOCIKS.L .1 #M,N,0) 
X(L)-X4J)*PC ^ Xdl^PS 


C-COHPUTE APLLITUDE,OF RESULTANT 

C „ 

sx»o.o 

00 10 L«1.N 

10 SX“SX4.X(L) 

SX*SX/FLOATIN) 

SX S=0 . O' 

00 20 -L = ltN 
20 SXS = SXS«^( XiL )-SX) 

AMP-SQRTI SXS/FLOAT(NM 

RETURN 

END 


TOTM 

JOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

TOTM 

_TOT^ 

f bt M 

TOTM 

TOTM 

TOTM 

TOTM 

tWm 

TOTM 
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APPENDIX 2 


Three Exanples o£ Use of SPECT 


NOTE: 


In Tfest 1, as in other tests, it will be noted that 
the final air 5 >litude spectrum has values around the maxima 
which are not zero. However, the t*spectra always have 
only two (or the correct number) of peaks. Amplitude spectra 
should be interpreted by considering the t ■'Values. 
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iHflTRix nF CORRtLATIONS 

VAPU'ME 1 IS TME dependent VARIABLE, ALL OTHERS 'ARE COSINE AND SINE WAVES OF VARIOUS PERIODS 
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VAAlAHt I IS tHE OEPENOENT VARIABLE. ALL OTHERS ARE COSfME AND SINE WAVES Of VARIOUS fERIOOS 
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32,000 

0.036 

1.014 

0.954 
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* • * * • • « PftOCAAN P*R4«fTEHS 

o' ■ 20 100 d I “o 



monte 


CARLO 


simulated data 


BAwn ? * ^*2?2?t FREQUEHCV* 0.0A50, AMRLTTunE* 

.2 r C 100* 25.0300t_ FRfOOENCY*. 0„.0^00i AMRL|Tljn£« 

.RANOON_CAijSSlAN NOl SE ANPUTUHE* 0.5000 


SDPR TIRE l:l2.6 





□ 


MATRIX OF CniRflATlONS 

VA»i i^'lE i is THE OEPENOENT VARUSLEt ALL OTHERS ARf COSINE AND SINE NAVES OF VARIOUS PERIODS 


iO 


u 


12 


11 


lA 


19 




1 

■ - I#- 
«■/> 

i.oro 

-0.1*2 l.OOO 
-n.^22 0.912 

-0, 1 79 0.936 

-0.0*6 

-n.237 O.fSR 

1.000 

0.320 1.' 

0.919 0. 

0.565 0. 

7 

-O.o*') -0.S25 

0.761 -0, 

R 

-0.316 0.il9 

0.709 0, 

9 

-0.967 -C.670 

0.504 -0. 

10 

-0.603 0.262 

0,739 0. 

ll 

-'l.C'n -0.697 

0.2?5 -C. 

12 

-0.6HS 0.032 

0.£,55 C. 

15 

Q.0H3 -0.622 

-3.01? -C. 

lA 

-0.S36 -0.129 

0.461 0. 

IS 

0.207 -0.692 

-0. 167 -0. 

16 

-0.S<.3 -0.190 

-0, 

17 

0. H)2 -J, 310 

-0.227 -0. 

18 

-0.501 -0.106 

0.162 -0, 

19 

0.525 -0.169 

-6.194 -P, 

20 

-'^,60a -3.U2 

0.062 -0, 

21 

3,6 7C -0.035 

-0.102 -C, 

22 

-0.262 -0,013 

0.P13 -C, 

23 

0,776 0.012 

0,005 -0. 

7 *> 

-0.078 0,076 

0.066 -0. 

25 

0,8 S3 -O.OCl 

0.090 0. 

26 

0.M5 0.127 

0,114 0. 

27 

0,850 “0.057 

0, 129 0. 

28 

0.216 0.131 

0.183 0. 

29 

0. 797 -0. 125 

0.il3 -0. 

30 

0,662 0.095 

0.229 0. 

91 

0.6 'j 6* -3.177 

0.057 -0. 

32 

0,962 0,036 

0.233 0. 

33 

0.513 -'J.195 

-O.OU -0. 

36 

0.577 -0.023 

0.193 0, 

39 

0, 360 -O.l 79 

-0,072 -0. 

36 

0,557 -0.069 

0.129 -0. 

37 

0.218 -0,133 

-6.101 -0, 

38 

0.503 -0.079 

0.067 -C. 

39 

O.lOO -0.072 

-0.090 -0, 

60 

0,614 -0.058 

0.024 -0, 

41 

0.021 -O.OU 

-0.049 -0 

16 

17 

Ift 

16 1 7 

l.OCO 

0^029 l.COO 
0,961 0.320 

Id r<; 

1,000 

19 

-0.282 0.9 i2 

0.008 1 

20 

0,774 0,550 

0.919 n 

21 

-0.533 0.746 

-0.298 0 

2 > 

0,523 0,701 

0.769 0 

25 

-0,678 0.499 

-0.537 0 

24 

0.244 0.725 

0.499 0 



l.OOO 

►0»00 l 

0.9U 

•0.299 

0.7A8 


1,000 

0.10T 

0,959 

0,960 


1,000 

0.015 

0,955 


1,000 

0.526 


1.000 



xxvl 



?*> 

- 3. 7 T9 

0,?9^ -n,#>77 

0.923 -3.530 

0.777 

-P. 2 79 

0,09 7 

0,0 3.? 

1 ,000 







2»- 

9 . f. r ) 

0,99 1 0, ? />#. 

r>. 7 H 0.406 

0. 710 

0. 756 

0.5 69 

0. 01 1 

0. 132 

1 .000 






? r 

-’.),<> JO 

O.Ol i -0. 700 

0.760 -0.690 

0.527 

-0.520 

0. 76 5 

-P. 7 7(3 

C.9 35 

0. 029 

1 .000 





?fl 

» ‘•f 

0,9OR 0.090 

0.969 (3.2 30 

0, 795 

0,519 

0.719 

0. 7 72 

0.5M 

0.991 

0. 320 

l.OOO 



s 

7«3 

-0.967 

-O.l W -0.600 

0.020 -0.69S 

0.298 

-0.666 

0. 5(35 

-0. S?2 

0. 799 

-0.282 

0.932 

0,000 

l.OOO 


10_ 

-n.? ta 

_p. 339 -0. 19 7 

P. 979 O.P’O 

0.6 7 6 

P.27’ 5 

0. 795 

0,537 

. 0.71(3 

0. 7 75 

C.559 

0.9)9 

0, 305 

1,000 


31 

-0,i72 

-3, >02 -0.993 

-0.19 7 '-0. 9 OP 

-P.OPO 

-0.691 

0,2 22 

-0. 672 

0,907 

-0.5 32 

0.797 

-0,297 

0.939 

0.000 


1? 

-0. laa 

o.ia* -0.200 

0.350 -0.129 

0.529 

0.036 

0.963 

0.769 

0.739 

0.525 

0.7 05 

0. 760 

0.556 

0.9 19 


33 

-0, l 09 

-0,1 -3.2M 

-0.723 -3.99<i 

-0. 198 

-0,61 7 

-13.016 

-0.7P9 

0.219 

-0.6 79 

0,500 

-0, 536 

0. 763 

-0.292 


I** 

-0.399 

0,079 -0.169 

C. 190 -0. 195 

«.3 30 

-0. 1 29 

0.996 

0.075 

C.693 

0,299 

0.720 

0,999 

0,706 

0. 79 7 


Ifi 

0.0C2 

-0. 113 -0. no 

-C.2I3 -0.283 

-0.290 

-0,990 

-0. 1 72 

-0,690 

-0.OO5 

-0. 7 13 

0.299 

-0.677 

0.527 

-P.526 


“ih 

0*0 la 

0.339 -o.ora 

0,061 -0. 162 

0. 1 55 

-0. 195 

0. 309 

--0.19 3 

0.9 82 

, 0.002 

,0.695 

.0,225 

0.737 

0.999 


37 

0,039' 

-0.02 7 -o', c 1(9 

-0.130 -0.137 

-0.220 

-0. 31?' 

-0.23? 

-0, 996 

-0,155 

-0.6 39 

0.019 

-0.709 

0.265 

-0.660 


3R 

0. K7 

0.096 0.C21 

P.'’U -0.099 

0.090 

-0, 1»9 

0. 192 

-0,217 

0.307 

-0.158 

o.soi 

-0.002 

0.66a 

0.23 7 


39 

-o,oci 

0.35S O.OlO 

-C.033 -0. 038 

-O.IJO 

“0. 151 

-0.199 

-0, 309 

-0.210 

-0.976 

-0,139 

-0,615 

0.023 

-0,690 


<»n 

0. 1 ■jO 

0.005 0.007 

0.029 -0,000 

-0.001 

-0.115 

0.092 

-0,203 

0,160 

-0.221 

0.337 

-0. 199 

0,529 

0.020 


9t 

-0,0 79 

O.lOl -0.022 

0,063 -C.002 

-0,013 

-0.036 

-0.109 

♦0,13} 

-0,160 

■0.276 

-0.205 

-0,996 

-0. 197 

-0.600 



31 

32 3 3 

39 15 

36 

37 

1« 

39 

6C 

91 






i.ccc 

J? O.JIO l.OOO 


33 

0,990 

0.019 l.OOO 










39 

0.565 

0,932 0.3)2 

1,000 









15 

0.776 

-0.»70 0,991. 

0.093 1.000 









36 

0. ru 

0. 755 0. 576 

0.933 0. 392 

1.000 








37 

0.532 

-0.519 0, 767 

-0.259 C.936 

0.090 l.OOO 








30 

0.798 

0.519 0.718 

0,779 0,575 

0,99? 0,327 

l.OOO 







19 

0.251 

-0,662 0,505 

-0.515 0.797 

-0.275 0.930 

-JbQU. 

_.1,.000_ 






90 

0.6 78 

0.267 0,796 

0.590 0.712 

0.778 0.562 

0.990 

0.306 

l.OOO 





91 

O.PCO 

-0.689 0.220 

-0.667 0.982 

-0.520 0.79? 

-0.296 

0,933 

-0.002 l.OOO 





S t N 

0 L F 

8 A N 0 

S P E C T « U 

H 





FRro»ff9Cr 

AHPL ITUOF 

T-VAtUE 

PHASE 




















0. 

035 

0.958 

9.770 

1.523 






, 


0, 

016 

i.ni 

5.652 

1.379 








0. 

0 37 

1.318 

6.995 

1.198 








0. 

0)0 

1 .999 

7.972 

1.019 








0. 

039 

3.699 

8,699 

o.e 1 ) 








0. 

09C 

1.758 

10.351 

0.508 








0. 

091 

1 .651 

12.970 

0,157 








P. 

09? 

1.8)1 

16.019 

'7,135 








0. 

P93 

1 .866 

15.533 

-(0.092 








0, 

099 

1.810 

12.567 

-0.125 








0. 

095 

1.729 

9,966 

-0.597 








0. 

096 

1.601 

8.383 

-0,762 








0. 

P97 

1.9)3 

7.96S 

-0.983 








0. 

09a 

1.26 7 

6,6 30 

-1,70? 








0. 

099 

1 ,065 

5,552 

-1.901 








0. 

PSP 

0.880 

9.161 

1,567 








0. 

051 

0.700 

). )10 

1.921 








0. 

052 

0.518 

2,983 

1.328 








0. 

053 

0.907 

1.859 

1.32) 








0. 

059 

0.31 7 

1.938 

1.619 
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FPfOUENCy 

f >.0 




0.0 

0.0340 

0.0170 


C.C1SQ 

0.0370 

0.0400 


0,0410 
0,C4?0 
0,0 MO 


0,'^440 
0, 04*50 
0.C440 


0»Q^3^ 

0.0490 


C..0500 

0,0910 

0,05^0 


O.OSIO 

0.0540 


AHPUTune FSIHATE* 


0.9454 


* • « 

> •• •« 


r-VALoe, 


1.6907 0.„0_ 


__9.0095 


*16,0190 
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♦ . 
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XXIX 


□ 

I 


FPt. ino/FPr c. 


cn»«fLAfpjN» 

ITiJDf: 

1. B6&0 


r>.7?lT 

T-VAIUF 


Of = 0.'»'?9933779Z3<J655f 00 

NlfLT|DLF ^CUARFO CORBr,L*TinN« 0,701*5 
PEfinO/FREO, AHPLITUOe T-VALUe 

O.OAaC u 6 SI0 12,9776 


4 



co“»*ijTrir, !.PtCTtiijH u«:iNr, T‘4F fniirwiNr, in mqdfi 

0, T^TC 


detfrhinakt • 


0. 151 12/2^l55lfi76F-0l 


HlJlTIPl.F 5CtlFRE0 CORRELATION* 0.887J 

JfU ^1'-^ 0* A>* Pl 1 TOOF T- VAHJE 

■"b.0<.60' 5^0615 13.<J587 

O.rXiRO 25.15«)8 



OFTf aNlNa Nt « 0 . 67 7 1 7 1 3»<?5 25452F- 0 1 


muitIPLF tclaheo 

PfR lOn/FOTO, 
0.04Bf 


CnP^ElATirJN* 
ampl 1 Tunp 
2.57J2 
6,9'%6 7 


O.0'7OO 

t-value 

26,9'726 
25, 1598 


' 6fTFP><JN46T « 0^179735681 739066E 00 

kUITIplF SCI AREO C0RRELAM0N= 0.8893 

PF9 ton /r9 f, P_, AMP L I T JOE T. VALUE 

b^C*4V0 1.765 7 19,2663 

0.0A8C 6,9667 25,1593 


2 PFTFRHINANT ». 0 , 1 53 5297 1 85 I 8 25 7E‘0 1 


multiple SClARED correlation* 0,8882 
PEPino/CRFO, AMPLITUOr ” T-VAlUE 

0,065r 2,5702 26,9926 

' b,b67C 1,8832 ■ ‘ 36.9735 


DFTERMtNANI • 0. 1 77 635967 7 3 1 6 76E 00 

multiple SCL'ARtO CnRRFLATION= 0.8893 
PEP imvFREc. amplitude T-VALUE 

0,065C‘ 2,5702 26.9926 

0.0690 2,6569 22,6255 











D 




« 


FUFOtJF.^CY 

- 6,2832 


PM4SE ESTIMftTf, 


0.0 


6.2832 


0,0^60 
o.o^sn 
n.03 70 


O.OliO 

0.0390 

0.0400 


0.0410 

0.0420 

0.0430 


,0.0440 

0.0450 

0.0460 


0, 0470 
_0._q4_80 
0.0490 


,o-q‘'QP 

0.0510 

0.0520 


0.0530 

0.0540 




I*ULT|PLE 5 A^ARED C0R8ELATinN« 0.6903 


T-VALUE 


t 


**• 

* 


0.0 


1 3.9 109 


JT,6618 









